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ABSTRACT 



EXPLORING THE DIFFERENT PHASE DIAGRAMS 
OF STRONG INTERACTIONS 

Leticia Faria Domingues Palhares 
Advisor: Eduardo Souza Fraga 

Abstract of the PhD thesis presented to the Graduate Program in Physics of the Physics 
Institute of the Federal University of Rio de Janeiro, as part of the requirements for obtaining 
the title of Doctor in Sciences (Physics). 

In-medium field-theory is applied to different effective models and QCD to describe mass 
and isospin effects, finite volume corrections and magnetic fields in the phase diagram of 
Strong Interactions, keeping close contact with experiments and lattice results. Findings 
range from a technical nonperturbative solution of the hot and dense regime of a general 
massive Yukawa Theory and the computation of the cold and dense chiral surface tension 
- a key quantity for supernovae explosions and compact star structure - to the proposal of 
a novel signature for the QCD critical endpoint in heavy-ion collisions based on finite-size 
scaling. The behavior of the deconfining critical temperature as a function of the pion mass 
and the isospin chemical potential is also addressed in an effective model and the description 
obtained is in agreement with lattice simulations, in contrast to what is found in (Polyakov- 
extended) chiral models. We also discuss the thermodynamics of QCD in the presence of a 
(Abelian) magnetic field to two-loop order and a Functional- Renormalization- Group analysis 
of relativistic Bose- Einstein condensation of pions in isospin-dense media. 

Key-words: Quantum field theory at finite temperature and density. Effective Models, 
Phase Transitions in Quantum Chromo dynamics 
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Chapter 1 




Introduction: the genesis of 
the different phase diagrams 



A broad comprehension of the physical phenomena governed by a given fundamental the- 
ory requires the description of the phase structure it predicts. Since the statistical-mechanics 
definition of the phase of a system is intimately connected to symmetry characteristics, the 
first step towards understanding the phase structure of any system consists therefore of an 
investigation of the symmetries that play an important role in its characterization. In the 
case of Strong Interactions, the corresponding fundamental theory. Quantum Chromody- 
namics (QCD), presents a very rich structure of symmetries in the vacuum [1]. In this vein, 
an intricate structure of phases is also expected. 

The different phases of Strong Interactions and the transitions between them play im- 
portant roles in different physical scenarios, with observational and experimental interest. 
In the primordial universe, for example, it is generally believed that there existed an era 
dominated by a hot soup of quarks and gluons at about 10~^^ seconds after the Big Bang 
[2] . As the universe expanded and cooled down, a primordial confinement phase transition is 
thus expected. Our universe today may also contain naturally occurring deconfined systems: 
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estimates point out that the core of compact stars could accommodate high enough pres- 
sure so that stable matter is not the conventional nuclear one, nor the hadronic dominated 
medium, but rather a system of deconfined quarks and gluons [3]. Moreover, the energy 
scales of the phase transitions of Strong Interactions are within the reach of laboratory ex- 
periments colliding ultra- relativistic heavy ions [4], such as RHIC-BNL [5] and LHC-CERN 
[6] and, in the near future, FAIR-GSI [7] and NICA-JINR [8]. 

Several studies and developments were accomplished in the last decades concerning the 
investigation of the phase structure of Strong Interactions. In special, the description of 
the phase diagram in the plane of temperature T versus baryon chemical potential fiB has 
received remarkable attention, with sensibly deeper advances on the temperature axis, due 
to the development of lattice QCD in this domain [9]. Experimentally, heavy-ion collisions at 
RHIC-BNL have established the existence with a new state of matter dominated by partonic 
degrees of freedom [10]. 

Quantitatively, however, the knowledge about the phase diagram of Strong Interactions 
is still extremely incipient. There are even some qualitative open questions, such as e.g. the 
existence of the second-order critical endpoint (CEP) as the final point of a presumed first- 
order transition line and what concerns the connection between chiral symmetry breaking 
and confinement, besides the structure of the region where color superconductivity should 
occur [11]. The phase structure of Strong Interactions is essentially, therefore, an open 
problem. 

This is due to various difficulties encountered in the description of Strongly interacting 
matter under extreme conditions. First of all, it is clear that the investigation of any phase 
diagram requires the consideration of the theory in a broad domain of energy scales. In 
quantum field theory, according to the renormalization group formalism, the study of the 
predictions of a theory for different energy scales must be implemented via running param- 
eters (couplings, masses, etc). In this vein, when large enough variations of the energies are 
considered, non-perturbative regimes arise naturally. 

Among the non-perturbative methods to approach the fundamental gauge theory in this 
regime, lattice QCD stands out, since it solves in principle the full fundamental theory. 
There are, nevertheless, technical obstacles even for simulations of QCD on a lattice. The 
limitation of maximum lattice sizes and minimum lattice spacings imposed by reasonable 
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computing times generates infra-red and ultra-violet cutoffs for the physical phenomena 
being analyzed. A deeper restriction, the notorious Sign Problem [9] renders impossible 
the utilization of conventional Monte Carlo methods (based on importance sampling) in 
the case with a finite baryon chemical potential, in which the weight of the path integral 
representation of the partition function becomes complex. Summing up, there exists no single 
method or theory that can be successfully applied to obtain the complete phase structure of 
Strong Interactions. The theoretical development in this area of research relies necessarily on 
the exploration of the complementarity of the various available methods, including different 
effective treatments. 

On the experimental side, the years that follow promise significant advances in this area. 
Different experiments and programs are planned to scrutinize a broad domain of intermediate 
chemical potentials, with the general goal of placing the first experimental point in the QCD 
phase diagram at high energies. Searches concentrate not only on critical correlations built 
up in a second-order phase transition (at the CEP), but also on possible signatures of the 
expected first-order line at higher baryon-antibaryon asymmetries. From the observational 
side, increasing improvements in astronomy recently point out to the possibility that precise 
data from astrophysical phenomena will soon provide constraints to ultra-dense Strongly 
interacting media. This experimental and observational enterprise will certainly contribute 
concretely to the knowledge of how the predicted critical region of the phase diagram of 
Strong Interactions is in fact realized in nature. 

The conjuncture of factors presented above implies an intense demand for theoretical 
studies of the phase diagram of Strong Interactions. During this challenging odyssey to 
overcome the obstacles in the description of the phase structure of Strongly interacting hot 
and dense matter, new control parameters arise naturally, springing new axes and new phase 
diagrams. The genesis of the different phase diagrams of Strong Interactions is therefore 
spontaneously triggered by the intricacy of the task of understanding Strongly interacting 
environments. 

Axes involving the numbers of colors and flavors of the theory or the masses of the 
fundamental fields are examples of phase diagrams that are built due to theoretical difficulties 
that drive the study of "cousin" theories forward. Even though variations of these parameters 
are not realized in nature, one gains in comprehension of their role in the predictions via the 
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study of theories which are not the original, physical one but "next of kin" to it. Mass effects 
are investigated in this thesis nonperturbatively in the context of a general Yukawa theory 
with massive fermions and bosons at finite temperature and density. This is a general result, 
with potential application in effective models of QCD and in different phenomena concerning 
media with interacting fermions. 

On the other hand, computational technical difficulties may also motivate the analyses 
of new control parameters: in order to circumvent the Sign Problem when simulating dense 
media on the lattice, one may study imaginary chemical potentials and the isospin chemical 
potential, for example. The latter should also occur in nature: the imposition of charge 
neutrality of compact astrophysical objects generates an asymmetry between quarks up and 
down in the ultra-dense core. We construct in this thesis an effective model inspired in chiral 
perturbation theory and perturbative QCD with the aim of investigating variations of masses 
and baryon and isospin chemical potentials. Findings agree with lattice QCD simulations 
at finite temperature and involve further predictions within the regime in which the Sign 
Problem is manifest, as the critical baryon chemical potential at low temperatures, which 
may affect compact star structure. 

Moreover, extra macroscopic parameters may arise from its relevance in a specific exper- 
imental setup. In heavy-ion collisions, the energy scales associated with typical sizes of the 
system created are of the same order as the energies within the critical region of the phase 
transitions, so that it is important to investigate finite-size effects. Indeed, we show using the 
linear sigma model with constituent quarks that the pseudo-critical phase diagram that is 
currently being probed in heavy-ion collisions may differ considerably from the expectations 
in the thermodynamic limit. We further propose a finite-size scaling analysis as an alterna- 
tive signature for the second order critical endpoint of QCD and, in collaboration with an 
experimentalist from STAR at the accelerator RHIC-BNL, demonstrate the viability of the 
application of this method to data analysis. 

In Nature and in the laboratory, the QCD phase transitions take place generally out of 
equilibrium. The phase structure probed is therefore that of dynamical critical phenomena. 
In this direction, we compute the chiral surface tension in the cold and dense regime within 
the linear sigma model, suggesting a much lower value as compared to previous estimates, 
possibly allowing for observable signals of a QCD phase transition in core-collapse supernova 
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explosions. Still in an out-of-equilibrium context, memory effects could be relevant to the 
time evolution of an order parameter (as can be seen in simpler models in quantum statistical 
mechanics [12]). 

The initial stage of noncentral collisions and the core of some ultra-compact astrophysical 
objects (the so called "magnetars") feature Strongly interacting media exposed to enormous 
magnetic fields, pointing thus to the study of background magnetic fields as a new axis 
of the phase diagram of Strong Interactions. We address this problem by analyzing QCD 
thermodynamics in the presence of a uniform and constant intense magnetic field. Questions 
like "what are the adequate quasiparticles?" and "how does the non-Abelian interaction 
affects the pressure in this extreme regime?" are investigated analytically within the well- 
defined context of the fundamental gauge theory. Results can be compared to recent lattice 
QCD simulations, and one expects in this way to gain information about different ingredients 
of in-medium QCD and about the behavior of perturbation theory in this domain. 

With the progress of the analyses of these new phase diagrams - initially considered as 
side effects of the original quest of predicting the temperature-baryon chemical potential 
phase structure -, it is not surprising that some of them develop independent relevance. At 
least academically, they are in principle part of the phase structure of the Standard Model of 
fundamental interactions. Besides that, such conditions may be very interesting phenomeno- 
logically, since new axes in general allow for the existence of new phases. The standard case is 
that of the baryon chemical potential, which gives rise to color superconducting phases [13]. 
Nevertheless, many other examples appear in the literature. The existence of a new confined 
phase, "quarkyonic" , has been proposed based on its formal prediction for the theory with 
infinite number of colors [14]. More recently, the consideration of the magnetic field axis 
has opened the discussion on the possibility of a superconducting lowest-energy-state due 
to p— meson condensation at intense magnetic backgrounds [15]. Other chemical potentials 
(isospin/strangeness) should also generate new phases, corresponding to Bose-Einstein con- 
densation of pions/kaons. In the final part of this thesis, pion Bose-Einstein condensation 
in isospin dense media is discussed. The effect of interactions has proven to be a non-trivial 
issue in non-relativistic Bose-Einstein condensation. Using the nonperturbative framework 
of the Functional renormalization Group, we address this question in the relativistic case, 
with the application for Strongly interacting media in mind. 
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Our goal in this thesis is therefore an exploration of some of these different phase diagrams 
through the inescapable utilization of various techniques, formal and phenomenological, in- 
volving effective models and the fundamental gauge theory, with a close connection with 
experiments and lattice simulations whenever possible. We hope that this study contributes 
to understanding the role played by different ingredients encountered in Strongly interacting 
environments and ultimately to the crucial and arduous task of showing that the underlying 
microscopic theory of Strong Interactions and our interpretation of it are consistent with 
nature. 

The thesis is organized as follows. After this initial introductory chapter, results on dif- 
ferent phase diagrams of Strong Interactions are exposed in Chapters 2-7. Mass effects are 
analyzed for the deconfinement phase transition in an isospin-dense QCD medium in Chap- 
ter 2. The pseudo-critical phase diagram of Strong Interactions, including the typical system 
size as an extra axis, is presented in Chapter 3. Consequences for observables in heavy-ion 
experiments, especially those related to the critical endpoint, are investigated. Chapter 4 ad- 
dresses how out-of-equilibrium phenomena may influence the QCD phase transitions actually 
probed in observations and experiments: nucleation parameters are computed for the cold 
and dense chiral phase transition, with applications in astrophysics. The thermodynamics 
of QCD in a magnetic background is analyzed in Chapter 5. Nonperturbative corrections 
to effective models are addressed by computing the pressure of a general massive Yukawa 
theory within Optimized Perturbation Theory at finite temperature and baryon chemical 
potential in Chapter 6. Chapter 7 contains a Functional Renormalization Group analysis of 
relativistic Bose-Einstein condensation in an interacting, charged scalar field theory, as well 
as in a chiral model for an isospin-asymmetric pion medium. Finally, the conclusions of this 
thesis are presented in Chapter 8. Appendices contain some technical issues tackled during 
the development of this work. 



Chapter 2 




Phase diagrams of Strong 
Interactions in the presence of 
various chemical potentials and masses 



Net baryon number is commonly considered when investigating Strongly interacting mat- 
ter. Nevertheless, there are other important conserved quantum numbers involving quark 
degrees of freedom, such as isospin and strangeness, which are usually set to zero. In this 
chapter we will analyze a more general type of QCD medium including the possibility of a 
nonzero net isospin density and investigate its interplay with finite quark mass effects. 

The specific problem we address is that of describing the deconfinement critical temper- 
ature as a multiple variable function, Tc{fiB, f^iif^Tr) [16], and its analogous at zero temper- 
ature, the critical baryon chemical potential fi'^^^^{fij,mT^) [17]. Due to the absence of the 
Sign Problem in isospin dense QCD, comparison of Tc{fiB = 0,/i/,m7r) with lattice results 
is possible, providing information on the ingredients in an effective theory that suffice to 
capture the role played by mass and isospin effects in the critical region. The model we 
construct using this criterium supplies predictions also within the Sign Problem domain; 
the zero temperature critical baryon chemical potential being relevant to the astrophysics of 
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compact objects. 

2.1 Introduction 

Small quark masses and a nonzero baryon chemical potential have always represented major 
challenges for lattice simulations. Presently, although viable lattice sizes still prevent ex- 
tensive and precise studies with realistic quark masses, and the Sign Problem considerably 
restricts the applicability of Monte Carlo simulations to the description of chemically asym- 
metric media. Lattice QCD is starting to provide results with smaller quark masses, and 
probing, with approximate methods, a larger domain at finite chemical potential [18, 19, 20]. 
For the latter, it is believed that much can be learned from simulations of realizations of 
QCD that avoid the Sign Problem, such as those with vanishing baryon chemical potential 
and finite isospin density, which has a positive fermionic determinant. 

QCD at finite isospin density is certainly, but not only, a playground to test numerical 
approaches to the case of finite baryon density on the Lattice [21, 22]. It is also part of the 
physical phase diagram for strong interactions, and exhibits a very rich phenomenology [23]. 
It has been under careful study during the last years, both theoretically and experimentally, 
with a clear identification of certain phenomena that depend directly on the isospin asym- 
metry of nuclear matter at intermediate-energy heavy ion experiments (see Ref. [24] and 
references therein). 

Nevertheless, theoretical and phenomenological studies often focus on the chiral limit of 
QCD, putting aside effects from finite quark masses, and on isospin symmetric hot mat- 
ter, mainly stimulated by the physical scenario found in current high-energy heavy ion col- 
lision experiments [25] and the quark-gluon plasma [11]. Some exceptions are, however, 
chiral model analyses of pion-mass dependence of the finite temperature transition (e.g., 
[26, 27, 28]), which fail completely to reproduce the lattice behavior, and results within 
Polyakov-loop models that are consistent with lattice via fitting procedures. At finite isospin 
chemical potential, though there are phase-diagram investigations in PNJL models [29], none 
of them has addressed the critical temperature dependence on the isospin chemical potential 
in particular. 

To investigate the effects of nonzero quark masses and isospin asymmetry on the de- 
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confining transition, we build an effective tlieory that combines ingredients from xPT in 
tfie low-energy sector witli tlie plienomenological fuzzy bag model at high energy. The 
high-energy regime is described perturbatively by two-flavor QCD with massive quarks and 
explicit isospin symmetry breaking. Nonperturbative (confinement) effects at this scale are 
incorporated through a fuzzy bag description [30] with coefficients extracted from lattice 
data. For the low-energy sector, we adopt an effective action inspired by xPT that exhibits 
the same structure of symmetries contained in the high-energy theory [31]. The quasi-nucleon 
degrees of freedom described in this regime seem to be crucial for understanding the pion- 
mass dependence of the deconfining temperature. The definition of parameters, as well as 
masses, is such that variations in the deconfined sector are totally consistent with variations 
in the confined one, which guarantees the bookkeeping in the different degrees of freedom 
that are present in the description. 

In this chapter of the thesis, after presenting the effective model, we use it to investigate 
the effects of finite quark masses and isospin number on the equation of state of hot and dense 
strongly interacting matter and on the deconfining phase transition. The setting we have 
proposed is simple and completely fixed by vacuum QCD properties (measured or simulated 
on the lattice) and lattice simulations of finite temperature QCD. More explicitly, there is no 
fitting of mass or isospin chemical potential dependence at all. Nevertheless, our findings for 
the behavior of the critical temperature as a function of both the pion mass and the isospin 
chemical potential are in remarkably good agreement with lattice data. 

It is crucial to note that several detailed studies of chiral models failed to describe Tc{mT^) 
[26, 27, 28], while Polyakov-loop models, whose results can be fitted to the lattice points for 
Tc(m7r) [27], cannot at the present form address isospin effects (nor any other observable 
related to fermionic chemical potentials). This suggests that the functional behavior of the 
critical temperature is a good observable of effects from mass variations, being sensitive to 
the physics entering in each case. The full function Tc (m^) is probably much more sensi- 
tive to the underlying physical phenomena than the actual value of the critical temperature 
in Nature. Putting together the fact that first-principle lattice QCD simulations can access 
these "cousin" theories, with different masses, a thorough analysis of the QCD critical temper- 
ature as a multi- variable observable might be a path worth following to identify the possible 
differences between chiral and deconfinement transitions and understand why they seem to 
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be connected in real- world QCD. These issues are also related to the question of whether 
QCD is more influenced by the pure gauge/large physics or by its chiral properties, the 
top and bottom regions of the cartoon phase diagram in Figure 2.1. 




T 

Figure 2.1: Cartoon of the QCD mass phase diagram. 

In our approach, one can also investigate in a straightforward manner the effects of quark- 
mass asymmetry and nonzero baryon chemical potential, physical cases in which the Sign 
Problem develops, constraining systematic lattice studies. The case of cold, dense matter is, 
of course, relevant for compact stars [32] and for the phenomenon of color superconductivity 
[13]. 

The chapter is organized as follows. Section 2.2 presents the effective model adopted, as 
well as a discussion of the approximations implied. In Section 2.3 we show our results for 
the effects of mass and isospin chemical potential on the critical temperature and critical 
baryon chemical potential, comparing to lattice data whenever possible. We also compare 
some of these results to those obtained in different phenomenological frameworks. Section 
2.4 contains final remarks. Some details of the technical derivation of the pion physical 
(dressed) masses are left for Appendix D. 
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2.2 Effective theory 

As discussed above, to investigate the effects of nonzero quark masses, baryon chemical 
potential, and isospin asymmetry on the deconfining transition, we build an effective theory 
that combines ingredients from chiral perturbation theory in the low-energy sector with the 
phenomenological fuzzy bag model at high energy. In this section, the model and its degrees 
of freedom will be presented explicitly, showing thus in a concrete setting the consistent 
connection between the parameters and masses of the two sectors. 



2.2.1 Low-energy sector 

The physical setting we adopt for the low-energy regime of strong interactions is that of a 
system of heavy nucleons in the presence of a hot gas of pions whose masses are already 
dressed by corrections from temperature, isospin chemical potential and quark masses. The 
effective Lagrangian has the form Ceff = + C„n + C-,,, with 

Cn = N (^t^- MN + ^fiHoT^ + ^fiBlo^ N , (2.1) 

9a— ,^ f ^_ 1 ... r 3 



C^N = yiiVz75 (^^v^--^,7o[r^7r]J AT, (2.2) 

= {d^ - ih5^^)Ti+{d, + ih5,,)Ti~ - m^TT+TT- + ^ [{dTiy - m^(7r°)2] , (2.3) 

where qa is the axial vector current coefficient of the nucleon, which accounts for renormal- 
ization in the weak decay rate of the neutron, /^r is the pion decay constant, and /i is a 
function of temperature and isospin chemical potential. Here the nucleons are represented 
by = (p, n) with p, n being the proton and neutron spinors, respectively, and have a mass 
matrix Mn = diag(Mp, M^) = diag(M - 6M, M + 6M). This is the 0{P) nucleon chiral La- 
grangian [33], but considering mass corrections at zero temperature and chemical potential, 
and coupling to the dressed pions. 

The effective (dressed) masses of the pions 7r° = and vr^ = ^(tt^ =i=i7r^), which depend 
on temperature T, isospin chemical potential fij and mass asymmetry 6m = [ma — m„)/2, 
are denoted, respectively, by mg = fnQ{T, fij,6m) and fn± = fn±{T, fij). The isospin and 
baryonic chemical potentials can be written in terms of the chemical potentials for the up 
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and down quarks. Here we adopt the following convenient definition^: fij = fiu — f^d and 
fJ'B = fJ'u + fJ'd- We also use the customary notation tt = vr^r'*, with r" being the Pauli 
matrices. 

The physical masses in the dressed theory defined above were computed in Refs. [31]. 
In the so-called first phase, a regime in which \fij\ < m-jr, they have the following form (see 
Appendix D for more details): 



mjrO = mo = m. 
m^ri = m± ^ h = m 



1 

1 + 2"^oo 



1 1 (To 

1 H — acTi ± -a — 

2 2 



(2.4) 

T/i/, (2.5) 



up to first order in a = {mT./AiifT.y. The functions (Too, cto; o"i and h are defined in the 
appendix. 

In the second phase {\fii\ > rriT^), a condensation of pions occurs, and a superfluid phase 
sets in^. In this new phase, in order to reestablish the vacuum structure, a chiral rotation is 
produced due to the isospin symmetry breaking. All this produces a pion mixing, and the 
nucleons also couple in a different way. The degrees of freedom do not correspond anymore 
to pions, but we can still call them quasi-pions since their masses in the two phases match at 
the transition point. The tree level masses do not have the shape as in the equations above. 
Instead, m^o = |yU,/|, m^- = 0, and = + "iivTiT^ j [iiY [23]. However, it is possible 

to treat this phase with a simple approximation near the superfluid phase transition. In the 
regime in which > m^, the natural expansion parameter is given by = 1 — m^/fij 
[31], after scaling all the parameters by The result of the first terms in this expansion 
(s^ = 0) provides the same equations as in the normal phase, just replacing by 



Strictly speaking, this is valid only for values of very close to m^, i.e. < a/8/7 772,^ 
(for a more detailed discussion cf. [31]). For simplicity, we also apply this results for slightly 
higher values of the isospin chemical potential^. These results in xPT are confirmed by a 
NJL analysis [34]. 



^The baryonic chemical potential is also frequently defined with an overall factor 2/3. 
^Here we consider the tt~ condensation taking ~ — 

■^For higher values of |^/| one needs to consider more terms, not only in the expansion of the pion 
Lagrangian but also in the coupling with the nucleons, due to the appearance of the pion condensate. These 
terms will also be present at the perturbative QCD level via the chiral rotation. The case in which 3> TOtt 
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The direct effect of the baryonic chemical potential in the pure pion quasiparticle gas is 
omitted since, without considering gluonic corrections, it appears as an anomalous term in 
the O(P^) chiral Lagrangian [35], and will be present only in two-loop corrections according 
to power counting. For very large values of fiB, one has to incorporate effects from the color 
superconductivity gap in the calculation of meson masses in an effective theory near the 
Fermi surface [13, 36]. For simplicity, we do not include these corrections, and restrict our 
analysis to more modest values of baryonic chemical potential. 

The nucleon masses, Mp = M — 6M and M„ = M + 6M, are dressed by leading- 
order contributions in zero-temperature baryon chiral perturbation theory. Using the results 
from Ref. [37], for the isospin symmetric case with explicit chiral symmetry breaking, and 
Ref. [38], which includes explicit isospin breaking effects, we have (neglecting terms ~ m^, 
~ log(mq), and of higher order in m^): 

M(m) = Mo + 2 7i m + 2^/^ ^^^^ ^3/2 ^ ^^2.6) 
5M{6m) = 2 7f^'""^5m, (2.7) 

Mq being the nucleon mass in the chiral limit, m = {rriu + md)/2 the average quark mass, 
and^ 



71 = -^{-m, (2i 

J n 



-.2 



71 - j^. (2.10) 

Here, all parameters and coefficients are fixed to reproduce properties of the QCD vacuum 
either measured or extracted from recent lattice simulations. Explicitly, {qq) = —(225 MeV)^ 
is the (1-flavor) chiral condensate in the chiral limit [39] and, from Table I in Ref. [37], 
Mo = (0.882 ± 0.003) GeV, ci = (-0.93 ± 0.04) GeV^\ gA = 1.267, and = 92.4 MeV. 



can be explored, nevertheless, using rn^/fij as an expansion parameter, even though the validity of the whole 

treatment is restricted to |^/| smaller than the rj or p masses [23]. 

^To relate the pion mass (in the isospin symmetric case) with the quark masses, we use the Gell-Mann - 

Oakes - Renner relation: iri^f^ = —{rriu + md){qq) = 2m{— (qq)) . 
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From Table 3 in Ref. [38], one can extract 

2a- f3 j (0.198 ± 0.093) (lattice units) , ©(m^) fit 

3 I (0.229 ±0.058) (lattice units) , 0(mg/^) fit 

where 1 (lattice units) = b = 0.125 fm. Therefore: 

2a -(3 f (0.12543 ± 0.0589) GeV"^ , 0(m,) fit 

(0.145064 ±0.036741) GeV"^ , 0{mf^) fit 



(2.11) 



(2.12) 



Evaluating the coefficients 7 exphcitly, we arrive at 71 = 4.9630±0.2135, 73/2 = -0.273424 MeV"^/^ 
and 

fo.l6734± 0.07858 , 0(m„) fit 
[0.19354 ±0.04902 , 0(mg/^) fit 

Finally, since the nucleon chemical potentials are given by 

3 1 

/ip = 2/i„ ± /irf = -yUfi ± -/i/ , (2.14) 
3 1 , , 

/^n = /i„ ± 2/id = -/iB - -/i/ , (2.15) 

the dispersion relation satisfied by the proton and neutron are, respectively: 

^proton(p) = ± (M - 5My ± ^flB + ^/i/ , (2.16) 

^ncutron(p) = ± (M ± SM)^ + ^fXs " ^/i/ , (2.17) 

where p is the 3-momentum and the antiparticle dispersion relations are obtained from the 
ones above by the substitution /ij i— )• — /ij. 

Using the dispersion relations for pions and nucleons, we can write the total pressure for 
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the gas of free (dressed) quasiparticles as Pchirai = -Ppions + -Pnucicons, with 

Ppions = / I log [1 - e^^I-(-±)^'^l] + log [1 - e-'^[-(-±)+^]] + 

+ log[l-e-^'"(™°)] i, (2.18) 



nucleons 



2T J ^ I log [1 + e-^mM-SM)+,,]^ ^ ^ ^.plEiM-SM)-,,]^ ^ 

+ log [1 + e-/3[^(M+5M)+^„]] ^ [1 ^ ^_;,[s(Af+5M)-M„]] 1 ^ (2.19) 



where we have defined a;(m) = y/^^Trn? and E{M) = a/p^ + M^. 



2.2.2 High-energy sector 

The Fuzzy Bag Model has been proposed by Pisarski [30] as a phenomenological parameteri- 
zation of the equation of state to account for the plateau in the trace anomaly normalized by 
T^, (e — 3p)/T^, observed in lattice results above the critical temperature. Besides the usual 
MIT-type bag constant, the total pressure for QCD in this model has also a non-perturbative 
contribution^ ~ T^: 

Pdcconf(T) ~ PpQCd(T) - 5fuzzy T'^ - BuiT ■ (2.20) 

The trace anomaly associated with this equation of state, assuming that PpQCD ~ 7"^, is then 

e - 3p = - 4p = 25f,,,y + ABuit • (2.21) 

In Ref. [19], a similar parameterization has been used to fit lattice results for the trace 
anomaly at high temperatures, T > 1.5 Tc ~ 300 MeV, yielding^ 

e — 3»\ 3 , . 6 c , , 

-f^] = i^o^?^ + ^ + ^, (2.22) 

/ high T ^ ^ ^ 



^Thc underlying theoretical framework is that of an effective theory of Wilson lines and their electric field 
[30, 40, 41]. 

^In the regime Tc < T < 1.5 Tc, the authors of Ref. [19] identify a strongly non-perturbative regime, 
which is not included in the fuzzy-type fit. 
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with (cf. Table VIII of Ref. [19]) 



(0.101 ± 0.006) GeV^ , K = 4: ; g"^ = 

(0.16 ± 0.06) GeV^ , iV^ = 6 ; g'^ = 2.3± 0.7 

(0.024 ± 0.001) GeV^ , iVr = 4 ; g^ = 

(0.013 ± 0.006) GeV^ , Nr = 6 ; g"^ = 2.3± 0.7 



(2.23) 
(2.24) 



Notice that the first term in Eq. (2.22) comes from a O(a^) perturbative contribution to the 
pressure and is important for the fit only at very high temperatures. We neglect this term 
in what follows. Comparing Eqs. (2.21) and (2.22), we extract the following values for the 
bag coefficients: 

' (0.051 ± 0.003) GeV^ , iVr = 4 ; g^ = 

-Bfuzzy = < (2.25) 

(0.08 ± 0.03) GeV^ , = Q ; g"^ = 2.3± 0.7 

(0.0060 ± 0.0003) GeV^ , iVr = 4 ; = o 
-Bmit = \ (2.26) 
(0.0033 ± 0.0015) GeV^ , iVr = 6 ; g'^ = 2.3± 0.7 

In our effective theory, we adopt, phenomenologically, a simple extension of the equation 
of state in Eq. (2.20) which includes the influence of finite chemical potentials in the per- 
turbative pressure, neglecting for simplicity non-perturbative contributions due to the finite 
quark chemical potentials /i/, so that: 

Pdeconf(7', /i/) ~ PpQCD{T,flf ) - ^f^^zy " ^MIT ■ (2.27) 



2.3 Results for the functional behavior of the decon- 
finement critical parameters 

With the explicit expressions for the pressure in the deconfined phase and in the low-energy 
(confined) phase, we can study the effects of varying the quark masses and the chemical 
potentials on the thermodynamics of the deconfining transition. In particular, we can study 
the behavior of the critical temperature, with the aim of comparing with available lattice 
results. The critical baryon chemical potential in the cold case is well within the region 
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where the Sign Problem sets in rendering lattice QCD essentially inapplicable, so that model 
predictions are very welcome. 

From our results for the massive free gas contribution of the pQCD pressure in the fuzzy 
bag model at finite temperature, isospin and baryon number, and the free gas pressure 
of quasi-pions and nucleons in the low-energy regime, the critical temperature and chemical 
potential for the deconfining phase transition are extracted by maximizing the total pressure. 
The validity of our approach is, of course, restricted by the scale of xPT: e.g. for m^r ~ 
nip 770 MeV, the expansion parameter in xPT is a ~ 0.45, so that the extension of the 
predictions to m^r < 1 GeV is in principle justified^. 

By matching the two branches of the equations of state, corresponding to the high and low 
temperature regimes, we of course obtain a first-order transition. Systematic Lattice QCD 
calculations, with realistic quark masses [18] (or almost realistic [19]), seem to indicate a 
crossover instead. On the other hand, from the experimental standpoint a weakly first-order 
transition is not ruled out, and in fact corresponds to the scenario adopted in very successful 
hydrodynamic calculations [25]. Although this is a crucial question for the understanding of 
the phase structure of QCD, it is essentially of no consequence to the analysis we undertake. 
The value of the critical temperature we obtain is ~ 5 — 10% different from current values 
extracted from the Lattice [19]*, which is always the case in effective field theory approaches 
and of no harm to our analysis, either. Our concern is with providing a good description of 
the behavior of the critical temperature for increasing values of quark masses and the isospin 
chemical potential. 

2.3.1 Case with = and comparison to lattice QCD 

Let us first concentrate on the results for a system with no baryon-antibaryon asymmetry. 
This is the domain for which lattice QCD simulations can be applied in a controlled and 
systematic way, providing a reference for comparison. 

rough estimate of the error (due to neglecting higher orders in chiral perturbation theory) in the 

effective masses of the low-energy sector of our approximation for m^r ~ 1 GeV yields ^ 30%. 

^For many years, there has been no agreement between different groups regarding the value of the critical 

temperature (cf. Ref. [18]). This important issue has been settled in 2010 and the results converged [42] to 

that of Rcf. [18]. 



CHAPTER 2. PHASE DIAGRAMS OF STRONG INTERACTIONS IN THE. 



28 







1 ■ 1 ' 1 


• \ • \ • y 






Fuzzy bag + Pion-Nucleon gas 


/ . 






Bag + Pion-Nucleon gas 


/ 






0(4) LoM 


/ 






RG improved LoM 


/ 

y 




■ 


Lattice = 2 

F 






• 


Lattice = 3 

F 

Lattice N„ = 2+1 

F 


y 

y 

y 

y 








y 














— ; — n 




, . ■ " * 

1 . 1 . 1 . 1 . 1 . 



)l 1 I 1 I 1 I 1 I 1 I 1 I 

0.5 1 1.5 2 2.5 3 



m /a 

Figure 2.2: Critical temperature as a function of m^r normalized by the string tension -y/o" = 
425 MeV. Our results are compared to lattice data [43, 19] and other approaches [27, 26]. 

The pion mass dependence, or equivalently the quark average mass dependence, of the 
critical temperature is displayed in Fig. 2.2 for vanishing chemical potentials with the 
temperature and the pion mass normalized by the square root of the string tension. Our 
curves stop at a point from which our "xPT approach clearly breaks down. The results of 
our framework are compared to lattice data from Refs. [43] {Np = 2, 3) and [19] {Np = 
2 + 1) and two other phenomenological treatments: the 0(4) linear sigma model [27] and a 
renormalization group improved computation [26] (cf. also Ref. [28] that discusses the quark- 
meson model using the proper-time renormalization group approach). The approximate mass 
independence observed in the lattice data is well reproduced within our framework, while 
the other descriptions tend to generate a qualitatively different behavior. 

This feature is yet another indication that the functional dependence of Tc(mjr) requires 
confinement ingredients to be reproduced, being incompatible with a phase transition dic- 
tated by pure chiral dynamics. This argument goes in the same direction of Ref. [27], 
the main difference here being the fact that we construct the mass dependence from the 
^quar-k = limit, with the heavy quasi- nucleons as the key new element at low-energies. 

Moreover, our results are not strongly sensitive to the choice between the fuzzy bag model 
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and the usual MIT bag model. The critical values for the MIT bag model are systematically 
lowered, but the qualitative behavior is not altered, as illustrated in Fig. 2.2. This indicates 
that a consistent treatment of the quark mass dependence connecting both perturbative 
regimes of energy is probably the essential ingredient to describe this observable. 
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Figure 2.3: Critical temperature as a function of fij for two different values of the pion mass. 
Our results are compared to lattice data [22] and other approaches [44]. 

In Fig. 2.3, the critical temperature is plotted as a function of the isospin chemical 
potential. The critical temperature is normalized by its value in the absence of chemical 
asymmetry, whereas the isospin chemical potential is normalized by the critical temperature 
itself. The curves obtained within our framework are represented by solid lines. Once again 
our results are in good agreement with lattice computations [22], even though the curve 
that is closer to the lattice points corresponds to a small vacuum pion mass, which is not 
the situation simulated on the lattice [22]. Our curves for different values of m^r start to 
disagree more appreciably for fij > 2Tf.. Recall that our treatment is valid up to isospin 
chemical potentials that are larger but not much larger than the pion mass, and contains the 
effects from pion condensation for /i/ > m^. Fig. 2.3 also exhibits, for comparison, results 
using the hadron resonance gas model [44] that appear to depart from the lattice data at 
a much lower value of the isospin chemical potential. The two curves are produced by two 
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different methods to determine the critical temperature: the dotted curve is obtained from 
the observation that the deconfined phase emerges at a constant energy density, whereas the 
dashed one uses the fact that the quark-antiquark condensate for the hght quarks almost 
disappears at the quark-hadron transition [44]. Similarly to what we observe for the quark- 
mass dependence of Tc, it is clear from Fig. 2.3 that plain chiral considerations render 
the largest discrepancies for the behavior of Tc(yU/) as compared to lattice data. It is also 
interesting to notice that contrary the mass- dependence of Tc, the behavior of Tc with fij 
discriminates between the phenomenological bag treatments adopted at high energies: the 
effective model with the usual bag pressure does not reproduce lattice data as well as the 
fuzzy one: even for m^r = 25 MeV, the bag model curve is still between the = 400 MeV 
and = 600 MeV fuzzy results. 
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Figure 2.4: Critical temperature as a function of fij for different values of 171^^. Lattice data 
from [22]. 



An up-down quark mass imbalance, characterized by the relative difference in quark 
mass, which is in principle not much smaller than one, tends to increase the critical tem- 
perature, though by a quantitatively small amount, as expected. This result is presented in 
Fig. 2.5 as a function of the isospin chemical potential. At finite fij, the finite difference 
between the masses of quarks up and down spoils the cancelation of imaginary parts that 
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caused the fermionic determinant in isospin dense QCD media to be positive definite. This 
situation suffers therefore from the Sign Problem even at zero baryon-chemical potential, 
being essentially inaccessible to current lattice algorithms. 




Figure 2.5: Critical temperature as a function of ni for different values of the quark mass 
difference dm. 

2.3.2 Including a finite baryon density 

Now that we have shown that the model yields reasonable predictions for the behavior of the 
critical parameters as compared to available QCD lattice data, one may explore the domain 
of the phase diagram in which a baryon- antibaryon asymmetry is present. 

The dependence of the critical temperature on the pion mass is shown once again in Fig. 
2.6, where the influence of a finite baryon chemical potential is also displayed. As discussed 
previously, the approximate mass independence observed in the lattice data [43] is well 
reproduced here, while previous treatments tended to generate a rather different behavior 
(cf. Fig. 2.2). The effect of a finite /is is completely imperceptible for /i^ < 100 MeV. We 
also computed the critical parameters by considering the usual bag model in the high energy 
regime, finding values systematically lower, but with the same qualitative behavior as with 
the fuzzy bag. 



CHAPTER 2. PHASE DIAGRAMS OF STRONG INTERACTIONS IN THE. 



32 



1.5 



1 - 



_ _ . Fuzzy bag + Pion-Nucleon gas (|Xg = 200 MeV) 

. . . Fuzzy bag + Pion-Nucleon gas (|Xg = 400 MeV) 

Fuzzy bag + Pion-Nucleon gas 
■ Lattice Np = 2 = 0) 
• Lattice Np = 3 (^^ = 0) 

Bag -I- Pion-Nucleon gas (p^ = 0) 



0.5 




0, 







0.5 



1.5 



2 



2.5 



3 



m /a 

7t 



1/2 



Figure 2.6: Critical temperature as a function of m^r normalized by the string tension ^/a = 
425 MeV, for different fiB (/i/ = 0). Lattice data from [43]. 

Fig. 2.7 displays the cold, dense case: the value of baryon chemical potential beyond 
which matter is deconfined increases with the vacuum pion mass. The outcome of a finite, 
small isospin density is to decrease the critical baryon chemical potential. This property 
could in principle play an important role in astrophysics. Interestingly, the critical baryon 
chemical potential seems to be significantly more sensitive to mass variations in comparison 
with the finite temperature deconfinement transition. The relevance of mass effects in the 
thermodynamics of cold, dense matter has already been shown in different related contexts, 
such as perturbative QCD [45, 46] and perturbative Yukawa theory [47] (cf. also Ref. [48]). 



The role played by physically reasonable values for nonzero quark masses, as well as for finite 
isospin and baryonic chemical potentials becomes more and more important nowadays in the 
understanding of the phase structure of strong interactions. Present lattice simulations are 
aheady capable of probing the region of nonzero (small) values of fis/T and more realistic 
quark masses. Since the isospin chemical potential does not suffer from the Sign problem, 
this road is open for detailed Monte Carlo studies. 



2.4 Final Remarks 
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Figure 2.7: Zero-temperature critical baryon chemical potential as a function of m^, for 
different values of the isospin chemical potential. 

On the experimental side, intermediate-energy experiments in nuclear physics are pro- 
viding data and observables that are sensitive to chemical asymmetry, whereas high-energy 
heavy ion collisions to start soon at FAIR-GSI will probe a region of the phase diagram of 
QCD where effects from fiB become important. Hence, effective theories that can encode 
the relevant characteristics of the rich phenomenology that emerges are called for. 

In this chapter we constructed an effective model profiting from results from chiral per- 
turbation theory, perturbative QCD, and the fuzzy bag model to build a gas of quasiparticles 
that exhibits most of the relevant ingredients, and yet remains simple. The model implements 
a clear and consistent connection between the high and low energy parameters, allowing for a 
systematic investigation of the influence of finite quark masses and isospin symmetry break- 
ing on the phase diagram for QCD at finite temperature and density. From a simple free 
gas calculation of the equation of state, we found surprisingly good agreement with different 
lattice data for the behavior of the critical temperature with masses and isospin chemical 
potential, indicating that the model captures the essential features brought about by the 
inclusion of those effects. Furthermore, this model can supply information about the cold 
and dense regime of strongly interacting matter, which is not yet fully covered by lattice 
simulations. 

As discussed in the beginning of the chapter, to predict, using effective models, the 
functional behavior of the critical temperature in agreement with lattice results is not at all 
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a trivial task. Several detailed studies of chiral models failed to describe Tc{mT^) [26, 27, 28], 
while Polyakov-loop models, whose results can only be fitted to the lattice points for Tc{m.^) 
[27], cannot at the present form address isospin effects (nor any other observable related to 
fermionic chemical potentials). 

Finally, the fact that different effective models seem to provide distinguishable behaviors 
for Tc(m^) together with the capability of first-principle lattice QCD simulations to assess 
"cousin" theories with various masses strongly suggests that the full function Tc(m^, {/ij}) 
is actually a good observable of the physics driving the phase boundaries. If thoroughly 
explored, it could even provide information on the puzzling connection between chiral and 
deconfinement transitions. 



Chapter 3 




PSEUDOCRITICAL PHASE DIAGRAM 

OF Strong Interactions in 

FINITE SYSTEMS AND CEP SEARCH 



In this chapter, we analyze how the QCD phase diagram is reahzed in real, finite systems. 
This amounts to considering a three-dimensional diagram, including the typical size of the 
system as an extra control parameter, and investigating how and where in the T — fiB plane 
we should expect a pseudocritical behavior to appear for each system volume. 

As the discussion that follows will show, this should be an important issue when assessing 
the QCD phase transitions via HICs. In particular, we address different aspects of finite- 
size effects in the experimental search for the QCD critical endpoint in HICs, ranging from 
the addition of obstacles [49] (further discussion in Refs. [50, 51]) to the provision of an 
alternative signature [52] (cf. also [53]). 
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3.1 The critical endpoint in the QCD phase diagram 

Before addressing the problem of assessing the QCD phase structure through real, compli- 
cated systems as HICs, let us discuss the general picture in the thermodynamic limit and 
point out why we expect a second order critical point to exist in the first place. For an 
extended discussion and further literature the reader is referred to Ref. [54] and references 
therein. 

The QCD phase transitions at high energies are usually described in terms of chiral 
and/or confinement properties. The symmetries associated with both of them (namely, 
SU{Nf) X SU{Nf) and Z{Nc)) are explicitly broken in Nature by the presence of quarks and 
their nonzero current masses, respectively. No exact order parameter can be thus defined in 
real-world QCD. Nevertheless, one clearly expects a drastic change in the medium of QCD at 
very high energies, in which a much larger number of degrees of freedom should be dynamical. 
This expectation has been intensively tested in lattice QCD simulations that indeed find a 
smooth but sensible increase in thermodynamical properties such as pressure and entropy 
for very hot media [55]. At zero baryon chemical potential, first-principle lattice results 
strongly indicate a crossover transition between hadron-chirally broken and quark- chir ally 
(approximately) restored phases. At large /ig, as discussed in the last chapter, the lattice 
QCD reference is lost^ and one has to resort to effective approaches. The great majority of 
(chiral or mixed chiral-confinement^) models predicts a first-order phase transition for cold 
and dense systems. In order to compatibilize these theoretical expectations for the phase 
diagram, a second-order critical endpoint has to be assumed at a finite /i^, below which the 
transition becomes a crossover. 

If we concentrate on chiral properties and the up-down sector, this CEP "topology" of 
the phase diagram may be viewed as reminiscent of the chiral limit, in which the small 
explicit symmetry breaking due to the current quark masses is neglected. In this case, with 
'^quark = 0, chiral symmetry is exactly preserved by the lagrangian and the chiral condensate 

^Some approximate lattice approaches to the baryon-dense domain also predict critical point [56] (cf., 
however, Refs. [57]). 

^It should be noted that currently known pure confinement models cannot describe baryon-antibaryon 
asymmetric media in a simple way. 
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is a well-defined order parameter. There is thus a sharp distinction between the two phases 
and a singularity (or finite discontinuity) must be reached when passing from one of them to 
the other; no crossover transition is allowed. Indeed, two-fiavor chiral QCD belongs to the 
same universality class as the Ising Model, featuring a second order phase transition at finite 
temperature (as also verified on the lattice^; cf. e.g. [55]). This second order line turns into 
a first-order one at higher /i^'s; the intersection being a tri critical endpoint. With this chiral 
picture of the phase diagram in mind, one can turn on the current quark masses mquark- 
The clear distinction between the phases is lost and all transitions get smeared: while the 
second-order one becomes a smooth crossover, the first-order transition remains, but with a 
reduced discontinuity. The CEP is the smoothed version of the tricritical point present in 
the chiral limit. 

Of course, these theoretical indications are not proof of existence. Any experimental 
guidance would be invaluable in this quest for mapping the true QCD phase diagram. That 
is the opportunity provided by the current and near-future heavy-ion collision experiments. 

After more than a decade of experiments colliding heavy nuclei to investigate the proper- 
ties of matter under extreme conditions of temperature and density, a clear signal of collec- 
tivity at the partonic level observed at RHIC-BNL together with other indications point to 
the formation of a new state of matter, even though its properties seem to be more intricate 
than those of the originally predicted quark-gluon plasma [10]. Current and near-future col- 
liders stretch the currently probed frontiers in two directions. The higher energies attained 
at the LHC provide the closest conditions to the primordial universe plasma and the rich 
statistics will allow for better studies of fiuctuations and the properties of the medium. In the 
other direction, different experimental programs are dedicated to probe ultra-dense media 
with higher baryon-antibaryon asymmetry: the on-going Beam Energy Scan (EES) program 
[59] at RHIC-BNL and the future colliders at FAIR-GSI and NICA-JINR. This region of the 
phase diagram is exactly the one expected to contain the richer part of the phase structure 
of strong interactions. 

■^Actually, there is still some controversy in the literature regarding this result. In Ref. [58], for instance, 
the authors argue that the transition may be dominated by confinement properties, rendering it of first order 
instead. This is yet another indication of the subtleties sometimes involved in extracting continuum results 
from lattice simulations. 
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In particular, the establishment of the existence and location of the critical endpoint of 
QCD (CEP) is a crucial step in the mapping of the phase diagram. The investigation of the 
experimental signatures of the CEP in the specific setup of heavy-ion collisions (HICs) is 
therefore of great importance. In this chapter, we discuss some of the proposed signatures for 
the CEP in HICs, concentrating on the possible role played by the finiteness of the system 
created in these collisions, as proposed in Refs. [49, 52]. 

3.2 From the thermodynamic limit to real systems 

In the idealized case of an equilibrium system in the thermodynamic limit, a critical endpoint 
is always associated with unambiguous, sharp signals due to the second order phase transition 
that characterizes it. Large fluctuations are expected at all length scales and the correlation 
length diverges, leading to the conformal invariance present in second-order phase transitions 
(cf. e.g. Ref. [60]). Direct consequences and signals are then divergences in correlation 
functions of the order parameter a (e.g. susceptibilities), that scale with the correlation ^ 



where pn is a positive exponent related to the specific correlation function and to universal 
aspects of the system under investigation. 

In any real system, however, even if we can assume equilibrium, the volume is finite, so 
that the partition function (and hence all other thermodynamic quantities) is always well- 
defined and nonsingular. Similarly, the correlation length is trivially bounded by the finite 
volume of the system^ and no criticality can be defined. Instead of the divergent critical 
behavior expected for correlation functions in the thermodynamic limit, the observables 
measured will show pseudocritical peaks, corresponding to a smoothening of the associated 
critical singularity. As illustrated in Figure 3.1, these peaks are actually shifted from the 
position of the original singularity by a size- dependent amount [60, 61]. 

''The correlation length is also bounded, through causality, by the finite lifetime of the system. Actually, 
in principle, one has to consider the growth rate of correlated domains in a fully dynamical approach, the 
speed of light (i.e. causality) being the ultimate bound in a relativistic system. 




(3.1) 
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Figure 3.1: The n-point correlation function of the order parameter a is shown as a function 
of the external parameter T. The singular (black) curve corresponds to the thermodynamic 
limit, while the behavior in a finite system is illustrated by the shifted peak (red curve). 

In this case, the only way to precisely establish the presence of a second-order phase 
transition is to verify that the correlations in the system satisfy a specific behavior as a 
function of the system size leading to the conformal invariance that characterizes the second- 
order phase transition in the thermodynamic limit. Such specific trend is the well-known 
Finite-Size Scaling (FSS) [61]. 

3.3 CEP signatures in HICs 

The connection between CEP phenomena in the thermodynamic limit and the asymptotic 
particles measured in HIC detectors is even more indirect. The medium created in such 
experiments is an extremely complicated and fastly evolving system that, as illustrated in 
Figure 3.2, goes through different stages that cannot in general be disentangled in the final 
observed spectrum. The QCD phase transition is only one of these stages, being followed 
by an interacting hadronic medium before the freeze-out that ultimately defines the final 
particle distributions detected. Moreover, it is probable that the phase conversion process 
occurs actually far from equilibrium, which could change completely the correlations in the 
system and even manifest features which are not consistent with the ones expected from the 
phase structure in the thermodynamic limit [63]. 
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Figure 3.2: [extracted from Ref. [62]] Schematic cartoon of the stages occuring during a 
heavy-ion colhsion. 

The problem of defining a clear and unambiguous signature for the CEP in HICs is 
therefore a nontrivial task that should be pursued with caution, always having the various 
caveats resulting from the specific experimental setup in mind. In what follows, we will 
disregard completely several of these caveats (e.g. out-of-equilibrium effects, critical slowing 
down near the CEP [64], hadronic medium modification, etc) and concentrate on the role 
played by finite-size effects. Besides discussing the caveats brought about by the finiteness 
of the system, we show how the fact that HICs provide data from media of different sizes 
may actually allow for a complementary alternative signature of the CEP based on FSS [49]. 

As discussed above, the CEP in a real system manifests as peaks in the cumulants of the 
distribution of the order parameter as functions of the distance to the criticality in the exter- 
nal parameter space. In this vein, it is natural to try to construct signatures of the CEP in 
HICs related to non-monotonic behavior in observables connected to the order parameter of 
the QCD transitions as one passes through or nearby the critical temperature Tc and baryon 
chemical potential fiBc- In HICs, however, one does not have direct experimental access to 
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neither of the (approximate) order parameters for the chiral and deconfinement transitions, 
namely the chiral condensate a and the expectation value of the Polyakov Loop, respectively. 
Moreover, even if we assume that the system is equilibrated, direct measurements of T and 
fiB are not available. 

The first concrete proposal of signature that tries to address these issues is the one by 
Stephanov et al [65]. The idea is based on the assumption that the critical correlations of 
the chiral order parameter will be transmitted to particles in the final HIC spectra that can 
effectively couple to the chiral condensate field a, such as pions (via the interaction Gavr^Tr") 
and nucleons {gN<yNN). Within this effective description, the critical contribution to the 
correlation of fluctuations in the particle numbers can be written in terms of the correlation 
length ^ (or the inverse of the a mass), that should drastically increase (or diverge in the 
thermodynamic limit) at the CEP. The leading term in two-particle correlations comes from 
the exchange of a a, as depicted diagrammatically in Figure 3.3, giving (for more details 
and notation cf. e.g. Ref. [66]): {6np6nk) ~ ^/f^a ~ Thus, a peak in this observable 
is expected in the vicinity of the CEP. Its height is related to the maximum correlation 
length, which in turn is limited by various factors: the size of the nuclei overlapping region 
(in practice, the effective system size can be considerably small, as for example is the case in 
core-corona scenarios [67]), the finite lifetime of the medium, critical slowing down [64], etc. 




Figure 3.3: Diagram of the leading contribution to two-point pion correlations featuring the 
transfer of critical correlations of the order parameter field a to pions. 

Besides the critical contribution to the particle number fluctuations, there are also differ- 
ent noncritical contributions (e.g. those coming from thermal decays) that form a background 
to the CEP signature which could consequently be hidden. With the aim of increasing the 
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sensitivity with respect to the critical behavior near the CEP, different observables that cou- 
ple to higher cumulants of the order parameter have been proposed. Respectively connected 
to the three- and four-particle correlators, the critical contributions to skewness co^ ~ ^^/^ 
and kurtosis U4 ~ have stronger dependence on the correlation length, yielding a higher 
chance that these peaks surmount the noncritical background [68]. 
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Figure 3.4: Measurements of second, third and fourth cumulants of fluctuations. Preliminary 
data from the Beam Energy Scan program at RHIC. 



Specific dimensionless ratios of these higher cumulants that are believed to be less model- 
dependent [66] have been compared to current available HIC data and, very recently, to the 
analysis of the first run of the RHIC BES program. Some results are shown in Figure 
3.4. Up to now no clear indication of criticality has been observed in this type of signatures. 
Instead of the expected increase in the fluctuations at lower energies, one has actually found a 
decreasing trend in the kurtosis observable. Recent work has been put forward that interprets 
this apparent drop as the behavior in the direction of negative kurtosis, which is claimed to 
be a signature of the vicinity of the CEP from the crossover side [69] . Nevertheless, this is 
not a sufficient condition for the CEP and can be a result of other constraints of the medium 
created in HIC, such as baryon number conservation 



3.4 Role of finite-size effects in HICs and CEP search 

First of all, the finiteness of the system will only play a significant role if it is comparable 
to the relevant scales associated with the QCD phase transitions under investigation. The 
typical critical temperatures for deconfinement and chiral symmetry restoration are predicted 
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by lattice QCD to be Tc ~ 150 MeV [55], which gives a "critical size scale" Lc ~ 1.3 fm. If we 
consider the early universe during the epoch of the primordial quark-hadron phase transition, 
one can estimate roughly its radius as given by the particle horizon in a Robertson- Walker 
space-time (in which the scale factor grows as a{t) ~ t"). Using a simple approximation for 
the equation of state (3P ^ e ^ (7^V30)A^T^ with being the total number of thermal 
degrees of freedom), one gets: 

, m~ 1 /" 1 V 45 y-"" Mpi.„d, ^ 1.45x10" ^ 

or Luniv ~ 10^® fm for QCD {n = 1/2, ~ 50 and T^ ~ 150 MeV). The primordial universe 
corresponds thus to a system of typical size several orders of magnitude bigger than the 
critical size scale being essentially in the thermodynamic limit as far as the QCD phase 
transitions are concerned. 

Although the process of phase conversion of hot hadronic matter that presumably hap- 
pens in ultra-relativistic heavy-ion collision experiments is often compared to the cosmo- 
logical quark-hadron transition in the early universe, the relevant space-time scales differ 
by almost twenty orders of magnitude. In HIC experiments, the volume of the system is 
not only finite but also centrality-dependent [10]. The typical linear sizes L are actually 
quite small {L < d^on, with the ion diameter dion ~ 10 — 15 fm), corresponding to energies 
> 13.1 — 19.7 MeV, not negligible in comparison with the expected scales in the critical 
region of the phase diagram < 150 MeV, especially in non-central collisions. This suggests 
that finite-volume effects may affect significantly the different physical phenomena occuring 
in a heavy-ion collision. 

In what follows, we address the role played by finite-size effects in the QCD (pseudocrit- 
ical) phase diagram as probed by HICs and in the context of the signatures of the CEP. In 
the next subsection, we show, using a well-defined and established effective model for the 
chiral transition, that indeed the expectation put forward by the rough estimates above is 
verified: for system sizes realized in current HICs, the (pseudo) critical region of the chiral 
phase diagram differs sensibly from the result in the thermodynamic limit [49]. Subsection 
3.4.2 presents then the consequences for the proposed observable signatures of the CEP in 
HICs; finite-volume effects representing in general an extra obstacle. In Subsection 3.4.3, we 
argue that the presence of large finite-size corrections in HICs data allows for an alternative 
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signature for the CEP search [49]. The practical implementation of the method of data 
analysis is demonstrated [52] in Subsection 3.4.4. 

3.4.1 Volume dependence of the chiral pseudocritical phase dia- 
gram 

Let us consider the pseudocritical phase diagram for the chiral transition in Strong Interac- 
tions, including an extra axis for the system size. 

Finite-size effects in the study of phase transitions via lattice simulations, seen as an 
inevitable drawback of the method, and the necessary extrapolation to the thermodynamic 
limit have been thoroughly studied for decades [60], the ultimate answer to the problem 
being given by the method of FSS [61]. Systematic calculations of finite- volume corrections 
can then be computed not only near the criticality of continuous (second-order) transitions, 
but also for first-order phase transitions [71], so that the thermodynamic limit can be taken 
in the calculation of properties of the phase diagram. 

The study to be presented here goes in the opposite direction: for natural systems that are 
truly small, one should study the modifications caused by its finiteness in the phase diagram 
before comparing to experimental observables. As discussed above, this is the situation in 
the study of the chiral transition in HICs. Nevertheless, contrasting to the large number of 
studies of finite-volume corrections in the computation of the chiral condensate and related 
quantities within chiral perturbation theory (see e.g. [72] and references therein), this issue 
is often overlooked in the case of the quark-gluon plasma. Among the exceptions, there is 
a lattice estimate of finite-size effects in the process of formation of the quark-gluon plasma 
[73], a few studies within the Nambu-Jona-Lasinio model [74], and an analysis of color 
superconductivity in dense QCD in a finite volume [75], besides investigations of finite-size 
effects on the dynamics^ of the plasma [77, 76] . The effect of finite volume on the pion mass 
and chiral symmetry restoration at finite temperature and zero density was investigated by 
Braun et al. using the proper-time renormalization group [78]. 

In this subsection, we present results [49] (cf. also Refs. [50]) for the modification of 

^For simplicity, we completely ignore the dynamics of the phase conversion. Nevertheless, it was shown 
in Ref. [76] that finite-size effects will play a relevant role in the case of HICs. 
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the pseudocritical line in the temperature-chemical potential phase diagram, as well as for 
isentropic trajectories, employing an effective chiral model: the linear sigma model coupled 
to quarks with two flavors, A'^^ = 2 [79], which is presented in detail in Appendix A. This 
effective theory has been widely used to describe different aspects of the chiral transition, 
such as thermodynamic properties [80, 81, 82, 83] and the nonequilibrium phase conversion 
process [84], as well as combined to other models in order to include effects from confinement 
[85]. We focus on the chiral phase transition, avoiding the inclusion of confinement ingredi- 
ents (although the linear sigma model coupled to quarks contains mesons and quarks that 
alternate dominance in the two different phases), and show that the volume-dependence of 
the phase diagram of the model in the regime of energy scales probed by current heavy ion 
experiments can be large. Even though the output numbers in an effective model should 
not be taken as accurate predictions, the quantitative question that matters here, and that 
can be answered within this approach, is how big the dislocations of the pseudocritical line 
and critical point, isentropic trajectories, etc are relative to the same phase diagram in the 
thermodynamic limit. As mentioned above they turn out to be non- negligible. 

Strictly speaking, only systems that are infinite in volume exhibit "true" spontaneous 
symmetry breaking [86]. Nevertheless, for finite systems of volume V the barrier penetration 
factor is already suppressed as exp(— aV^), a being a positive constant determined by the 
microscopic features of each system. Of course, if one waits long enough, the system will be 
"distributed" between the two minima. If this suppression is very effective, the two minima 
are physically well separated and one can consider them as two almost different phases of 
the system, even if there is no true phase transition. 

Keeping this caveat in mind, our phenomenological analysis is motivated by the fact 
that heavy ion experiments strongly suggest a new phase of deconfined matter [10] and the 
system is obviously finite, and by the enormous success of this methods in the description 
of {pseudo)ph.ase transitions in small systems [60, 61]. For simplicity, we disregard inhomo- 
geneity corrections. Even though we are aware that these are relevant (see Refs. [83, 87]), 
they should not modify the general picture we present or diminish the significance of the 
(large) finite size effects we obtain in this simpler description. 

The effective theory adopted is described in Appendix A. In the case of a finite system of 
linear size L, the momentum integral from the one-loop quark contribution to the effective 
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potential (cf. Eq.(A.8)) is substituted by a sum 



^ = TTFW E H (1 + e-^^^-'^y^) + log (1 + e-(^^^>^y^)] , (3.3) 




k 



where Ei 




'nT'eff, and Tn(.ff — g\(T\ is the effective mass of the quarks^. It is clear that 



finite-size effects in a quantum field theory are, in a way, very similar to thermal effects with 
1/L playing the role of the temperature [88] . However, in the latter boundary conditions are 
determined by the spin-statistics theorem, whereas in the former there is no clear guidance 
(see, e.g., discussions in Refs. [74, 89]). Choosing periodic boundary conditions (PBC) are 
useful to focus on bulk properties of the model undisturbed by surface effects [88, 90, 91], 
but one is always free to choose anti-periodic boundary conditions (APC) or any other, 
instead. Depending on the physical situation of interest, a boundary condition may be more 
realistic than others. To illustrate the dependence of the pseudocritical phase diagram on the 
boundary conditions, we show results for PBC and APC''. For PBC, e.g., the components 
of the momentum assume the discretized values ki = 27r£j/L, ii being integers, and there is 
a zero mode (absent for APC, for which ki = n{2ii + 1)/L). At zero temperature, this zero 
mode in the case with PBC will modify the classical potential, generating size-dependent 
effective couplings and masses. 

In what follows, we analyze systems of sizes'^ between 10 fm and 2 fm which relate to 
the typical linear dimensions involved in central and most peripheral collisions of Au or 
Pb ions at RHIC and LHC, respectively. Figure 3.5 displays the shift of the pseudocritical 
transition lines and their respective endpoints as the size of the system is decreased. The 
transition lines represent pseudo-first-order transitions, characterized by a discontinuity in 

^In general, the effects of finite size are not fully accounted for by the replacement of continuum integrals 
by discrete sums over momentum states. In certain cases, the importance of inhomogeneities can be enhanced 

in a small system, especially near its surface. For simplicity, we neglect these effects here. 

^Lattice studies [92] of the pure-gauge SU(3) deconfining transition in a box with wall-type boundary con- 
ditions show sensibly larger effects as compared to the analogous case with PBC. We also do not investigate 

the influence of shape variations, choosing a cubic system. 

®Our range of values for the linear size L is motivated by the estimated plasma size presumably formed 

in high-energy heavy ion collisions at RHIC [10]. The upper limit is essentially geometrical, provided by the 

radius of the nuclei involved, whereas the lower limit is an estimate for the smallest plasma observed. 
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Figure 3.5: Pseudocritical transition lines and endpoints for different system sizes within the 
hnear a model with periodic (left) and antiperiodic (right) boundary conditions. 

the approximate order parameter, the chiral condensate a, and the production of latent heat 
through the process of phase conversion. We find that those lines are displaced to the region 
of higher fi and shrinked by finite-size corrections. The former effect is sensibly larger when 
PBC are considered, indicating that the presence of the spatial zero mode tends to shift 
the transition region to the regime of larger chemical potentials. Both boundary conditions 
reproduce the infinite-volume limit for L > 10 fm. Figure 3.6 shows the corresponding 
displacement of the pseudocritical endpoint, comparing PBC and APC: both coordinates 
of the critical point are significantly modified, and ficEP is about 30% larger for PBC. For 
/i = 0, the crossover transition is also affected by finite-size corrections, increasing as the 
system decreases, as shown in Figure 3.7. Again, PBC generate larger effects: up to ~ 80% 
increase in the crossover transition temperature at yU = when L = 2 fm. 

Results for the isentropic trajectories are shown in Figure 3.8, comparing the infinite- 
volume limit with the finite system with L = 2 fm in the cases with PBC and APC. For 
sufficiently high temperatures, the isentropic lines in the thermodynamic limit are repro- 
duced, while large discrepancies are found around and below the transition region. The 
remarkable variations at low temperatures are due to the shift of the pseudocritical line and, 
for PBC, to the zero-mode-induced modifications of the vacuum properties, especially the 
vacuum constituent quark mass. Although there is no strong focusing effect around the 
critical endpoint (as observed previously for this chiral model [81] and similar ones [93]), 
it is clear that the density of isentropic trajectories around the pseudocritical endpoint is 
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Figure 3.6: Displacement of the pseu- 
docritical endpoint in the T — fi plane 
as the system size is decreased for dif- 
ferent boundary conditions. 



Figure 3.7: Normalized crossover 
temperature at /i = as a function 
of the inverse size 1/L for the cases 
with PBC and APC. 



sensibly higher than in the thermodynamic limit. 



3.4.2 Consequences for signatures of the CEP 

As already stated, the finiteness of the system created in HICs imposes a natural limitation 
on the growth of the correlation length. Besides this trivial constraint that smoothens out 
CEP divergences into peaks, the finiteness of the system created in HICs also guarantees 
that non-monotonic signals will actually probe pseudocritical quantities which might be 
significantly shifted from the genuine criticality, as illustrated in Fig. 3.1. Our results for 
these shifts in a chiral model clearly indicate that the displacement of the (pseudo-)CEP and 
the 

(pseudo-) first-order lines can be large for the size scales present in current HIC experiments 
[49] (cf. Fig. 3.5). If this is indeed the case in QCD, the non-monotonic behavior signaled 
by the fluctuation observables will not only be blurred in a region and the effects from 
criticality severely smoothened, but also this region will be shifted with respect to most of 
the effective-model and lattice predictions, that are done in general in the thermodynamic 
limit. Moreover, once finite-size corrections are important, the observables will be also 
sensitive to boundary effects. Smoothening of the fluctuation peaks will be further increased 
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Figure 3.8: Isentropic trajectories, labeled by the respective value of entropy per baryon 
number, in the thermodynamic limit (solid, red lines) and for L = 2 fm with PBC (left) and 
APC (right). The red dot is the genuine critical endpoint, while the square is pseudocritical 
one. 

by the averaging procedure within not sufficiently small centrality (i.e. system size) bins, 
as illustrated in Figure 3.9, tending to hide them underneath the background. Of course, 
the width of the centrality bins in data analysis is also bounded from below so that the 
statistics is enough to guarantee reasonable statistical errors. Figure 3.8 indicates further 
that signatures based on the focussing of isentropic trajectories could probe a quantitatively 
different scenario due to the finiteness of the systems created at HICs. 

superposition of shifted peaks 
broadening of si g nal 

(a") 




T 



Figure 3.9: Cartoon illustrating the broadening of a pseudocritical signal due to the averaging 
over a centrality interval. 
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3.4.3 Finite-size scaling as a tool for the CEP search 

Now that we have argued that finite-size corrections to critical phenomena in the context of 
HICs can be sizable, we can approach this issue from a different perspective with the aim of 
taking advantage of this feature in the experimental search for the CEP. 

As is well-known [60, 94], the second-order transition that occurs at the CEP in the 
thermodynamic limit is characterized by a divergent correlation length and by the property 
of scale invariance. These constraints in the thermodynamic limit are translated into real, 
finite systems as the existence of finite-size scaling (FSS) [61, 90, 91] in the vicinity of 
criticality. The phenomenon of FSS can be rigorously proved through a renormalization 
group analysis and is extensively studied and successful in condensed matter physics (cf. 
Ref. [95] for an example within the context of spin glass transitions in disordered Ising 
systems). Since the conformal invariance expected at the CEP in the thermodynamic limit 
leads to strong constraints for the behavior of finite systems, large finite-size effects as those 
encountered in Subsection 3.4.1 may, instead of only bringing obstacles for the CEP search, 
be turned into a tool. 

FSS is a powerful statistical mechanics technique that prescinds from the knowledge of 
the details of a given system; instead, it provides information about its criticality based 
solely on very general characteristics. Furthermore, it is valid also when one is not so close 
to the critical point^, provided that we use full scaling plots. This technique is predictive 
even for tiny systems in statistical mechanics (see, e.g., Ref. [95]), and, since it can be used 
in various different regions of the phase diagram not far from the critical point, it can in 
principle provide enough statistics for data analysis. 

And since the thermal environment corresponding to the region of quark-gluon plasma 
formed in heavy ion collisions can be classified according to the centrality of the collision, 
events can be separated according to the size of the plasma that is created. Heavy ion 
collisions indeed provide an ensemble of differently-sized systems, calling for a finite-size 
scaling analysis. In the following, we suggest an application of this method to search for the 
critical point of QCD (and possibly determining its universality class), assuming that we 

^In a system with size L and temperature T, the region where the scahng becomes important is charac- 
terized by X ~ '?oo(r) = (length sca\e)t~'^. 
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can build out of heavy ion data a set of systems of different temperature and sizes via the 
incident energies and the distribution of coUisional centrahty. 

The FSS hypothesis [61] was conjectured before the development of the renormalization 
group (RG). However, it can be derived quite naturally by applying RG techniques to critical 
phenomena [91, 94]. For a system with typical linear dimension L, the singular part of the 
free energy density scales as 



when the lengths are reduced by a factor £ in a given RG transformation, {g} are couplings 
and d is the space dimension. The fact that the system is finite is irrelevant for this imple- 
mentation, since the RG transformations are local. Near a RG fixed point, one can write 
(3.4) in terms of the right eigenvectors of the linearized RG transformation: 



so that one notices that L~ behaves like a relevant eigenvector with eigenvalue Al = i, 
and yi = 1- Here, t = [T — T^/T^. is the reduced temperature, being the temperature 
associated with the critical point. The reduced temperature t is a dimensionless measure 
of the distance to the critical point when no other external parameter is considered. In the 
presence of other external parameters, such as the baryonic chemical potential yU, one should 
redefine this distance accordingly to include the dependence on r = {fi — fic) / fic, which plays 
a role analogous to the magnetic field variable in a Ising system. and fic are defined in 
the thermodynamic limit, and the true second-order phase transition occurs for vanishing t, 
r and L^^ (and other couplings reaching their value at the fixed point, {g = g*}). 

For finite L, crossover effects become important. If the correlation length diverges as .^oo ~ 
t~'^ at criticality, where u is the corresponding critical exponent, in the case of L~^t^^ ^ 1 
the system is no longer governed by the critical fixed point and L limits the growth of the 
correlation length, rounding all singularities [60]. If L is finite, ^ is analytic in the limit 
t — )■ 0, etc and one can draw scaling plots of L/,^ vs. a given coupling g for different values 
of L to find that all curves cross aX g = g* in this limit, which is a way to determine g* . The 
critical temperature can also be determined in this fashion, since the curves will also cross 
at t ^ tc- 



mg},L-')=r'M{g'}JL-'), 



(3.4) 



Ut,r,...,L-') = r'f,{uy\riy^...,iL-') , 



(3.5) 
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This scaling plot technique can be extended, taken to its full power for other quantities, 
such as correlation functions. An observable X in a finite thermal system can be written, in 
the neighborhood of criticality, in the following form [91]: 

X{t,{gy,i;L) = L^-/''f{tL'/n, (3-6) 

where '^x is the bulk (dimension) exponent of X and {g} dimensionless coupling constants. 
The function f{y) is universal up to scale fixing, and the critical exponents are sensitive 
essentially to dimensionality and internal symmetry, which will give rise to the different uni- 
versality classes [60, 94]. To simplify the discussion, we ignore for the moment the chemical 
potential variable r. In principle, / is a function of two scaling variables, though. Again, if 
one plots the scaled observable versus t, the curves should cross at t = tc- Moreover, using 
the appropriate scaling variable, instead of t, all curves should collapse into one single curve 
if one is not far from the critical point. So, this technique can be applied to the analysis 
of observables that are directly related to the correlation function of the order parameter of 
the transition, such as fluctuations of the multiplicity of soft pions [65]. 

This is the idea of using full scaling plots to search for the critical endpoint of QCD in 
heavy ion data. Here, as in statistical mechanics searches for critical points, one can treat 
Tc, fic and the critical exponents as parameters in a scaling plot fit, the system size being 
provided by different centrality bins or by the number of participants in the collision. In 
fact, the scaling variable is defined up to L- and t-independent multiplicative factors, so that 
the knowledge of the actual size of the system is not needed. 

As mentioned previously, the correct scaling variable should measure the distance from 
the critical point, thereby involving both temperature and chemical potential, or their re- 
duced versions t and r. This produces a two-dimensional scaling function and makes the 
analysis of heavy ion data highly nontrivial. Phenomenologically, we adopt a simplification 
motivated by results from thermal models for the freeze-out region, connecting temperature 
and chemical potential. These models, which fit ratios of particle multiplicities surprisingly 
well, suggest that the temperature and chemical potential characterizing the data measured 
in a heavy-ion collision are constrained to a "freeze-out curve", as depicted in Figure 3.10. 
Since this curve is parameterized by the center-of-mass energy of the collision (a/stvat), it is 
reasonable to build our one-dimensional scaling variable from this quantity and the size of 
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Figure 3.10: [extracted from Ref. [96]] Values of temperature and net baryon density at 
chemical freeze-out obtained with statistical model fits to hadron abundances. 

the system. 

Actually, the reduction in the number of dimensions in the external parameter space 
seems very natural also from the standpoint of hydrodynamical descriptions of the evolution 
of the system. As suggested by the success of ideal hydrodynamic models, the expansion 
of the system may well be considered adiabatic. Therefore, we expect that the point of the 
emission of soft pions should be on this adiabatic line, characterized by the initial energy 
density and naturally parameterized in the (T — fi) plane by ^/sNN. 

The range of sizes that can be accessed in heavy ion collisions over which FSS can be 
tested is limited. It is reasonable to assume that for a locally thermalized quark gluon plasma 
to form, the system needs to be several times the hadronic size (several fm). The largest 
possible system will have a diameter of approximately 15 fm. Estimates from an analysis 
of HBT data from STAR indicate that from peripheral to central collisions, the system 
size changes by a factor of 3 — 4 [97]. We plot the scaled observable vs [(yU — ji^)/ l^c\ L^^" ■ 
The value of L for the various centralities can be estimated from a Glauber Monte-Carlo 
model. For our study, the overlap area is calculated from 5*^ = -Ra«(® " sinO) where 
6 = 2 cos~^(6/2i?Aii)5 where Rau is the nuclear radius for Au and h is the collision impact 
parameter. Then the length is taken to be L = 2a/S^x/V. 

A system of size Li can be compared to a system of size L2 only when both settings 
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This constrains which /i values can be compared when testing for scahng by directly compar- 
ing data rather than by extrapolations. Since a central Au + Au collision is about four times 
larger than a peripheral one (L2 = 4Li) and taking u = 2/3, we find that a measurement 
in peripheral Au+Au collisions at /ii can be compared to a central collision at fi2 when 

Hi- flc^ 8(/i2 - /ic). 

3.4.4 FSS analysis of preliminary results from RHIC Beam Energy 
Scan 

Data on various fluctuations have been measured at RHIC and the SPS over an energy 
range from ^/snn = 5 GeV to 200 GeV and the centrality dependence of px fluctuations 
has been measured at RHIC for energies from 19.6 GeV to 200 GeV. These data provide an 
opportunity to look for evidence of finite-size scaling. 

To search for scaling, we consider the correlation measure ap^,/ (pr) [98] scaled by L""'''^'^, 
according to Eq. 3.6. We consider the Pt fluctuations a^j, scaled by (pt) to obtain a 
dimensionless variable We use the correlation data measured in bins corresponding to the 
0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, and 60-70% most central colhsions. 
We estimate the corresponding lengths L to be 12.4, 11.1, 9.6, 8.0, 6.8, 5.6, 4.5, and 3.4 fm. 
The exponent z/ = 2/3 is determined by the Ising universality class of QCD and we consider 
values of jx around 1 (ignoring small anomalous dimension corrections). We also varied the 
value of 72; from 0.5 to 2.0 and found that changing 7^ within this range does not improve 
the scaling behavior. Using a smaller value of 7^ simply causes apj,/{pT) scaled by L~"''^/^ 
to drop more slowly with increasing centrality and increasing 7^ causes this drop to happen 
more quickly. 

In Figs. 3.11 and 3.12 we plot cTpj,/ {pr) scaled by L'^^^l^ vs /i for different system sizes, 
using data extracted from collisions at y^s/vTV = 19.6, 62.4, 130, 200 GeV. We use the param- 

^'^Thc STAR collaboration lias also published fluctuations using an alternative variable AtTp^ [99]. The 
results found using this variable are qualitatively equivalent to the ones for a^^ j {pr) ■ 



correspond to the same value of the scaling variable, i.e. 




(3.7) 
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Figure 3.11: Scaled o'pj./ipT) vs /i for different system sizes, and with u = 2/3 and •y^ = 1- 
Data extracted from RHIC collisions at energies y^sjyN = 19.6, 62.4, 130, and 200 GeV (linear 
fit, see text). 

eterization from Ref. [100] to convert from ^/sNN to fi. If there is a critical point at /i = /icrit; 
the curves for different sizes of the system should cross at this value of fi. However, since 
the currently available data is restricted to not so large values of the chemical potential, one 
has to perform extrapolations using fits. The scaling function / in Eq. 3.6 is expected to 
be smoothly varying around the critical point, so we fit the data corresponding to a given 
linear size L to a polynomial, but constraining the polynomials to enforce the condition 
that all the curves cross at some /i = ficnt, where /icrit is an adjustable parameter in the fit. 
This clearly assumes the existence of a critical point. In Fig. 3.11 we use a linear fit. The 
approximate energy independence of cTp^,/ (px) along with the linear fit, leads to a very large 
/i value where the curves can cross (/i ~ 3 GeV). There is no reason however to assume a 
linear fit function, so in Fig. 3.12 we also try a second order polynomial. Using a second 
order polynomial function for / allows the curves from different system sizes to cross at a 
much smaller value of fi. Based on this fit, we find that the data is consistent with a critical 
point at /i ~ 510 MeV corresponding to a yj/vw of 5.75 GeV. 
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Figure 3.12: Scaled dp^/ (px) vs /i for different system sizes, assuming ficrit = 509 MeV, which 
corresponds to a critical point at ^/sNN = 5.75 GeV. Again, z/ = 2/3 and 7^ = 1. Data 
extracted from RHIC collisions at energies ^/snn = 19.6, 62.4, 130, and 200 GeV (second 
order polynomial fit, see text). 

The value estimated for ficrit based on finite-size scaling of current data is highly depen- 
dent on the assumed functional form of /. The approximate energy independence of ap^/ (pt) 
for a given L, however, already indicates within the finite-size scaling assumption that the 
critical point should be at fi values well above those currently available. Based on finite-size 
scaling of o"pj,/ (pr), one would not expect a critical point at /i < 400 MeV. Having data for 
lower values of ^/s, i.e. higher values of chemical potential, as expected from the analysis of 
the Beam Energy Scan program at RHIC [59], one should be able to study full scaling plots 
of aprj,/ (pt) scaled by L""'''^'^ vs. without the need of long extrapolations. For the 

current set of data, these full scaling plots would not be very enlightening. 

RHIC has also run at lower energies in order to search for a critical point in the Beam 
Energy Scan program. That data is currently being analyzed. Here we use the quadratic 
polynomial fit of STAR data (Fig. 3.12) and assume the critical point is at 509 MeV to 
make predictions for ap^/{pT) at lower energies. The finite-size scaling scenario along with 
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currently available data, allows us to predict the energy and system size dependence of 
fluctuations for any given values of /icrit? Ix and u. We show this expectation as a function of 
the number of participants, A'part, for three proposed beam energies: 11.5, 7.7 and 5 GeV in 
Fig. 3.13. Notice that the centrality dependence changes once one moves to the other side of 
the critical point. This is a condition enforced by finite-size scaling which provides a generic 
signal for having reached the first-order phase transition side of the critical point. 
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Figure 3.13: The expected measurement of ap^/{pT) as a function of the number of par- 
ticipants at lower energies assuming the critical point is at 509 MeV as extracted from the 
quadratic polynomial fit of STAR data. 

3.5 Final remarks 

Searching for the the critical endpoint in the phase diagram for strong interactions is a 
remarkably challenging task. On one hand, lattice simulations can not provide the sort of 
guidance that is possible at zero density, where there is no sign problem, and effective models 
point to its existence, qualitatively, but yield very different quantitative predictions. On 
the other hand, heavy ion collision experiments probe limited regions of the phase diagram. 
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Besides, they come with non-trivial background contributions that tend to blur the signatures 
provided by the non-monotonic behavior of observables built from correlation functions of 
the order parameter. 

In this chapter, we point out that this task should be pursued within an extended phase 
diagram of Strong Interactions: the one containing an extra axis associated with the volume 
of the system under consideration. Finite-size effects most probably play an important role 
in HIC experiments and particularly in the QCD CEP search. As a caveat: restricting the 
access of non-monotonic signatures to pseudo-critical quantities that are severely smoothened 
and might be significantly shifted with respect to the genuine CEP (in the thermodynamic 
limit). But also as a complementary tool: FSS analysis is perfectly feasible and simple to be 
implemented in the context of HICs, yielding predictions for comparison of data for different 
center-of-mass energies. 

The fact that finite-size scaling prescinds from the knowledge of the details of the system 
under consideration, providing information about its criticality based solely on its most 
general features, makes it a very powerful and pragmatic tool for data analysis in the heavy 
ion collision search for the critical point. From a very limited data set in energy spam, we 
have used FSS to exclude the presence of a critical point at small values of the baryonic 
chemical potential, below 450 MeV. We have also used the scaling function to predict the 
behavior of data with system size at lower energies. We are looking forward to compare our 
predictions to the outcome of data analysis from the Beam Energy Scan program at RHIC. 



T 



Chapter 4 




NONEQUILIBRIUM PHENOMENA AND 
PHASE STRUCTURE IN THE COLD 
AND DENSE DOMAIN 



In this chapter we consider the phase diagram of Strong Interactions from a dynamical 
point of view. Although a very challenging enterprise, this is crucial when investigating 
phase transitions at high energies, due to the current lack of well-controlled experiments in 
this domain. In Nature or in the laboratory, the structure of phases of Strong Interactions 
is generally probed via indirect, out-of-equilibrium conversion processes. 

The particular problem being addressed here is the dynamics of nucleation for the chiral 
phase transition in the cold and dense regime. Estimates of relevant time and energy scales 
are computed within the linear sigma model with quarks [101] and the presentation focus is on 
the possible consequences of these results for different astrophysical applications [101, 102]. 
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4.1 QCD matter and astrophysical phenomena 

The thermodynamics of strong interactions for cold matter under extremely high densities 
is of utmost importance for the understanding of the structure of compact stars [32]. Since 
quantum chromodynamics (QCD) is asymptotically free, it is believed that for high enough 
densities quarks will be in a deconfined state, the quark-gluon plasma [11]. Moreover, due to 
the approximately chiral nature of the QCD action, one also expects quarks to be essentially 
massless above a sufficiently high value of the chemical potential^. Depending on the location 
of the critical density, one might find several sorts of condensates and even deconfined quark 
matter in the core of neutron stars [105]. Furthermore, the order and strength of the chiral 
transition are crucial features in establishing the existence of a new class of compact stars 
[106, 107]. 

On the other hand, going in the opposite direction, in which observations provide insight 
to theory, astrophysical phenomena might shed some light in a mostly obscure region of the 
QCD phase diagram: the domain of low temperatures and high baryon-antibaryon asymme- 
try. This region is exactly where one expects the richest structure in the high energy QCD 
phase diagram to be and where theorists have less guidance for what QCD predicts, since 
first-principle lattice calculations are still not reliable (due to the Sign Problem, cf. e.g. Ref. 
[55]). Moreover, despite the great experimental efforts to create matter with a finite baryon 
chemical potential, accelerators cannot reach the domain of very low temperature and high 
chemical potential. 

Since extremely high densities and relatively low temperatures are believed to exist in 
the core of ultracompact objects, it is probable that the QCD phase transitions play an 
important role in the structure and dynamics of formation of such systems. Even though in 
general observations are much less controlable than laboratory experiments, the astronom- 
ical measurement techniques and instrumentation have much developed in the last decade 
and this might be the dawn of a precision era in which astrophysics really constrains non- 
perturbative QCD phenomena. One concrete example of astronomical observation that has 
reached the capability of constraining models for strong interactions in the nonperturbative 

^Nevertheless, these features and other thermodynamic properties of dense media were also shown to be 
significantly affected by nonzero quark masses [103, 104]. 
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Figure 4.1: [extracted from Ref. [108]] Neutron star mass-radius diagram. The horizontal 
red band shows the observational constraint from the J16142230 mass measurement. Any 
equation of state whose line that does not intersect it is ruled out by this measurement. 

regime is that of Ref. [108], shown in Figure 4.1, in which the error bars of less than 5% in 
the 2 solar masses measurement rule out several equations of state. 

4.1.1 QCD surface tension and observable implications in astro- 
physics 

The surface tension between partonic and hadronic phases of Strong Interactions seems to be 
a key quantity in connecting microscopic modeling of cold and dense QCD with astrophysical 
applications: recently it has been put forward that a low enough surface tension could allow 
for possible observable consequences in different phenomena related to compact objects. 
It was shown for instance that deconfinement can happen during the early post-bounce 
accretion stage of a core collapse supernova event, which could result not only in a delayed 
explosion but also in an antineutrino burst that could provide a signal of the presence of 
quark matter in compact stars [109]. In Figure 4.2, the output of a core-collapse simulation 
shows the emission of the extra antineutrino burst more than 0.2 second after the standard 
neutrino signal of supernovae explosions. 

However, as was discussed in detail in Ref. [110, 111] (see also Ref. [112] ) those pos- 
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Figure 4.2: [extracted from Ref. [109]] Neutrino luminosities obtained in a core-collapse 
simulation assuming a hadron-quark transition. 

sibilities depend on the actual dynamics of phase conversion, more specifically on the time 
scales that emerge. In a first-order phase transition, as is expected to be the case in QCD 
at very low temperatures, the process is guided by bubble nucleation (usually slow) or spin- 
odal decomposition ("explosive" due to the vanishing barrier between phases), depending 
on how fast the system reaches the spinodal instability as compared to the nucleation rate 
[113]. Nucleation in relatively high-density, cold strongly interacting matter, with chemical 
potential of the order of the temperature, can also play an important role in the scenario 
proposed in Ref. [114], where a second (little) inflation at the time of the primordial quark- 
hadron transition could account for the dilution of an initially high ratio of baryon to photon 
numbers. Moreover, as shown in Figure 4.3, significantly different compact star structures 
are obtained if one considers the possibility of a layer of a quark-hadron mixed phase in the 
core of compact stars [115]. A key ingredient in all these scenarios is, of course, the surface 
tension, since it represents the price in energy one has to pay for the mere existence of an 
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interface between quark and hadron phases. If the cost in energy is too high, quark-hadron 
mixed phases will not be favourable and similarly nucleation time scales will be too long for 
these astrophysical phenomena to take place. 




J , I , I , I , I , i_ 

6 8 10 12 14 16 
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Figure 4.3: [extracted from Ref. [115]] Mass-radius relation obtained from equations of state 
using perturbative QCD results at high energies are compared to different observations. The 
left-most magenta band is the result in the presence of mixed phases, energetically viable 
only if the surface tension is low enough. 

However, the surface tension magnitude and the whole process of phase conversion via 
nucleation in a cold and dense environment are not known for strong interactions. In fact, 
the mapping of this sector of the (equilibrium) phase diagram is still in its infancy [116], let 
alone dynamical processes of phase conversion. 

Computing the surface tension for the QCD transitions with first-principle methods is a 
prohibitively complicated task. It is a genuine nonperturbative quantity, in the sense that it 
cannot be obtained by a naive expansion around one minimum of the effective action, and is 
also inaccessible by current first-principle lattice QCD simulations, since the Sign Problem 
at finite chemical potential brings about major technical difficulties for performing Monte 
Carlo lattice simulations. One has therefore to resort to effective models. Estimates for 
the surface tension between a quark phase and hadron matter were considered previously in 
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different contexts. In a study of the minimal interface between a color-flavor locked phase 
and nuclear matter in a first order transition, the authors of Ref. [117] use dimensional 
analysis and obtain S ~ 300 MeV/fm^ assuming that the transition occurs within a fermi 
in thickness. Taking into account the effects from charge screening and structured mixed 
phases, the authors of Ref. [118] provide estimates in the range of 50 — 150 MeV/fm^ but 
do not exclude smaller or larger values. These values of the surface tension are too high, 
rendering the nucleation of quark matter bubbles improbable in the context of compact 
objects. 

In what follows, we present our computation [101] of the surface tension for the cold and 
dense chiral phase transition within a well-established chiral model (namely, the linear sigma 
model with constituent quarks (LSMq)) and assuming the homogeneous nucleation scenario 
[119] (cf. also [113]). Results indicate that the surface tension in the cold and dense regime 
can be considerably smaller than previous estimates, possibly allowing for the interesting 
observable astrophysical phenomena described above. 

The remainder of this chapter is organized as follows. In Section 4.2, we delineate how, 
within the framework of homogeneous nucleation, it is possible to derive the nucleation 
parameters starting from an effective potential for the chiral order parameter V"efr(a'). Using 
the results for the effective potential V"cff(o") of the linear sigma model with quarks from 
Appendix A, we compute the nucleation parameters for three different cases, applying the 
method described in Section 4.2.1. The cold and dense result is presented in Section 4.3, 
while Section 4.4 contains the discussion on the competition between thermal effects and 
vacuum logarithmic corrections. Application of the results to the context of supernovae 
explosions is discussed in Section 4.5. Final remarks are found in Section 4.6. 

4.2 How to estimate the surface tension? 

Let us now consider the formation of droplets of quark matter at high density that happens 
via homogeneous nucleation [113]. Dynamically this process will occur either via thermal 
activation of droplets (thermal nucleation) or quantum nucleation. The physical setting 
we have in mind is the one that gives the best chances for thermal nucleation of quark 
droplets in cold hadronic matter found in "hot" protoneutron stars. As discussed in Ref. 
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[110, 111], and previously in Refs. [120, 121], that corresponds to temperatures of the 
order of 10 — 20 MeV [122]. At these temperatures, and in the presence of a barrier in the 
effective potential, thermal nucleation dominates over quantum nucleation. As soon as the 
barrier disappears, the spinodal instability is reached and the mechanism that takes over 
is the explosive spinodal decomposition. The range of temperatures under consideration is, 
then, high enough to allow for thermal nucleation (quantum nucleation being comparatively 
negligible [110, 111, 112, 123]) and low enough to justify the use of the zero-temperature 
effective potential. Temperatures of a couple of tens of MeV will not modify appreciably 
the equation of state, bringing corrections 0(T^/ /i^) ~ 1% for the typical values of chemical 
potential for the system under consideration. Nevertheless, as will be shown later, thermal 
corrections can be important for the process of nucleation. 

As stated above, with the aim of quantifying the surface tension associated with the cold 
and dense chiral phase transition in QCD, we will adopt the framework of homogeneous 
nucleation proposed by Langer [119]. It is based on the construction of a coarse-grained 
free energy functional J^[cr] for the configuration of the order parameter field a from a given 
effective potential Vcs as follows: 



Figure 4.4: Cartoon of the temperature versus chemical potential phase diagram with a first 
order critical line (solid) and the spinodal curves (dashed) defining the metastable region. 




(4.1) 



where spherical symmetry has been assumed. 



T4< spinodal lines 




Around the domain of temperature T and chemical potential /i in which a first-order phase 
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transition takes place, a metastable phase appears. More precisely, metastability manifests as 
the existence of (at least) two minima of the effective potential with a maximum in between 
inside the region defined by the two spinodal lines, as depicted in Fig. 4.4. Within the 
so-called metastable region, the coarse-grained free energy functional in Eq. (4.1) develops 
a nonuniform extremum configuration of the bubble type 0fe(r; i? = Rc), connecting the 
metastable minimum to the absolute one (for an illustration cf. Fig. 4.5). This unstable 
nonuniform extremum is called the critical bubble and is intimately related to the dynamics of 
phase conversion via bubble nucleation. Differently from the uniform configuration sitting at 
the absolute minimum, the critical bubble is unstable, corresponding actually to a maximum 
of the free-energy J-". 

If the system goes through an out-of-equilibrium stage in which the external parameters T 
and /i are modified fast enough as compared to the time scales of the microscopic dynamics, 
then the system may get trapped in the metastable phase. In this scenario, fluctuations 
(brought about by the interaction with a thermal or particle reservoir) may generate bubbles 
of the absolute minimum inside the metastable medium. These fluctuation bubbles will then 
evolve according to the competition between volume and surface terms, as dictated by the 
free-energy difference between the configurations with and without a bubble of radius R: 



whose general form is sketched in Fig. 4.6. The volume term encodes the gain in energy 
resulting from the conversion to the true equilibrium state and is proportional to the pressure 
difference AP > between the phases. On the other hand, the surface term accounts for 
the energy spent in the construction of a surface between the phases, being given essentially 
by the surface tension S. 

As illustrated in Fig. 4.6, the free-energy difference AJ-" shows a maximum associated 
with the critical bubble {R = Rc) and fluctuation-generated configurations containing bub- 
bles with radii R < Rc will tend to the uniform metastable vacuum in order to minimize the 
energy. For radii bigger than the critical radius, the minimization of energy implies that the 
droplet will grow and eventually complete the phase conversion. The critical radius goes to 
zero at the spinodal curves, where the metastable false vacuum becomes unstable (the bar- 
rier disappears) and the phase conversion occurs explosively via the spinodal decomposition 




(4.2) 
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Figure 4.6: Sketch of the difference 



Figure 4.5: Example of bubble pro- 



of coarse-grained free energy be- 



file. 



tween the configurations with and 



without a bubble of radius R. 



process. 



Finally, from Eq. (4.2) it is clear that the surface tension can be obtained via the 
knowledge of the effective potential V^s and the critical bubble profile 0;,(r; R = Rc). Since 
the critical bubble itself is an unstable extremum of the coarse-grained action J-" in Eq. 
(4.1), it is ultimately obtained from the effective potential Vcs as well. In the thin wall 
approximation (valid when the critical bubble has a thin interface region as compared to its 
radius, i.e. ^ <^ Rc, cf. Fig. 4.5) for a quartic effective potential, the solution 0b(r; R = Rc) 
can be written analytically in terms of the Taylor coefficients of the potential. 

4.2.1 Extracting nucleation parameters from the effective poten- 
tial 

Since our framework is an effective model, we can only aim for reasonable estimates and 
functional behavior, not numerical precision. Therefore, it is convenient to work with ap- 
proximate analytic relations by fitting the relevant region of the effective potential by a 
quartic polynomial and in the thin-wall limit approximation for bubble nucleation. Follow- 
ing Refs. [82, 83], we can express the effective potential over the range between the critical 
chemical potential, fj,c, and the spinodal, figp, in the familiar Landau-Ginzburg form 



4 




(4.3) 



n=0 
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Although this approximation is obviously incapable of reproducing all three minima of Vcs, 
this polynomial form is found to provide a good quantitative description of this function in 
the region of interest for nucleation, i.e. where the minima for the symmetric and broken 
phases, as well as the barrier between them, are located. 

A quartic potential such as Eq. (4.3) can always be rewritten in the form 



V{ip) = a {if'-a'Y+jip 



with the coefficients above defined as follows: 



(4.4) 
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(4.5) 
(4.6) 

(4.7) 
(4.8) 



The new potential V{(p) reproduces the original Vcs{4>) up to a shift in the zero of energy. 
We are interested in the effective potential only between fic and figp- At fic, we will have two 
distinct minima of equal depth. This clearly corresponds to the choice j = in Eq. (4.4) 
so that V has minima at (/9 = ±0 and a maximum aX (f = 0. The minimum aX (f = —a 
and the maximum move closer together as the chemical potential is shifted and merge at 
Hsp- Thus, the spinodal requires jjoic? = —8/3^/3 in Eq. (4.4). The parameter jjac? falls 
roughly linearly from 0, at /i = /ic, to — 8/3v^ at the spinodal. 

The explicit form of the critical bubble in the thin- wall limit is then given by [124] 



1 — tanh 



Rr 



(4.9) 



where (ff is the new false vacuum, Rc is the radius of the critical bubble, and ^ = 2/m, 
with = V"{ipf), is a measure of the wall thickness. The thin- wall limit corresponds to 
$,/Rc <^ 1, which can be rewritten as (3|j|/8ao^) ^ 1 [124]. Nevertheless, it was shown in 
[82, 125] for the case of zero density and finite temperature that the thin-wall limit becomes 
very imprecise as one approaches the spinodal (This is actually a very general feature of this 
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description [113]). In this vein, also as remarked above, the analysis presented below is to 
be regarded as semi-quantitative, and we aim for estimates. 

In terms of the parameters a, a, and j defined above, we find [82, 83] 



^tj ^ ±a-^ , (4.10) 

1/2 



in the thin- wall limit. The surface tension, S, is given by 



oo 



(4.11) 



and the critical radius is obtained from Rc = (2S/AV), where AV = V{(j)f) — V{(j)t) ~ 2a|j|. 
Finally, the free energy of a critical bubble is given by Fi, = (47rS/3)-R^, and from knowledge 
of Fh one can evaluate the nucleation rate T ~ e"^''/^. In calculating thin-wall properties, 
we shall use the approximate forms for (pt, (pj, S, and AV^ for all values of the potential 
parameters. 



4.3 Results for nucleation in the cold and dense LSM 

In what follows, we consider the LSM with quarks in the absence of vacuum corrections. 
The effective potential up to 1-loop order is then given by Eq. (A.5) with = n^^d- 

This 

corresponds to the standard case, adopted frequently in the literature (cf. e.g. Ref. [81]). 

Using the method described in the previous section, we characterize quantitatively the 
nucleation process predicted within this case, calculating different nucleation parameters. 
In this section, we present results for both metastable regions, above and below the critical 
chemical potential /icrit = 305.03 MeV, which correspond respectively to the nucleation of 
quark droplets in a hadronic environment and the formation of hadronic bubbles in a partonic 
medium. 

In Fig. 4.7, the critical radius, namely the radius of the critical bubble, is displayed. 
Any bubble created in the system via external influences or thermal fluctuations will either 
grow or shrink unless its radius equals the critical value, which settles the threshold between 
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Figure 4.7: Critical radius as a function of the quark chemical potential in the two metastable 
regions (between the spinodals). The arrows indicate the critical chemical potential. 

suppressed and favored bubbles. For radii bigger than the critical radius, the minimization 
of energy implies that the droplet will grow and eventually complete the phase conversion. 
The critical radius goes to zero at the spinodal curves, where the metastable false vacuum 
becomes unstable (the barrier disappears) and the phase conversion occurs explosively via 
the spinodal decomposition process^. On the critical line the vacua become degenerate so 
that both of them are stable yielding no nucleation, or equivalently, a divergent critical 
radius. In the case of the chiral transition in the LSM with quarks, we obtain that the 
critical radius is i?c < 10 fm over about ~ 90% of the metastable regions. 

The correlation length ^ which provides a measure of the size of the bubble wall is plotted 
in Fig. 4.8, in units of the critical radius. As discussed in the previous section, the thin-wall 
approximation relies on the assumption that the ratio between the surface correlation length 
^ and the radius of the critical bubble is small. Fig. 4.8 shows clearly that this is a reasonable 
condition in the vicinity of the critical line, away from the spinodals. 

The surface tension is a key parameter in quantifying the nucleation process in a given 
medium since it measures the amount of energy per unit of area that is spent in the construc- 
tion of a surface between the phases. In Fig. 4.9 we show our result for the surface tension, 

^In the thin-wall approximation, which is a poor description of the regions close to the spinodals, the 
critical radius does not vanish, only becomes very small. The same is true for other quantities, as the surface 
tension. 
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Figure 4.8: Ratio between the correlation length ^ and the critical radius Rc as a function 
of the quark chemical potential in the two metastable regions. 

E, for the chiral transition in the cold and dense LSM with quarks. The surface tension 
assumes values between ~ 4 and ~ 13 MeV/fm^, being, throughout the whole metastable 
region, of the order of magnitude that renders the formation of quark matter viable during a 
supernova explosion, according to Ref. [109]. The biggest values occur near criticality, since 
this domain is characterized by large barriers and a small free energy difference between 
the true and false vacua. It should be noticed that those values are well bellow previous 
estimates of this quantity, as will be discussed in the conclusion of this chapter. 




Figure 4.9: Surface tension as a function of the quark chemical potential in the two 
metastable regions. 
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We have also estimated the nucleation rate for homogeneous nucleation [119, 126] as F ~ 
Tj e~^*>/Tf ^ ^Y^QYe Fi, is the free energy of the critical bubble configuration and Tf = 30 MeV 
is an ad hoc temperature. Applying the expression above for the nucleation rate in this 
case means that we are neglecting the temperature dependence of the critical-bubble free 
energy in the exponent. The results are shown in Fig. 4.10 for both metastable regions. The 
nucleation rate falls abruptly as the chemical potential approaches its critical value. 

It is interesting to point out that the difference in the size of the metastable regions above 
and below criticality might play an important role in determining whether nucleation is a 
viable process of phase conversion for a given system, as compared to spinodal decomposition 
scenarios or even if the lifetime of the system represents enough time for nucleation to take 
place. It is clear in Fig. 4.10, for instance, that the domain in chemical potential for which 
the nucleation rate assumes sizable values is larger at the metastable region corresponding 
to bubble formation of hadronic matter in a partonic medium. 




Figure 4.10: Nucleation rate as a function of the quark chemical potential in the two 
metastable regions. 
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4.4 Effects from vacuum terms versus thermal correc- 
tions 

Having presented in the section above the important parameters for nucleation and discussed 
their results for the zero-temperature chiral transition in the LSM, let us now analyze the 
influence of quantum logarithmic terms in the vacuum effective potential as well as thermal 
corrections. 

In this section we focus on the metastable region above the critical chemical potential, 
aiming for thermal nucleation of quark matter droplets in compact objects and especially in 
supernovae explosions. In this vein, we also keep relatively low temperatures when including 
thermal corrections, up to T = 30 MeV. The systematic inclusion of thermal corrections 
provides a measure of the validity of the zero-temperature approximation for the thermo- 
dynamics and phase structure in low-temperature scenarios such as those found in compact 
objects. We find that results for the nucleation parameters up to T = 10 MeV present 
variations within ~ 10% of the zero-temperature values, displayed in the previous section. 

In Figures 4.11-4.14, the role played by temperature corrections and quantum vacuum 
effects is presented. More specifically, we compare results for the three following situations: 
(a) classical vacuum effective potential with thermal corrections for T = 30 MeV, (b) quan- 
tum vacuum effective potential at zero temperature, (c) classical vacuum effective potential 
at zero temperature, i.e. the results shown in the last subsection. Each case is associated 
with a different critical chemical potential (indicated by arrows in the plots; the normaliza- 
tion /icrit.o = 305.03 MeV corresponds to the critical chemical potential of the case (c), of 
the last subsection) and also a different spinodal chemical potential (denoted by the vertical, 
dashed lines). 

Concerning the modification of the metastable regions, we show that the quantum vacuum 
corrections increase significantly the domain in chemical potential in which metastability 
persists; the metastable region in the presence of quantum vacuum terms is ~ 40% bigger 
than the one associated with the classical vacuum effective potential. This feature indicates 
that such corrections might generate large differences in the dynamic evolution of the system, 
even though the absolute value of the critical chemical potential itself shifts only ~ 2%. 
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Figure 4.11: Density in units of the nuclear saturation density uq = 0.16 fm~^ as a function 
of the quark chemical potential in the cold and dense LSM (curve (c)), including vacuum 
terms (curve (6)) and thermal corrections (curve (a)). 

During the dynamical process of phase conversion, the density is increased (e.g. via 
gravitational pressure in the supernova explosion scenario) and therefore the free energy 
required for including a new quark in the system also augments. The mapping between 
density and quark chemical potential within the cases we consider is plotted in Figure 4.11. 
The discontinuities signal the first-order phase transition that restores chiral symmetry at 
high energies. As the evolving system reaches a high enough density, it penetrates the 
metastable region above the critical chemical potential, where higher densities are associated 
with higher super compression. 

It is clear that the values of density predicted below the transition are irrealistically 
low. This is a consequence of the fact that the LSMq fails to describe nuclear matter. 
Being an effetive theory constructed with the purpose of describing the chiral symmetry 
features of strongly interacting matter, the LSMq does not predict the low-energy nuclear 
matter properties nor the gas-liquid transition that forms it. However, since our ultimate 
interest is investigating the formation of quark matter in astrophysical processes involving 
ultra-compact objects, the framework most suited is exactly one that describes the high- 
energy transitions of strong interactions, namely chiral restoration and/or deconfinement. 
We choose therefore the LSMq and analyze the nucleation of chirally-symmetric droplets in 
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a strongly interacting medium. 

In order to interpret the results for the nucleation parameters, one should keep in mind 
the variation of the metastable regions and its implication on the value of chemical potential 
that needs to be reached above the respective criticality, i.e. the amount of supercompression. 




Figure 4.12: Critical radius (left) and correlation length (right) in the nucleation region of 
quark droplets as a function of the quark chemical potential in the cold and dense LSM 
(curve (c)), including vacuum terms (curve (6)) and thermal corrections (curve (a)). 

In Figure 4.12, the left panel shows the variation of the radius of the critical bubble, 
while the right one displays the correlation length (or equivalently the approximate width of 
the bubble wall) for the three cases. A direct consequence of the shifting of the metastable 
domain is the fact that, for fixed chemical potential, the critical radius decreases signifi- 
cantly both when one includes quantum vacuum corrections (curve (6)) or in the presence of 
thermal effects (curve (a)). However, the amount of supercompression needed for reaching 
a certain viable value of critical radius clearly increases when one includes quantum correc- 
tions and decreases in the case with thermal corrections. Therefore, there is a nontrivial 
competition between those effects: while quark matter nucleation is disfavored by the inclu- 
sion of quantum vacuum corrections, it is favored by thermal corrections to the nucleation 
parameters. 

The competition between vacuum and thermal corrections is more explicitly shown by 
our findings for the surface tension, displayed in Figure 4.13. As compared to the zero- 
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Figure 4.13: Surface tension in the nucleation region of quark droplets as a function of the 
quark chemical potential in the cold and dense LSM (curve (c)), including vacuum terms 
(curve (b)) and thermal corrections (curve (a)). 

temperature, classical computation (curve (c)), the establishment of a interface between 
quark matter droplets and the hadronic medium costs more with the incorporation of quan- 
tum corrections in the effective potential. On the other hand, finite temperature terms tend 
to push the surface tension down, reaching a minimum of only ~ 3 MeV/fm^, facilitating 
the surface formation. Therefore, if the environment in the core of a compact object is 
cold enough and if the effective theory including quantum corrections in the vacuum is the 
most suited for describing the chiral transition, then the surface tension can be as high as 
~ 17.5 MeV/fm^ near criticality rendering quark matter nucleation much slower and even 
an improbable phenomenon. 

Analogous results for the nucleation rate F can be seen in Figure 4.14. Once again, the 
amount of super compression associated with a given chemical potential should be kept in 
mind. The nucleation rate falls abruptly for chemical potentials approaching the respective 
critical values. In the vicinity of the spinodals, the nucleation rate is dominated by the 
pre-exponential factor and reaches sizable values, around ~ 0.1 MeV/fm^. 
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Figure 4.14: Nucleation rate in the nucleation region of quark droplets as a function of the 
quark chemical potential in the cold and dense LSM (curve (c)), including vacuum terms 
(curve (b)) and thermal corrections (curve (a)). 

4.5 Consequences for the quark-matter-induced super- 
nova explosion scenario 

To contribute to the investigation of the possibility of nucleating quark matter droplets 
during the early post-bounce stage of core collapse supernovae, we follow Ref. [110, 111] and 
define the nucleation time as being the time it takes for the nucleation of a single critical 
bubble inside a volume of Ikm^, which is typical of the core of a protoneutron star, i.e. 



Fig. 4.15 shows this quantity as a function of the normalized chemical potential for the 
three scenarios considered in this analysis. The relevant time scale to compare is the time 
interval the system takes from the critical chemical potential to the spinodal during the 
supernova event, if it ever reaches such high densities in practice. Implicit in the definition 
above is the approximation of constant density and temperature over the core, which is fine 
for an estimate, since density profiles are quite flat in this region of the star. The typical 
time scale for the early post-bounce phase is of the order of a fraction of a second, so that 
the time within the metastable region is smaller, in the ballpark of milliseconds. 




(4.13) 
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Figure 4.15: Nucleation time defined for a typical volume of 1 km^ in the cold and dense 
LSM (curve (c)), including vacuum terms (curve (b)) and thermal corrections (curve (a)). 

Keeping these simplifications in mind, together with the caveat that the LSMq does not 
describe nuclear matter besides the chiral properties of strong interactions, we show in Fig. 
4.15 that vacuum corrections tend to increase the density depth required for efficient nu- 
cleation whereas thermal corrections push in the opposite direction, consistently with the 
results discussed in the previous subsection. Our results for the nucleation time, together 
with those for the surface tension, tend to favor the best setting for nucleation of quark 
matter in the supernova explosion scenario considered in Ref. [110, 111], especially when 
thermal corrections with physical temperatures are included. To a great measure, this hap- 
pens because the nucleation time, and the whole process of phase conversion, depends very 
strongly on the surface tension, since it enters cubed in the Boltzmann exponential of the 
rate F. 



4.6 Final remarks 

In this chapter we have explored the interplay between the dynamics of cold and dense phase 
transitions of Strong Interactions and astrophysical phenomena involving ultra-compact ob- 
jects. Since these environments contain QCD matter under extreme pressures, much beyond 
any current laboratory can achieve, their understanding requires a theoretical description of 
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Strongly interacting cold and dense media. For the same reason, with the development of 
observational techniques, astrophysical observations may become an invaluable experimen- 
tal channel to assess the low-temperature high-density domain of the QCD phase diagram 
as realized in Nature. Taking thorough advantage of this channel demands quite often the 
description of dynamical critical phenomena within Strong Interactions. 

As emphasized throughout the chapter, a striking example of this connection of out-of- 
equilibrium QCD phase transitions with observable astrophysical phenomena is the surface 
tension between the hadronic and the chirally restored deconfined phases in the cold and 
dense regime. If its value is small enough, quark-matter formation could occur during core- 
collapse supernovae explosions, providing an alternative dynamics and even an observable 
signal of the QCD phase transition within compact objects. Moreover, a reasonably small 
surface tension between partonic and hadronic phases could also contribute to allow for 
different compact star structures including mixed phases. 

Modelling the chiral phase transition via the linear sigma model with quarks, we have 
studied homogeneous nucleation in a framework that allowed for analytic calculations, given 
by a potential fit by a quartic polynomial and the thin-wall approximation. All the relevant 
quantities were computed as functions of the chemical potential (or the baryonic density), 
the key function being the surface tension. 

In particular, we have showed that the LSMq predicts a surface tension of S ~ 5- 
15 MeV/fm^, much lower than previous estimates (e.g. [117] and [118] ) and possibly allowing 
for the interesting observable phenomena described above. Including temperature effects and 
vacuum logarithmic corrections, we find a clear competition between these ingredients in 
characterizing the dynamics of the chiral phase conversion, so that if the temperature is low 
enough the consistent inclusion of vacuum corrections could help preventing the nucleation 
of quark matter during the early post-bounce stage of core collapse supernovae. 

As discussed in Subsection A. 4, effects from interactions between the sigma field and the 
quarks that come about at two-loop order could in principle contribute to this competition 
as a third sizable modification and should be investigated as well. 

The linear sigma model, however, does not contain essential ingredients to describe nu- 
clear matter, e.g. it does not reproduce features such as the saturation density and the 
binding energy. Therefore, the results obtained in this chapter should be considered with 
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caution when applied to compact stars or the early universe. It is an effective theory for a 
first-order chiral phase transition in cold and dense strongly interacting matter, and allows 
for a clean calculation of the physical quantities that are relevant for homogeneous nucleation 
in the process of phase conversion. In the spirit of an effective model description, our results 
should be viewed as estimates that indicate that the surface tension is reasonably low and 
falls with baryon density, as one increases the supercompression. First-principle calculations 
in QCD in this domain are probably out of reach in the near future. Therefore, estimates 
within other effective models would be very welcome. 



Chapter 5 




Thermodynamics of QCD 
in a magnetic background 



In this chapter, a hot QCD medium exposed to a constant and uniform external magnetic 
background is considered. The magnetic field becomes therefore a new axis in the phase 
diagram of Strong Interactions to be explored, giving rise, for example, to the temperature 
versus magnetic field plane, which we address in particular. 

Some interesting phenomena have been predicted and explored recently regarding Strongly 
interacting matter under extremely intense magnetic fields, but we will be interested here 
in posing a well-defined question within the fundamental gauge theory: what happens to 
the thermodynamics of QCD in the presence of such an intense magnetic background? At 
sufficiently high energies, the answer to this question should be given by perturbative QCD 
in a nonperturbative magnetic background. That is the approach we adopt in this on-going 
project [127] to compute the QCD pressure and the thermal self-energies, on the path of 
identifying the adequate quasiparticles in this domain of the QCD phase diagram. 
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5.1 Motivation and framework 

The area of in-medium Strong Interactions under extreme magnetic fields has gone through 
a fast development due to the combination of three features realized recently: experimental 
relevance, capability of being investigated on the lattice and a rich phenomenology. 

Hot QCD matter exposed to sufficiently high magnetic fields is currently being generated 
in laboratory experiments. The setup of ultra-relativistic heavy ion collisions is such that 
extremely intense magnetic fields are produced in peripheral collisions, with a nonzero impact 
parameter. A beam of heavy ions like Au accelerated to velocities very close to that of the 
light creates an enormous electric current. Parallel, non-coinciding beams going in opposite 
directions induce, through a basic classical electrodynamics process, a maximum magnetic 
field at the center, exactly where the medium forms in such collisions. At RHIC the intensity 
of the magnetic fields created may reach 6m^ ^ 10^-10^ MeV^ ~ 10^^ Gauss [128]. 

There have been many developments in the direction of using this extreme field to pro- 
duce signatures to assess the properties of the new state of matter created in these exper- 
iments. For instance, the Chiral Magnetic Effect [129] takes advantage of the capability 
of magnetic fields of spatially separating different electric charges to render measurable the 
possible formation of sphaleron-induced CP-odd domains in the medium. Currently under 
investigation are also the modifications of the QCD phase structure in the presence of this 
new external parameter. Sufficiently large magnetic fields could allow for drastic phenomena, 
from affecting the chiral and the deconfinement phase transitions [130, 131, 132, 133, 134] to 
transforming vacuum in a superconducting medium via p-meson condensation [15]. Although 
complemented by lattice results [135, 136, 137, 138], most of these analyses were done within 
effective frameworks, as one is usually obliged to do when treating the many ingredients that 
are relevant to approach experimental conditions or describing phase transitions. 

It is important, however, to address the high-temperature, high-magnetic field regime of 
the QCD phase diagram (semi-) analytically from the point of view of the fundamental gauge 
theory. Since the T—B plane of the phase diagram is amenable to lattice QCD simulations, it 
represents a new open channel for comparison between theory and numerical nonperturbative 
approaches. Our aim in this on-going project [127] is to make proper use of this opportunity 
by providing perturbative QCD predictions for the thermodynamics of a hot environment 
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immersed in an intense magnetic field. In this way, one gains further analytical understanding 
of lattice results on one side and tests the validity of the perturbative expansion in this new 
regime on the other. Furthermore, this first-principle investigation within the fundamental 
gauge theory shall also furnish guidance for effective-model building to go beyond the domain 
of applicability of perturbation theory. 

Interesting phenomena in Strongly interacting media arise only when magnetic fields 
are sufficiently large as compared to relevant scales such as m^. Therefore, the specific 
physical system we address in this chapter is hot QCD matter under an intense constant and 
uniform magnetic background at very high energies. Asymptotic freedom then suggests the 
applicability of perturbation theory, while the (classical) magnetic sector must be treated 
nonperturbatively. In this vein, the framework adopted is described by the QCD Lagrangian 
in the presence of a classical electromagnetic vector potential A (V x A = B = Bz; Ad = 
(0,1)): 

£qcd = -\TT[G'^'^G,,] + J2^f[z$-mf-qf4ci]'^f-9j2^f0^f, (5.1) 

/ / 

where G^ = G^A" is the gluon field (A" are the Gell-Mann matrices) and G^j^ = — 
duG'^^ + gf^^'^G'lj^Gl is the gluon field tensor in a standard SU{Nc) gauge theory. Here, the 
quark-gluon interaction is treated perturbatively up to 0{g^), while the magnetic coupling 
is considered as a nonperturbative dressing of the fermion propagator. In this vein, our 
description is that of a medium of magnetically-dressed quarks perturbed by the coupling 
with gluons. 

In what follows, the framework is concretely built by the derivation of the exact fermion 
propagator in a magnetic background in Section 5.2 and then used to compute the one- 
loop gluon selfenergy in the vacuum (Section 5.3) and the QCD thermodynamical potential 
to two-loop order (Section 5.4). Remarks and perspectives of this work in progress are 
summarized in the end. 

5.2 Exact Fermion propagator in a magnetic background 

The most crucial nontrivial ingredient of our approach to QCD in an intense magnetic 
background is the nonperturbative description of the propagator of the fermions dressed by 
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magnetic effects. This propagator corresponds to the Green's function associated with the 
Dirac equation in the presence of a classical electromagnetic potential A'^y{x) = {0,A{x)), 
such that the magnetic field is produced B = Bz = V x A: 

+ ^q4M)) - ^/] So{x, y) = (27r)^5^(x - y) , (5.2) 

or, in general operator form, using Il'^{p,x) = {p^ — qA'^^{x)) and {x\So\y) = So{x,y), 

[]/L{p,x)~mf]So = l, (5.3) 

where it is clear that the presence of a space-dependent classical potential A^i (essential for 
generating the constant and uniform magnetic field) renders the Dirac operator a mixture 
of momentum and coordinate operators. The matrices in both of these representations are 
therefore nondiagonal and infinite-dimensional, so that their diagonalization is an intricate 
task, as well as computing their determinant. 

One may still formally write the solution of this equation as 



So = ^ = J, ^[yL{p,x) + mf]. (5.4) 



lZt(p, £) — m/ [lZl(p, x)]^ — 



In practice, one needs to transform this formal solution into an explicit matrix element 
(e.g. in coordinate space: {x\So\y) = So{x,y)) which can be used directly in diagrammatic 
computations. Originally, this task was accomplished by Schwinger [139], who introduced 
the now called Schwinger formalism to find an exact explicit solution which we shall sketch 
in what follows. 

First, it is convenient to re- write the inverse operator appearing in Eq. (5.4) as an integral 
of a unitary operator, getting the following Schwinger integral form for the propagator: 

POO 

Sq = i j (is exp {— is [lZl(p, x)]^ — mj] } [lZl(p, x) + m/] . (5.5) 
Jo 

We may therefore re-interpret the unitary operator U{s) = exp[— is [lZt(p, x)]^] as an evolution 
operator in an auxiliary time s governed by an auxiliary hamiltonian l-L = li(p, x)]^. Within 
this picture, the propagator becomes an average in the time s of the operator [lZl(p, x) +Tnf] 
evolved according to the auxiliary hamiltonian l-i: 



poo 

So = i ds e'^'^^f U{s)[t{p,x)+mf] 
Jo 



(5.6) 
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or, in coordinate space, 

POO 

So{x,y) = I dse'''^Ux\U{s)mP,i)+mf]\y) (5.7) 
Jo 

= i rfs e^"""/ = x(s)|l!(p,x)|?/ = x(0)) +m/(x = x(s)|?/ = x(0))}(.5.8) 

Solving the evolution in the auxiliary time s is therefore equivalent to solving the operatorial 
Dirac equation in the presence of a classical external field A'^^i^x). 

The detailed computation of the evolved matrix elements (x = x{s)\y = x(0)) and 
{x = x{s)\^{p,x)\y = x(0)) can be found in Ref. [139] and the results are in Eqs. (3.14) and 
(3.20) of this reference. They have the form: 

{x = x{s)\yL{p, x)\y = x{0)) = f^{x-y){x = x{s)\y = x{0)), 

{x = x{s)\y = x{0)) = f2{x -y)<!>{x,y) , (5.9) 

where ^{x,y) is the only contribution that breaks translational invariance, being known as 
the Schwinger phase: 

rdx'^A^,{x') . (5.10) 



$(x,2/) = exp 



Ultimately, this phase is responsible for the non-invariant form of the fermion propagator 
dressed by magnetic effects: 

So{x,y) = ^x,y)So{x-y) (5.11) 

POO 

= ^x,y)t dse'^"'" f2ix-y) [Mx-y)+mf]. (5.12) 

^0 

Notice that this apparent breaking of translational invariance is not a physical property (all 
points in space-time are completely equivalent), since the propagator itself is not physical, 
being gauge dependent as well. In fact, this phase evaluates to 1 in a specific nontrivial 
gauge choice [140]. 

Choosing the electromagnetic gauge as Ad = {0,0, Bx,0), the final form for the exact 
fermion propagator in a uniform and constant magnetic background {B = B z) is [139] (cf. 
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also [140, 141]): 



exp <,'-^K- + y) \ So{x - y) 



exp<j|[x^-y^]A-*(a: + y) 



iq 



- y^^]A'^^\x + y) 



exp , 2 L- ^ 



X 



d'^P 
{2nY 



Q-iP-i^-y) ^ 



(5.13) 
(5.14) 



POO 




/ ds exp 


is ^ 


Jo 





2 2 tan(gfis) 



X 



: { [Pl ■ 7l + "^/] [1 + tan(g55)] - ■ 7t [1 + tan^ (g5s)] } 



(5.15) 

where we have used a compact notation for the transverse and longitudinal quantities: = 

(pi,P2), 7t = (7\7^)> PL = (pcPs), and 7l = (7°,7^)- 

This propagator presents consistently some desired properties. First, the free Dirac 
propagator is recovered in the absence of magnetic fields. Indeed, this is obtained directly 
from limB_5.o tan(gi?s) = qBs — t- and the basic property of the Schwinger integral (with a 
regularization of the exponential implicitly assumed): 



ds exp [isF] 



F 



(5.16) 



Moreover, it can be shown that the poles in Eq. (5.15) occur at (and only at) the 
energies for which the Landau level dispersion relation is satisfied: — + + '^nB {n = 
0, 1, 2, ■ ■ ■ ). The demonstration [140] (cf. also [141]) goes through rewriting the Schwinger 
integral into a sum over Landau levels. Using {z = [— exp{—2iqBs)]) 



itsin{qBs) 



1 - 2- 



2)-("+^)exp 



xz 



n=0 



X Z" 



(5.17) 
(5.18) 



where -C/"(x) are the generalized Laguerre polynomials, the exact result for the propagator. 
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Eq. (5.15), becomes a sum over Landau levels^ 

So{P) = iexp 
with (defs.: L„ = L° and L^i = 0) 



Pt 

qB 



^ Dn{qB,P) 

' pl-mj- 2nqB 



(5.19) 



n=0 



DniqB, P) = (ki ■ 7L + m) (1 + i^'Y)Ln 2 



+4(kT ■ 7t)^^_i (2^) • 



(1 - tYY)Ln 



(5.20) 



5.2.1 Lowest Landau Level approximation 

This expansion over Landau levels is particularly convenient when one is interested in phe- 
nomena at very intense magnetic fields, because large fields should guarantee the fast con- 
vergence of this series, allowing for approximations that take into account only the first 
levels. 

When all mass and momentum scales are indeed much smaller than the considered mag- 
netic field, one may in principle take the Lowest Landau Level (LLL) approximation^, which 
is equivalent to substituting the sum over Landau levels for its value at n = 0. Up to 
regularization differences, the LLL propagator reads: 



exp 



iq 



x'' - y^']A''^\x + y) 



d'^P 



e ^'^■^^ ^^zexp 



\qB\J PL-lL-rUf 
(5.21) 



which is equivalent to the result used in Ref. [142], obtained via a method of constructing 
the projectors for the different Landau levels from the exact solution of the Dirac equation 



^This is nevertheless a subtle transformation that relies on assuming a nonzero B (so that, in the final 

expression, the free Dirac result cannot be recovered) and regularization of the Schwinger integral (and of 

the sums appearing in the relation in Eq. (5.18)). 

^From the form of the propagator in Eq. (5.20), and the behavior of the Laguerrc polynomials when their 

argument vanishes, the connection between the LLL approximation and the large B limit becomes clear. 

This equivalence between the leading-order 1/B expansion and the LLL approximation is not always true 

(cf. the computation of the free pressure in Subsection 5.4.1). 
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[143]: 

cLLL/ N f dpodp2dp3 \ ■ ( \ ^ ■ ( ^ , • ^ ^\ 

-^o K^^y) = j (-27r)'^^i y—ew[-m{xo-yo) + tp2{x2-y2) + tP3{x3-y3)jx 

with Vq = diag(l,0, 1,0) = (1 + i^^^'^)/2 being the projector over the spin states which 
are present physically in the LLL. The peculiarity of the LLL approximation should be 
noted: this propagator is still a 4 x 4 matrix, but it describes only two physical propagating 
modes. Two eigenvalues of Sq^^ go to zero in the vicinity of the n = Landau level 
pole [pq = mj + pg), posing problems for the computation of the free pressure in this 
approximation (cf. Section 5.4.1). 



5.3 Extreme magnetic fields and gauge boson selfen- 
ergy 

As a first step towards understanding the quasiparticle content of a Strongly interacting 
system exposed to an intense constant and uniform magnetic field, let us compute the self- 
energy of the gluon corrected by one loop of dressed quarks in the vacuum. 

This is a simple calculation that shall illustrate basic features of perturbation theory over 
the nonperturbatively dressed magnetic fermions at work, such as the effective reduction to 
(1 + 1) dimensions for very intense backgrounds. In Ref. [142] the analogous QED case 
in the LLL approximation was considered and complemented by an approximate spectral 
construction to show the appearance of a peculiar mode of the corrected photon propagation, 
similar to the zero sound mode in Fermi liquid theory [144]. 

The actual computation for the one-loop gluon selfenergy is very similar; the only differ- 
ence resting on the extra color tensor structure in the vertex {e'-f'^ in QED is mapped into 
Q'-f'^X"' in QCD, where A" are the Gell-Mann matrices). Indeed, in the general situation, in 
which translational invariance is not assumed, the selfenergies for the photon ^photoni^^l) 
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and the gluon Il^'^{k,q) may be written as: 



d^xd^ye^^-^+^^-^e^Tr [YSo{x, y)YSo{y, x)] , (5.23) 
d^xd'^ye'''-''+''''-yiW{x, y) (5.24) 
rf^xrfV^-"+*''-^^7'Tr [7^A'^S'o(x, y)^" X'' S^{y , x)] , (5.25) 



so that the Ilphoton maps into the gluon selfenergy if the following substitution is imple- 
mented^: (— (^^TicoiorA^A" = g'^{N^ — l)/2. Therefore, the presentation of this section 
follows closely that of Ref. [142]. 

In accordance with the framework adopted, Eq. (5.25) is exactly the perturbative ex- 
pression for the gluon selfenergy at one loop, but with the fermion propagator being dressed 
by (nonperturbative) magnetic effects. 

As discussed above, we adopt the LLL approximation, valid for extremely intense fields, 
much larger than any other mass or momentum scale in the problem. Taking the LLL quark 
propagator, Eq. (5.22), into the definition of the gluon selfenergy (Eq. (5.25)), one gets: 

dpodqodp2dq2dp3dq3 



m^"'{k,k') 



2 



(5.26) 



where the phase ip^y carries the x and y-integrals from the Fourier transform. 



■'xy 



d^xdSe'''"'^'^'-ye^Y>{ - «(Po - go)(xo - yo) + i{P2 - q2){x2 - 2/2) + ^(ps - q3){x3 - y^) 



-exp|--[(x,--j ^[y^-^J U' 

which is solved in Appendix E, yielding Eq. (E.9)). The delta functions from this phase are 
then used to perform the qi {i = 0, 2, 3) integrals, furnishing momentum conservation for the 



eB 



(5.27) 



•^In this section, we shall consider the contribution to the selfenergy coming from only one fermion with 
electric charge e in both QED and QCD cases. Nevertheless, in the full QCD gluon selfenergy quarks carry 
color and flavor indices and their charge is flavor-dependent, so that Ilg'^j-, = ^cJ2f ^^"{Qf)- 
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external legs in directions i = 0,2, 3: 

1 



iW{k,k') 



-9 
exp 
Tr 



iV2 



(2^)3^(0,2,3)^^/^^) 



dpodp2dp3 
(27r)3 



AeB L 



k'l + kl + 2i{2p2 + k2){k[ + ki] 
1 ., _ 1 



2ik2y 



(5.28) 



The last integral that does not involve the trace in Eq. (5.28) is then the integral over p2: 

j^exp[-^[2i(2p2 + k2)ik[ + k,)]] = eBSik[ + k,), (5.29) 

where we have absorbed k2 via a shift before integrating. 

Finally, the expression for the gluon selfenergy showing explicitly the momentum conser- 
vation of the external legs reads: 



iW{k,k') 



pB 

{27r)^S\k' + k)—exp\ - 



Tr 



2ti 
1 



2eB 



(27r)2 



(5.30) 



VL-lL-mf (pl + ■ IL - ruf. 

It is interesting to note that the effect of an extremely intense magnetic field on the 
gluon (or photon) propagation is two-fold: (i) the B dependence is restricted to an overall 
decaying exponential factor on the transverse momentum and (ii) the nontrivial dynamical 
contribution is exactly the vacuum gluon selfenergy in (1 + 1) dimensions. 

To verify this latter property it is necessary to examine further the projected trace ap- 
pearing above, i.e. 

1 ., _ 1 



Ta = Tr 



(5.31) 



PL-lL-mj (pi + ki) ■ 7i - m/. 

where Vq = (1 +27^7^) /2 is the projector over the degrees of freedom actually present in the 
LLL (only one spin state). Its properties 7^'Po = 'Pol'^ /^^'^ g g^^j-an^^gg fQ^- gxample that 
U^'^ = whenever /i = l,2orz/ = l,2. Thus, effectively the Dirac indices run only through 
(1 + 1) dimensions (/i = 0,3). One can show that (cf. Appendix E): 



Tr [VoM] 



^Tr [M] 



(5.32) 



for any given matrix Ai that is a linear combination of 1,7^ and '^^'^^ , with ^ 1,2. 
The trace T4 entering the gluon selfenergy computation is therefore equivalent to a trace in 
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dimension c? = 2, in which the Dirac indices assume only two values (but the matrices are 
still 4x4): 



^Tr,=2 



1 



1 



(5.33) 



Using this result, the nonzero components of the one- loop gluon selfenergy in a magnetic 
background within the LLL approximation may be rewritten in terms of the vacuum gluon 
selfenergy in dimension (1 + 1): 

^1 ^2 1 ^ 



iU^''{k,k') "'"^''^ (27r)V(A;' + A;)^exp 



2eB 



with 



U^^^iko, ks 



^9 



^Trj=2 



(2vr)^ 



1 



r 



Vl-Il- ruf 



-1" 



[Vl 



7l - m/J 



(5.34) 



, (5.35) 



or, solving the Dirac trace in the standard way [145], 



{ko,k^ 



dpodps 
(2vr)2 



4 K - PL ■ (PL + ^lW-" + p^^ik + pY 



+p''{k + py 



(5.36) 



Notice that the theory in dimension 2 is in general much better behaved in the UV as 
compared to c? = 4, requiring no renormalization procedure. The integral above, however, is 
not convergent, but, as we shall see below, has no poles in 1/e in dimensional regularization, 
being thus of the form: J d'^pp"'. The absence of poles in dimensional regularization guar- 
antees that no nontrivial subtraction is needed, so that the final result does not depend on 
a renormalization scale at all, as expected. 

Let us now obtain an explicit form for Il^'^{k, k'), Eq. (5.34). The result for the one-loop 
vacuum selfenergy in d dimensions is of course well-known and can be obtained straightfor- 
wardly using dimensional regularization and Feynman parameters [145]: 

N?-l fdpodps 



4{ [mj - PL ■ (PL + + 'P^{k + pf + 



-^p^ik + p) 



e->0 




k^'k'' 

^L 




(5.37) 
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with 



(5.38) 



'0 yiX-y)-x 1 

Therefore, the final result for the gluon selfenergy corrected by a loop of fermions^ dressed 
by the magnetic field in the LLL approximation is: 



I{^^-{k,k') ^'"^^'^ (2vr)V(A;' + fc)— exp 

271 



(5.39) 



k i + \g'Ni-l 



2eB 



71 



ki'k' 



m 



f, 
(5.40) 



reproducing the analogous photon result [142] for g^^N^ — l)/2 ^ e^. 

As discussed in the QED context, the function f{x) presents a divergence at the threshold 
for decay of the photon into an electron-positron pair and an imaginary part signaling this 
instability for the high momentum modes. Since the position of the threshold is of course 
dependent on the fermion mass mj, the chiral limit is subtle and is not interchangeable with 
the zero- momentum limit. This feature may be problematic for the thermal version of this 
computation, since the integrals over pq and p are naturally separated by the discretization of 
the zeroth component into Matsubara frequencies. The computation of the thermal one-loop 
gluon selfenergy in the presence of an intense magnetic background is work in progress. 



5.4 QCD pressure 

After the experience acquired in the computation of the gluon selfenergy in our framework of 
perturbative QCD in a nonperturbative magnetic background, let us address the thermody- 
namics. It is interesting to recall that these results may be directly compared to high-quality 
lattice QCD simulations with a magnetic field that are currently underway, with some results 
already reported in Ref. [138]. 

Once again, we shall use perturbative QCD with the dressed fermion propagator from 
Section 5.2. The standard diagrammatic expansion for the thermodynamic potential of QCD 
reads: 

"'in the case of quarks, the electric charge is flavor-dependent: e i— > gj. 



CHAPTER 5. THERMODYNAMICS OF QCD IN A MAGNETIC BACKGROUND 93 




+ [diagrams with counterterms] + 0(3 loops), (5-41) 

where full lines are fermions, dressed by the magnetic field, curly lines are gluons and dashed 
lines represent ghosts that will cancel the spurious degrees of freedom in the gluonic pressure. 

To investigate the role of color interactions, our computation will go up to two-loop order. 
In this approximation, the thermodynamic potential may be separated into a fermionic part 
{Qq(^j^{B)), that carries all the magnetic- field dependence, and a pure glue contribution 
{^qcd) ^he following way: 

^QCD = ^qcd{B) + ^QCD + [diagrams with counterterms] + 0(3 loops) , (5.42) 

where 




(5.43) 
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The gluonic part is equivalent to the usual hot perturbative QCD result and is, therefore, 
well-known (cf. Appendix F for references and details on the collection of renormalized 
results) : 

n^CD = -2(iVc-l)^ + (iV.^-l)iV./T^^. (5.45) 

In the presence of a constant and uniform magnetic field, however, one needs to compute 
the nontrivial contribution Qqqj^{B), which contains two terms: 0{g^) and 0{g'^). 



5.4.1 Pressure of free magnetically-dressed quarks 

The first diagram in Eq. (5.43) corresponds to the free Fermi gas pressure in the presence 
of the magnetic background. In the free theory, the action becomes quadratic, so that the 
partition function can in principle be solved exactly, in terms of the determinant of the 
inverse propagator which in this case is the exact fermion propagator including the magnetic 
dressing. The computation of this determinant, however, turns out to be prohibitively 
complicated due to the absence of translational invariance of the propagator. 

For the LLL approximation, the problem is worsened by the presence of vanishing eigen- 
values associated with the spurious degrees of freedom with spin anti-aligned with the mag- 
netic field, which are physically absent from the LLL. It is crucial to realize, however, that 
the LLL approximation for the free gas pressure is not equivalent to the leading order of 
a large magnetic field expansion at all. As it shall become clear in what follows, for the 
zero-temperature, finite-S term of the pressure, the LLL is the energy level which less con- 
tributes in the limit of large B\ the result being dominated by high values of n. Nevertheless, 
the equivalence between the LLL approximation and the large B limit remains valid for the 
temperature dependent part of the free pressure. 

This free contribution to the fermionic pressure has been considered in different contexts 
(usually, in effective field theories [131, 132, 146, 147, 148]) and computed from the direct 
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knowledge of the energy levels of the system, i.e. the Landau levels^ 

E'^{n,p3) = pl + mj + 2qfBn, (5.46) 

and their degeneracies q/B/ (27r) for n = and qfB/{7i) for n = 1, 2, ■ ■ ■ . The exact result, 
including all Landau levels, is: 

= -NcVY^^-^{l - <5„o/2) J ^|/3i5(n,p3) + In (1 + e-/^[^("'^3)-M,]) + 

+ ln(l + e"'^[^("'^'^)+'^^]) |, (5.47) 

where the first term is a clearly divergent zero-point energy term and the other two are 
particle and anti-particle finite-temperature contributions. Since E{n,p-i) grows with B, the 
largest the n labeling the Landau level considered the larger the zero-point energy term 
becomes, being minimal for the LLL. Thus, in the limit of large B, the LLL approximation 
is inadequate. The decaying exponential dependence of the finite-temperature term on the 
energies E{n,p3) guarantees, however, that the LLL dominates indeed the result for intense 
magnetic fields. 

To obtain a good approximation for the large B limit of the free pressure, we choose to 
treat the full exact result and take the leading order of a m'j/qfB expansion in the final 
expression. Let us then discuss the treatment of the divergent zero-point term. Despite be- 
ing a zero-temperature contribution, the first term in Eq. (5.47) cannot be fully subtracted 
because it carries the modification to the pressure brought about by the magnetic dress- 
ing of the quarks. One needs therefore to address its renormalization. Using dimensional 
regularization and the zeta-function representation^ for the sums over Landau levels and 



^It is equivalent to describe the Landau levels in terms of a orbital number I = 0, 1, 2, • • • and the spin 

quantum number s ~ ±1 as E'^{1, s,p3) = P3 + ni'j- + qfB{2l + s + 1). In this representation the degeneracy 

is the same for every state 77;. s = qfB/{2TT). 

^Notice that this is also a type of regularization. 
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subtracting the vacuum term in (3 + 1) dimensions, one gets (def.: Xj = m'j/{2qfB)): 

' ^ IT 



6no/2) I ^E{n,ps] 



27r2 
1 



f 

2 



12 



(2/e + ln(AV2g/5) + 1) , (5.48) 



where a pole ~ (g/i?)2[2/e] still remains. This infinite contribution that survived the vacuum 
subtraction can be interpreted as a pure magnetic pressure coming from the artificial scenario 
adopted, with a constant and uniform B field covering the whole universe. In this vein, one 
may neglect all terms (~ QfB"^) and independent of masses and other couplings (as done, 
e.g. in Refs. [131, 147, 148]), concentrating on the modification of the pressure of the QCD 
matter under investigation. 

The final exact result for the free pressure of magnetically dressed quarks is therefore 
(pressure = —(1//?^) [diagram]): 



pF 

free 



E 



27r2 



c {-i,xf) - c (-1, 0) + 2 (^/ - 4) ^^^f + ^ 

+ Tj2^a - Sno/2) / ^1 In (1 + e-/^[^("'^'3)-/^/l) + In (l + e-/5I^("'^'3)+M/]) 



(5.49) 



In the limit of large magnetic field (i.e. Xf = m'l/{2qfB) — )■ 0), one obtains: 



free largo B 



27r2 



Xf 



+ 



qfB f dp: 



2n 



2n 



In (1 



^/3[E(0,P3)-M/: 



) +ln(l 



-/3[£;(0,p3)+M/r 



(5.50) 



Adding the 0{g^) piece of the gluonic contribution in Eq. (5.45), the free pressure of the 
quark-gluon system in the presence of an intense magnetic field reads: 



pF+G 
free 



2(iV2 _ 1^ 



90 



f 



(5.51) 



Figure 5.1 illustrates the behavior of P^j.^ {T^B) for two-fiavor QCD with m„ = md 
5 MeV. 
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Figure 5.1: Free pressure of a gas of gluons and magnetically-dressed up and down quarks 
as a function of temperature. Different curves represent different magnetic background 
intensities: eB = {O (black, solid, bottom), 20 m^, 40 ml., 60 (magenta, dash-dotted, 
top)}, where = 138 MeV is the vacuum pion mass. 

5.4.2 Turning on quark-gluon interaction: the exchange diagram 

in a magnetic field 

The second diagram in Eq. (5.43) is referred to as the "exchange energy" and corresponds 
here to the first nontrivial contribution that perturbs the decoupled system of gluons and free 
dressed quarks by turning on the SU{Nc) gauge interaction. Its computation and analysis 
is our main goal in this on-going project. 

In terms of the propagators in coordinate space (notice that logZexch = (—1/2) [diagram]), 
the exchange diagram reads: 

er rjdr^M^ r Adjy- -iK-{y-x) 

i;f (5.52) 

where Aq are Gell-Mann matrices, with AqAq = the trace Tr acts over Dirac indices 

and the 4-momentum is given in terms of the Matsubara frequencies {uf = 21tiT) and of 
the 3-momentum as: K = [k^ = iuf , k) . 
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Notice that if we assume translational invariance and the free Dirac propagator for the 
fermions (with P = {jp = iu^ + fif , pj and = {2n + l)vrT), 



d'^P 
(27r)° 



^dp Q-iP-ix-y) 



m 



(5.53) 



/ 



this expression reduces to the usual one (cf. e.g. Eqs. (5.36) and (5.49) in Ref. [149]): 

1 1 




(3V 



g^uA {2nff3 6^'\K-P, + P2)TT 

JPl,P2,K 



(5.54) 



where we have defined = Nc{N^ — l)/2. 

In the presence of a uniform and constant magnetic background (B = B z), however, the 
fermion propagator becomes dependent on x and y in a nontrivial way due to the Schwinger 
phase, as discussed in Section 5.2. 

Introducing the LLL form of the fermionic propagators, Eq. (5.22), in the exchange 
contribution, Eq. (5.52), we obtain: 

" d'^xd'^y f d'^K q-'^^'^v-^) qj-B f dpodp2dp3 f dqQdq2dq3 




f3V V, 



(3V J {2ttY 



(27r)"'-i J {27rY~^ 



exp\ - i{po - qo){xo - yo) + i(p2 - q2)ix2 - 1/2) + i{P3 ~ (lz){x^ - y^) 



exp 
Tr 



q^B 



_P2_ 

q^B 



Xi 



12 

qjB 



yi 



P2 

qfB 



yi 



J2_ 

qjB 



7m 



T- 



(5 



As argued in Appendix E, from the properties of Dirac gamma matrices and the defini- 
tion of the projector Vq, one concludes that the trace appearing in Eq.(5.55) is essentially 
occurring in a (l+l)-dimensional subspace (/i = 0, 3). Using d = 2, Eqs. (E.l) and (E.3) and 
the standard relations for Dirac traces (Tr[(7^)"] = 0, with n odd; g'^^ = d; Tr[7^7,^] = 4 g^^), 
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we get (defining the 2- dimensional slash as: = ■ 7^): 

Tr [7^Po(/' + m;)7^Po(^ + rn^)] = Tr [Po7^(/' + my)7^(^ + 

= 2{2-d) pL-(lL + 2dm}, (5.56) 
so that the exchange contribution for the LLL becomes 

LLL 




d'^K f d^Pd^Q 2 (2 - d) pL-(iL + 2dmj q^B f dp2dq- 



{27tY J (27r)2'i" K^IpI - m'j][ql - mj] ir J (27r) 

^/ f ^ y f^-iK-iy-x)-i{po-qo){xo-yo)+i{p2~q2){x2~y2)+i{Pa-q3){x3-y3) 

2 / \ 2 



] } ■ (5.57) 



, 2 / \ 2 

P2 \ , / q2 



The integrals over x and ?/, in the second line above, give (def.: PV = PL1L2L3): 

^'^^^'^y ^-iK-{y-x)-i{po-go)(xo-yo)+iip2-q2){x2-y2)+i{P3-q3)(x3-y3) _ f ^i^l^„ifci(s/i-a;i) 



PV ^ J /3V " 

. j dx^dy^dx^dy^dx^dy..-^^^^^^^^^^^ 

= ^(27r)5(po - go - h){2'n)6{p2 - q2 - h){2n)S{ps - q, - h) j ^^^^^^.iv^-^r) 
= ^(2nf6{pL -qi- kL){2n)6{p2 - q2 - h) J ^^^e^^'^^^-^^) . (5.5^ 
We are then left with: 
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LLL 



d'^K r SPd^Q2{2-d) pL-qL + 2drr 



-(27r)'^5(pi-qz.-ki) 



exp 



TT 



(27r) 



-(27r)(5(p2 - 92 - /i;2) 



dkidk2 
(27r)2 



_P2_ 

9/5 



92 



P2 



yi 



92 

9/5 



]} 



d^Kd^Pd^Q 2{2-d) pL-qL + 2d 
dxidyi if,^^y^_^^)qfB f dp2dq2 



m. 



-i2n)%pL-ciL-kL) 



(5.59) 



TT 



(2vr) 



exp 



_P2_ 

9/5 



(27r)(5(p2 - 92 - h) 

2 



92 

9/5 



P2 

9/5 



92 

9/5 



]} 



(5.60) 



Defining the function tliat comes from the integrations related to the transverse directions 
as 



I{ki,k2,qfB) 



^e*(--)^ / #%(2vr)5(p2 - 92 - A.^) 



X exp 



TT 7 (27r)2 

\ 2 

P2 



+ X 



92 

9/5 



2/1 



_P2_ 

9/5 



92 

9/5 



]} 



(5.61) 



and reahzing that the hne labeled as Eq.(5.59) is nothing but the 2-dimensional result for 
the exchange contribution including a "massive" gluon {K'^ = k| — ml; ml = kf + k2), we 
have: 




LLL 



dkidk2 
(27r)2 



I{ki,k2,qfB) g {kl + klm}) . 



(5.62) 



CHAPTER 5. THERMODYNAMICS OF QCD IN A MAGNETIC BACKGROUND 101 



with Q (mljrnj^ being the 2- dimensional exchange diagram: 




d=2; ml=kl+k? 



2 



SkSpSQ 2{2-d) pL-c[L + 2dm 



/J 



(5.63) 



Finally, taking the result in Eq. (E.17) for the integral over the transverse directions 
J(/ci, /c2, g/-B) (Eq. (5.61)) obtained in Appendix E, 




I{h,k2,qfB) = . (5.64) 

ZTT 

into Eq. (5.62) for the exchange contribution to the QCD pressure, one obtains: 

e-47#e(,;^4»?) (5.65) 

This expression realizes concretely the intuitive expectation that the nontrivial dynamics 
in an extremely intense magnetic field should be restricted to (1 + 1) dimensions. Since 
the gluons do not couple directly to the magnetic field, their dispersion relation maintains 
its (3 + 1) -dimensional character (w^ = k\ + k\ + A;|), what effectively results in a "mas- 
sive" gluon (m^ = kl + k^) in the reduced (1 + 1) -dimensional diagram. In the end the 
exchange contribution to the QCD pressure in the lowest-Landau level approximation for 
the fermion propagation is essentially an average over the effective gluon transverse mass 
m\ = k\ + k^ of the exchange diagram in (1 + l)-dimensions with the Gaussian weight 
(g^5/27r)exp[-m2/2g^5]. 
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Let us now compute the 2-dimensional massive exchange diagram. Departing from Eq. 
(5.63), 

g{m,^m,) = pVg v. j 5{p, - q. - k,) _ _ _ 



y --^^ {2nr ^5(pL-qL-k^ 



v-^ 2{2-d) pL-ciL + 2dm} 

4^ ir tr ^ """^^'[ki-m^][pi-m|][qi-m}] 



(5.67) 



Apart from the tensorial structure, this diagram corresponds to the 2-dimensional version of 
the exchange diagram of in-medium Yukawa theory with both fermions and bosons massive 
that was computed originally in Ref. [47] (cf. also [48]) and is also used in Chapter 6. In 
dimension (1 -|- 1), however, UV divergences are much less severe, being absent in this case, 
and one expects renormalization to be trivial. Our concern here will be directed towards the 
IR domain, which should be subtle as proved to be the case for the selfenergy above. 

In principle the diagram computation in dimension 2 is fully convergent so that one may 
set d = 2: 

g{mlm)) = /3Vg'u,l^^0^i2n)6{p,-q,-k,) 



[ki - ml] [pI - mj] [q2 - mj 



(5.68) 



It is interesting to note from this result that in dimension 2 the exchange contribution to 
the QCD pressure is exactly zero in the chiral limit — t- 0. The magnetic QCD case is, 
however, more intricate and the extra ki, ^2 integrals in Eq. (5.66) could effectively produce 
a different result for the chiral limit or even bring IR divergences that should be treated 
carefully. It is not clear, for example, whether one should evaluate the Matsubara sums 
before the momentum integrations or not (as discussed in Ref. [142]). 

Let us proceed in the usual way, evaluating the Matsubara sums first, via the standard 
contour integration procedure [149]. 
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Using Kapusta's representation for the Kronecker Delta [150]^, 

JO Jo 

we can rewrite Eq. (5.67) as: 



(5.69) 



where 



(5.70) 



S{p,,q,,k,) = 77-2 ■^o(PO'PO'^o) (5_71^ 

^-^ [k^ — m^J[p^ — m^J[q^ — 'm^-l 



/J 



so(po,go,A;o) = [2{2-d} pl ■ cIl + 2d mj] . (5.72) 

Po ^ Qo ^ 1^0 

Notice that Sq{pq, Qq, k^) is an analytic function, even in the limit Pq — (go + ^o)- 

hm So{po,qo,ko) oc e^h+.o-.,] ^ - + ^(Po - - ko)] ^ ^,[.,+,o-.,] (5,73) 
po->{(jo+fco) J'o ~ 90 ~ kQ 

Therefore, solving the Matsubara sums using the standard results 



^Epf^ = ^{^b>H/-/^/)}i:::\ (5.75) 

with g(x) analytic, = -{uff - k^, = {lu^^ + ^ijf - p^, = {luj^^ + - eg 
u"^ = ml + k^, E^ = pI + mj, nsiu}) = [e^" — 1]^-^ and np^E) = [e/^^ + we arrive at: 



^This substitution actually allows for the computation of the sums in any desired order, solving a subtle 
problem related to the absence of unicity of the analytical continuation of a function defined in a discrete 
domain [151]. 
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>S(p^,q^,k^ 



I ni 712 



I 

1 



— K"^ ^ — m?f ^ — 77^2 



ScuEpEq 



{pQ,qo,ko)={Ep,Eq,uj) 



{pQ,go,ko)={-Ep,-Eq,-Lu) 



(5.76) 



\\ |X=a "I \y=b 



where we adopt the notation {/(x, 2;)}||^;^j5=[!_t^fe,-c) ^ i_ lv^ v-, ix=-aj ij,=_fe 

In order to further simphfy Eq. (5.76), we can rewrite the function So{po, Qq, /cq) through 
the identity 

1 111 



^/3[fco+'Jo-M/] _ p/3[P0-A'/] 



npiqo - IJ'f)nB{ko) ^^(fco) npiqo - fJ'f) npipo-jJ'/) 



(5.77) 



Finally, using these results in Eq. (5.70), after a long but straightforward algebra, we 
obtain: 



(27r)3('^-i) 



(27r)'^-i5(k, - p, + q. 



StuEpEq 



J_ E. - J, E_ 



J+{E_ +uj)- J.{E+ - Lo)\ Nf{1) -2J_E+ UBiu 
-J.{E+-uj)\, 



(5.78) 



where 



J± 



cl 9 
-2 ^ f 



1 (±EpEq-p,-q, 



2m) 



E^-u^ 



E%-uj^ 



(5.79) 



and 
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u 

Nf{1) 
Si 



m, 



(5.80) 
(5.81) 
(5.82) 
(5.83) 
(5.84) 



In d = 2 explicitly, we have: 

m - 1 



e;(mtm^) = 13V g'N, 
1 



X 



2 

u Si 



(27r)3 



(27r)(5(/c3 -P3 + 93) X 



wEpEq \El-uj'^ El- u;2 



^4 



2{E^ + Lo) 



_{E^ -L0){E++L0) 



Nf(1)-2 



El-u^ El-ul 
E+ 



El-u^ 



(5.85) 



E+ + UJ 



Finally, our preliminary result for the exchange contribution to the pressure of perturba- 
tive QCD in a magnetic background is®: 




LLL 



m. 



2 / f 
qfB\ f dkidk2 -^tf f dpsdqsdh 



2tt 
1 



e 



(27r) 



E, 



2iE^ + u) 



{E_ - io){E+ + io) 



Nf{1)-2 



El - W2 ^2 _ ^2 

^^(cu) 



El-u^ 



E.+u 



(5.86) 



Interestingly, the result seems to be proportional to the quark mass, although the re- 
maining momentum integrations may still bring l/m^ factors. If not, this means that in the 
chiral limit the gauge interaction does not correct the pressure up to 0{g^). 



*For the full result in QCD, summation over the flavor index / is still required. 
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To obtain a final result, the momentum integrations need to be carried out explicitly. 
However, the possibility of divergences, especially in the IR domain, is not excluded and a 
regularization procedure may be required, as usually happens in the description of a magnetic 
background. 

5.5 Remarks and perspectives 

In this chapter, wc have considered an external magnetic field as an extra axis in the phase 
diagram of hot QCD matter. As usual, the new axis opens up novel channels of comparison 
between different theoretical approaches, but in this case also between theory and experi- 
ment. As discussed in the first section, several recent developments have shown that hot 
QCD matter under intense magnetic fields has an exciting phenomenology, closely connected 
with on-going heavy-ion collision experiments at RHIC and LHC, and is amenable to first- 
principle lattice QCD simulations. 

Our approach here [127] goes in the direction of furnishing predictions for the thermo- 
dynamics of QCD in a magnetic background within the fundamental gauge theory. In this 
vein, we adopt an approach that combines perturbative QCD with a nonperturbative treat- 
ment of the coupling of quarks to the intense magnetic field. After the building blocks of 
the framework were constructed, the application to the computation of the one-loop vacuum 
selfenergy for the gluon and the preliminary results for the thermodynamics have revealed 
interesting features and subtleties intrinsic to media under extreme magnetic backgrounds. 

The effective dimensional reduction to (1 + 1) dimensions was verified in different con- 
texts, in agreement with the intuitive physical picture of Landau orbits in the limit of intense 
magnetic fields. Although the dimensional reduction tends to soften possible UV divergences, 
we have seen in several examples that regularizations are still required. The computation of 
the free pressure of magnetically-dressed quarks revealed an (infinite) pure magnetic contri- 
bution due to the artificial picture of an omnipresent constant and uniform field. Moreover, 
we have shown that the connection between the LLL approximation and the large B limit 
is not universal. A lot, however, is still to be understood and the subtleties encountered so 
far point out that the remaining momentum integrations in the exchange contribution, Eq. 
(5.86), should be carefully analyzed. Near future perspectives of this on-going work [127] 
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include as well the computation of the one-loop quark and gluon thermal selfenergies within 
this framework and comparison to lattice data. 



T 



tiean-field 



Chapter 6 




m 



Masses and condensates in 
nonperturbative thermodynamics of 
Strongly interacting matter 



Towards the goal of building a robust description of the phase diagram of Strong Inter- 
actions via effective models, one should pursue an approach that can accommodate finite 
temperature and chemical potentials, different masses, condensates and nonperturbative in- 
teractions. Since order parameters are often related to condensates, which in turn commonly 
appear in effective masses, it is desirable to have a framework that combines a nonperturba- 
tive method with the capability of investigating the role played by different massive fields. 

This chapter presents a nonperturbative computation of the thermodynamic potential of 
a general Yukawa theory within Optimized Perturbation Theory to two loops, allowing also 
for the presence of a scalar condensate [152]. We shall keep the discussion in general grounds, 
concentrating on the analysis of the nonperturbative method in question, as compared to 
naive perturbation theory and mean-field approximation. The results are applicable to a 
wide range of interacting fermionic environments; not only Strong Interactions, but also 
the Higgs Sector of the Standard Model and Condensed Matter systems described by the 
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Hubbard Model, as argued in the introduction. 

6.1 Introduction 

Matter under extreme conditions often requires the use of effective field theories in the 
description of its thermodynamical properties, independently of the energy scale under con- 
sideration. For instance, in the phenomenon of spontaneous symmetry breaking and tem- 
perature dependence of antiferromagnetic order in high Tc superconductors described by a 
field-theoretic version of the Hubbard model [153] at relatively low energies, one has to incor- 
porate the interaction of fermions with the crystal lattice via phonons, in a coarse graining 
of quantum electrodynamics. Moving up in energy scale, the phase structure of strong in- 
teractions and the quark-gluon plasma [11], investigated in high-energy heavy- ion collisions 
experiments [25] and in the observation of compact stars [32] , often demands simplifications 
of quantum chromodynamics, producing a variety of low-energy effective models [116]. 

Besides the major role it plays in the mechanisms of spontaneous symmetry breaking 
and mass generation in the Standard Model of particle physics, the Yukawa interaction 
stands out as one of the main ingredients in the construction of simplifying effective field 
theories to study the thermodynamics of systems under extreme conditions, especially if one 
imposes renormalizability. Although an effective theory does not require renormalizability 
to be consistent, a physically- motivated cutoff being usually more than satisfactory in this 
case, this feature can prove to be useful in the study of scale dependence and running via 
renormalization group methods [47]. 

In this chapter we investigate the full nonperturbative thermodynamics of the Yukawa 
theory at finite temperature and density by computing its pressure using the optimized 
perturbation theory (OPT) framework. In our evaluations we consider contributions up to 
two loops, which include direct (Hartree like) as well as exchange (Fock like) terms, so that 
the Yukawa thermodynamics can be investigated in the presence of condensates and also in 
their absence. In the first case, besides OPT and ordinary perturbation theory (PT) we shall 
also perform a mean field (MFT) evaluation. This is an important step in establishing the 
OPT reliability since the reader will see that when exchange contributions are neglected OPT 
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exactly reproduces MFT results which can be considered "exact" in this (large A^) limit^. We 
present results which are valid for arbitrary fermion and scalar masses, temperature, chemical 
potential, and coupling. The region of large values of the coupling is particularly interesting 
and useful, since several effective field theory models in particle and nuclear physics exhibit 
Yukawa coupling constants that are much larger than one, as is the case, for instance, in the 
linear sigma model [79], described in Appendix A and frequently used in the description of 
the chiral phase transition (e.g. in Refs. [80, 81, 82, 155, 156, 157, 158] and in Chapters 3 
and 4 of this thesis) and in pion-nucleon models extracted from chiral Lagrangians [159, 160]. 

When condensates are present we show that, as expected, PT has a poor performance 
since it cannot resum direct (tadpole) contributions associated with symmetry breaking. 
On the other hand, OPT improves over MFT by also incorporating exchange terms in a 
nonperturbative fashion. 

In the absence of condensates and in the regime of very small coupling our numerical 
findings can be verified in the limit of vanishing temperature by comparison to exact analytic 
two-loop perturbative results previously obtained for the equation of state of cold and dense 
Yukawa theory within the MS scheme [47] (cf. also Refs. [48, 161]). In Ref. [47], the two-loop 
momentum integrals were computed analytically for arbitrary fermion and scalar masses, the 
final result being expressed in terms of well-known special functions, which provides us with 
a solid and clear reference in this limit. 

Furthermore, this comparison also provides another way of testing the idea that pertur- 
bation theory at high density and zero temperature in the symmetric phase is much better 
behaved than its converse [47, 106, 162, 163, 164], which has well-known severe infrared 
problems [149]. Our analysis shows that, surprisingly, second-order PT agrees quite well 
with the OPT nonperturbative results up to values of the coupling of order one in the case 
of cold and dense as well as hot and dense Yukawa theory. 

The framework of OPT [165], also known as the linear delta expansion, is an example 
of a variational method that implements the resummation of certain classes of Feynman 
diagrams, incorporating nonperturbative effects in the computation of the thermodynamic 
potential (for related methods, see Refs. [166]). It has been successfully applied to the 

^Thc reliability of the OPT framework applied to symmetric and broken phases was studied previously 
in the eontext of (/)^ scalar theories (see for instance Ref. [154]). 



CHAPTER 6. MASSES AND CONDENSATES IN NONPERTURBATIVE. 



Ill 



study of many different physical situations, such as mapping the pliase diagram of the 2+1 
dimensional Gross-Neveu model [167], where a previously undetermined "liquid-gas" phase 
has been located, and determining the critical temperature for Bose-Einstein condensation 
in dilute interacting atomic gases [168]. Some early applications at finite temperature can 
be found in Refs. [169]. 

The chapter is organized as follows. In Section 6.2 we present the in-medium Yukawa 
theory and set up the notation to be used in what follows. Section 6.3 contains a sketch 
of the derivation of the perturbative thermodynamic potential at finite temperature and 
chemical potential. In Section 6.4 we apply the OPT machinery to the evaluation of the 
nonperturbative thermodynamic potential in the Yukawa theory. The thermodynamics of 
the Yukawa theory, in the presence of condensates, is considered in Section 6.5 where the 
results from PT, OPT and MFT are contrasted. In Section 6.6 we consider the symmetric 
case, where condensates are absent. Comparing the results from OPT and PT, we show 
that the latter performs well also for extreme temperature values up to two-loop order. 
Section 6.7 contains our conclusions. Technical details involved in the calculation of the 
vacuum contributions and direct terms in the two-loop thermodynamic potential are left for 
appendices B and C^. 

6.2 In- medium Yukawa theory 

In what follows, we consider a gas of Np flavors of massive spin-1/2 fermions whose inter- 
action is mediated by a massive real scalar field, 0, with an interaction term of the Yukawa 
type, so that the Lagrangian has the following general form: 



^The technicalities concerning Matsubara sums, renormalization, and the analytic evaluation of in- 
medium momentum integrals in the calculation of the medium contributions up to two loops were extensively 
addressed in the appendix of Ref. [47] (cf. also [48]). 



■int 1 
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where 

Np 

= ^ip^{i(^ - m)ipc,, (6.2) 
^ = l{d,<PW<P)-^ml<f>'-Xs<l>'-X<l>\ (6.3) 

Np 

l^int = ^9 • (6-4) 

a=l 

Here, m and are the fermion and boson masses, respectively, assuming all the fermions 
have the same mass, for simplicity. The Yukawa coupling is represented by g; A3 and A are 
bosonic self-couplings allowed by renormalizability. Here we choose A3 = A = 0, disregarding 
bosonic self-interactions. 

We work in the imaginary-time Matsubara formalism of finite-temperature field theory, 
where the time dimension is compactified and associated with the inverse temperature (3 = 
1/T [149]. In this approach, one has to impose periodicity (anti-periodicity) for the bosonic 
(fermionic) fields in the imaginary time r, in order to satisfy the spin-statistics theorem. 
Therefore, only specific discrete Fourier modes are allowed, and integrals over the zeroth 
component of four-momentum are replaced by discrete sums over the Matsubara frequencies, 
denoted by = 2miT for bosons and = (2n+l)-7rT for fermions, with n integer. Nonzero 
density effects are included by incorporating the constraint of conservation of the fermion 
number via a shift in the zeroth component of the fermionic four-momentum p° = iu^ 1— ?• 
p° = iu^ + /i being the chemical potential. 

Figure 6.1: Perturbative thermodynamic potential for the Yukawa theory to two loops. 
Solid lines represent fermions and dashed lines stand for bosons. Here we omit the diagrams 
containing counterterms. The third (tadpole) contribution can be neglected when we consider 
the symmetric case. 

From the partition function written in terms of the euclidean action for the lagrangian 
(6.1), Zy(T, /i) = Tr exp(— Sy), one derives the perturbative series for the thermodjTiamic 
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potential Qy = -{l/PV) In Zy [149]: 

~' (-1) 



(6.5) 



where V is the volume of the system, Zq is the partition function of the free theory and 
Sint represents the euclidean interaction action. Notice that Wick's theorem implies that 
only even powers in the above expansion survive, yielding a power series in ay = g'^/^n. 
However, at finite temperature the perturbative expansion also contains odd powers of g 
coming from resummed contributions of the zeroth Matsubara mode for bosons, such as 
in the case of the plasmon contribution [149]. In the zero-temperature limit this is not 
the case and, even including hard-dense-loop corrections, only g"^^ terms are modified [163]. 
Since we restrict our analysis to two loops (the resummed OPT calculation also departs 
from a two- loop perturbative setting), we are not concerned with this issue. Omitting the 
diagrams representing counterterms, the perturbative thermodynamic potential is shown, 
diagrammatically, in Fig. 6.1. The first two diagrams correspond to the free gas, while 
the last two represent interaction terms: the third diagram is a contribution of the direct 
type and the fourth is of the exchange type. Notice that the third diagram, which contains 
tadpoles, belongs to the one-particle-reducible class, that does not contribute in the absence 
of condensates. It will therefore be neglected in Section 6.6, where only the symmetric case 
is considered. 



6.3 Perturbative Thermodynamic Potential 

Considering the thermodynamic potential only up to the first non trivial (two-loop) contri- 
bution, which corresponds to direct plus exchange terms, one must evaluate the following 
0((7^) quantity: 

= ^ ^In [K' - mj] + z ^Trln [f - m] - ^ ^Tr 

where the trace is to be performed over the Dirac structure, and the 4-momenta are given in 
terms of the Matsubara frequencies and the 3-momenta for fermions, P = (p*^ = iu^ + /i , p) , 
and bosons, K = = iojf , k) . We choose the metric tensor signature g^'^ = diag(-|-, — , — , — ), 



Sdir(O) 



Scxc(-P) 



(Z'-m) 



(6.6) 
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and the following short-hand notation for sum-integrals: 

(27r)3 ■ 



The direct contribution is given by 

Sdir(O) = 



Tr, 



(6.7) 



(6.8) 



and, the exchange self-energy is given by the following sum-integral over fermionic momenta: 

1 



Scxc(^') = ig 



Q[(?-m][(P-g)2- 



(6.9) 



Computing the traces at /i 7^ and T 7^ (see Ref. [47] for details concerning the 



exchange term), one can write the thermodynamic potential as VLy 
where the vacuum contribution prior to renormalization has the form 



(2^ 



(2 



3^p 



2 r 

F 



2mi 



2m 



P 1 



(27r)3£;p 



d^Pi d^p2 J^{E+ - LO12) 



4 J (27r)3(27r)3 oo^E^.E, 



P2 



where Wk = (k^ + "t,^)^/^, Ep 
and we have defined 

J± = -2 



(p2 + ^2)1/2^ ^ ^ ^^^^ = (|p^ 



P2I 



m 



(6.10) 

2N1/2 



m 



Pi • P2 ± Ep^Ep^ 



Am? 



??2 , ,2 



El 



(6.11) 



Notice that Eq. (6.10) contains zero-point energy divergent terms which are T- and fi- 
independent only within the PT approach in which case they can be conveniently absorbed 
by the usual vacuum subtraction which normalizes the pressure so that it vanishes at /i = 0. 
In practice this means that, as done in Ref. [47], one does not have to care to their explicit 
evaluation and renormalization. However, within the other approaches (OPT and MFT) 
considered here the PT bare mass, m, is replaced by effective (T- and fi- dependent) masses 
within the same diagrams. Therefore, we must renormalize those contributions appropriately. 
Within the MS subtraction scheme, the fully renormalized vacuum term can be written as 
(see Appendices B and C for details) 



-5^™(m<^) + 4iVi.5^^^^(m) + a 



REN/ 



^dir,REN, 



m, m^) + n 



exc,REN, 



m,m^) 



(6.12) 
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where we have defined the following functions of the masses: 



647r2 



- + log ( 



V M2 



n 



dir,REN, 



cxc,REN , 



2mi 



I (27r)^ 



1 + ln 



A2 



(6.13) 
(6.14) 



2 4 

g m 
^ 2"64^ 



^^1 



4m2 



7 + log 



V2 



4m2 



+ 



1 

4 



7 + log ^ 



^^3 



4m2 



+ 6 a (m^) 



1 - 6- 



(6.15) 



Here 7 is the Euler constant, A is the renormalization scale in the MS scheme, and the 
functions Vi{z) and a (m^) are defined in Appendices B and C. 

Adding the direct term to the free gas plus exchange medium-dependent terms evaluated 
in Ref. [47] one obtains 



27r2 



z dz\og[l — e~ 



_ T^f!£ / z'dz {log[l + e-(^-^/^)] + log[l + e-(^^+^/^)] } 
^ Jo 



2 2 ^1^^ 
- g^T^ ^ 



(47r4)m2 



, /A2 
1 + In — 



z dz- 



z^+m^/TY'^ 



X 



g'T 



1 



2^4 ^'f^' 



(27r4)m2 



2; (i^- 



2Tn2 2 

— g 1 m 



'4> vJQ 



+^2/2^2] 1/2 



+ 



(2vr) 

- 9^T^77^ ("2 + 3^3) 



/■oo 


r^/(i)i 


/o 





2Tn4 



+ g'T 



(2vr)4 
2(2^ 





nb(w^) 


/o 





z^dzy'^dy l du 



"zy 



EyE> 



z dzx dx I du 



1 



1C^E+ - iC+E^ 



nb{uJx)Nf{l) 



(6.16) 
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+ m'/T^ E^ = E,± E,,, E^ 



where, in order to perform numerical investigations, we have defined the following di- 

mensionless quantities: oo^ = z'^ + m^^/T^, _ ^2 , ^2 /t^2 _ z? -u z? tp2 
+ 2^ + 2xz Uzx + m? /T"^., and 

^ , 4(m/T)2 - {m^/Tf 



J± 



/C 



± 

Nf{l) 
Si 
S2 



1 + 



4(m/r)2 - (m^/T)^ 



nj(E, + /i/T)+nj(E,-/i/T), 

72/(^2 + /i/T) + /i/T) + ^/(E^ - /i/T) - /i/T) , 

nf{E, + fi/T) nj{Ey- fi/T) + nf{E, - ^/T) nf{Ey + ^/T) , 



-4 



1 - 



2x1 



7 

'2 



4m^ 
3 



tan ^ 



\ 



4m2 



- 1 



tan 



£_ 





4m? 



2m? 



^ / 3 _ 
m? \2 Am? J 



log 



mi 



0:2 



"3 



— m. 



5 2 




1 


— ""^^ 




-m 


12 ^_ 




3 


■m?\ 






^ Off 




6 / 







4m 



m" 



- 1 tan"^ 



mi 



(6.17) 

(6.18) 

(6.19) 
(6.20) 
(6.21) 



(6.22) 
(6.23) 



(6.24) 



with uly = z^ + y'^ + 2zy U;,y + m//T'^, x = k/T, z = pi/T, y = P2/T, nfe(x) = [exp(x)-l] \ 
and nf{x) = [exp(x) + 1]^"*^. 

Given the perturbative expressions obtained in this section, we can compute the thermo- 
dynamic potential Qy in the OPT framework, which corresponds to resumming all dressed 
diagrams of the direct and exchange types, generating in practice an effective mass for the 
fermions. 
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6.4 Nonperturbative thermodynamic potential 

Following the standard procedure, the OPT framework [165] can be implemented in the 
Yukawa theory as 



£oPT = ip^ [i^ - (m + ?7*)] ipa + 5 



S g^ipo^il^a , (6.25) 



where a sum over flavors is implied, and rj* = rj^l — 6). As one can easily notice from the 
deformed Lagrangian above, at 5 = the theory is an exactly solvable theory of massive 
fermions even in the case when, originally, m = 0. In this particular case, the arbitrary mass 
parameter, rj, works as an infrared regulator which proves to be very useful in studies related 
to chiral symmetry breaking. We should remark that, due to our choice, the meson sector 
disappears at 5 = since there is no need for their mediation in this case. Finally, when 
6 = 1 the original, interacting theory is recovered, so that our £opt interpolates between a 
free (exactly solvable) fermionic theory and the original one. 

The Feynman rules generated by the interpolating theory are trivially obtained from the 
original ones: m h-^ m + 1]* , g 6g, and the bosonic propagator receives a. 1/6 factor. It is 
important to notice that the interpolation does not change the polynomial structure of the 
original theory and hence does not spoil renormalization, as proved in Ref. [170]. Now, a 
physical quantity such as fly^^ is perturbatively evaluated in powers of the dummy parameter, 
6, which is formally treated as small during the intermediate steps of the evaluation. At the 
end one sets 5 = 1, in a procedure which is analogous to the one used within the large- 
approximation, where A^ is formally treated as a large number and finally set to its finite value 
at the end. The result, however, depends on the arbitrary parameter, 77, which can be fixed 
by requiring that Qy^"^ be evaluated at the point where it is less sensitive to this parameter 
[171] (Principle of Minimal Sensitivity, PMS). This can be accomplished by requiring 



= . (6.26) 

r?=r?, (5=1 



dr] 

In general, the optimum value fj becomes a function of the original couplings via self- 
consistent equations, generating nonperturbative results. 

Now we can apply the method to the case of the two-loop perturbative Qy given in Eq. 
(6.6). Using the OPT replacements for the Feynman rules and expanding rj* to order 6, one 
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obtains the first nontrivial result 



K 



In [K^ - ml] + ij^ Trln [f - (m + r])] + 6 Ti 

Sdir(0,^) + Scxc(i^,^) 



where 



and 



Sdir(0,?7) 



Tr- 



[m + ?7)] 



To 



(6.27) 



(6.28) 



(6.29) 



tQ[Q-{m + r^mP-Qy~ml]- 
Notice that, when compared to the second-order perturbative result, Eq. (6.27) displays an 
extra contribution given by the third term on its right-hand side which represents a one-loop 
graph with a 6t] insertion. Now, using Eq. (6.26) one arrives at 







-z>^ Tr 



-6 Tr 



Sdir(0,r7) + Sexc(P,r/) 



+ 5 Tr 



_[/'- (m + r/)]^ 



(6.30) 



where, to obtain the last term, we have used a redefinition of momenta p i— )■ g, g i— )■ p at an 
intermediate step. Setting 5 = 1, one has a nontrivial, coupling-dependent self-consistent 
integral relation for the optimum mass parameter involving the self-energy given by 



lyj^ Tr 



7] - Sdir(O) - 'Lc^c{P,ri) 



0. 



(6.31) 



[f- (m + r/)]2 

Since the OPT propagator has an infinite number of rj insertions, one can see that the op- 
timization will resum exchange graphs in a nonperturbative way. It is very interesting to 
notice that when exchange terms are neglected f] = Sdir(O) = —g{<f))o, where (0)o represents 
the scalar condensate, satisfying the MET self-consistent relation for the effective mass (ob- 
viously, when (7 — the OPT results agree with the free gas case). This type of result is 
consistent with applications of OPT to different types of theories [172] and illustrates the 
way OPT works. Notice that within the Hartree approximation one also adds and subtracts 
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a mass term which is determined self-consistently. The basic difference is that the topology 
of this term is fixed from the start: direct terms in the Hartree approximation, direct plus 
exchange in the Hartree- Fock approximation. Within OPT the effective mass m + 77 is arbi- 
trary from the start, its optimum form being determined by the topology of the contributions 
considered in the perturbative evaluation of a physical quantity such as the thermodynamic 
potential considered here. 

In the next two sections we shall study the pressure numerically, using different approx- 
imations in several different situations. With this aim we set m = 0.1 Aj^, m<^ = m/2 in 
our numerical routines'^. The temperature and chemical potential ranges considered cover 
—7- 10m while the coupling values cover — t- tt. 

6.5 Results in the presence of a scalar condensate 

In this section we compare the results generated by PT, OPT and MFT when the scalar 
condensate represented by the direct (one-particle reducible) terms are considered. This 
case is also interesting because one can use the well-established MFT to analyze the results 
provided by OPT and PT. As already mentioned, when exchange contributions are not 
considered in OPT, MFT results are exactly reproduced since in this case both theories 
employ the same effective mass M^s = ''^ + Sdif"'"(0) = fn + fj. Each approximation considers 
different two-loop contributions: PT takes all graphs shown in Fig. 6.1 into account, while 
OPT considers all plus the extra one-loop fermionic graph with the 5ri insertion. In practice, 
as shown in the previous section, the MFT result is quickly recovered from the OPT one by 
neglecting the exchange term. Figure 6.2 shows the pressure as a function of yU at T = 0. As 
one can see, PT predicts very high values for the pressure as /i increases, in disagreement 
with MFT and OPT results. Since OPT agrees exactly with MFT when exchange terms are 
neglected, one can also see in this figure the effects of resumming exchange contributions: 
they yield slightly higher values of the pressure for increasing /i. Figure 6.3 shows the same 
situation but at a high temperature. In this case, the OPT-predicted pressure values are 

■^The value = 0, which is of particular interest in situations motivated by QCD, will also be considered 
in subsection 6.6.1. 
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smaller than the MFT ones as /i increases. 




Figure 6.2: The pressure, P/A^, as a function of /i/A^ at T = 0. Top: results from PT 
(continuous line), MFT (dashed line) and OPT (dot-dashed line). Bottom: differences, due 
to exchange terms, between OPT and MFT. 

In Figure 6.4 we analyze the pressure as a function of the coupling for low (0.5m) and 
high (5m) values of fi and T. As expected, all methods agree with the free g when 
(? — )■ 0. Also, at high T, fi values MFT and OPT results tend towards the free gas since 
one approaches the Stefan-Boltzmann limit, while PT has a completely different behavior in 
this situation. Finally, let us study the behavior of the OPT and MFT effective masses as 
functions of T and /i, as shown in Figure 6.5 for g = ir. Both quantities have a quantitatively 
as well as qualitatively different behavior at small T, /i values but, as seen in the previous 
figures, this effect does not manifest itself in the pressure, probably being compensated by 
the presence (absence) of exchange terms within the OPT (MFT). As T or /i increases 
the qualitative behavior of both effective masses becomes the same although there are still 
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0.0 0.2 0.4 0.6 0.8 1.0 

Figure 6.3: The pressure, P/A^, as a function of /i/Aj^ at T = l.OAj^. Top: results from 
PT (continuous line), MFT (dashed line) and OPT (dot-dashed line). Bottom: differences, 
due to exchange terms, between the OPT and MFT. 

quantitative differences. These results clearly illustrate how PT is not appropriate to deal 
with this kind of situation since it cannot resum the condensate which arises from the 
interaction between the scalar field and fermions and which is related to symmetry breaking. 
Our application nicely illustrates the reliability of the OPT results since they exactly agree 
with MFT at large- (when exchange terms are neglected), while allows us to improve over 
this approximation by resumming exchange contributions. 

6.6 Results in the absence of a scalar condensate 

Let us now follow Ref. [47] and neglect the one-particle-reducible two-loop diagrams con- 
taining tadpoles. In this case, which is relevant for the situation where symmetry breaking 
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Figure 6.4: The pressure, P/A^, as a function of for /i = T = 0.05A^ (top figure) and 
= T = O.SAj^ (bottom figure). PT corresponds to the continuous hne, MFT to the dashed 
hne and OPT to the dot-dashed hne. 

is not present, one expects that PT wiU perform better than in the previous case where 
the tadpole contributions have been considered, especially in the zero-temperature limit. 
Clearly MFT is not applicable in this situation, and we will restrain our analysis to OPT 
and PT results. At this stage, the reader should be convinced that, as explicitly shown in 
the previous section, the former method is able to generate nonperturbative results which 
can be used to access the eventual breakdown of PT. 



6.6.1 Cold and dense case with = 

In this subsection, the order-(5 OPT results will be compared to the second-order perturbative 
predictions of Ref . [47] . The case = is of particular interest since in this situation one 
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Figure 6.5: The OPT (dot-dashed hne) and MFT (dashed hne) effective masses (m + f/ and 
m + Sdir(0)) in units of m for g = n. Top: Mcs/m as a function of /i/Aj^ for T = 0. Bottom: 
Mes/m as a function of T/A^ for /i = 0. 

can write the thermodymanic potential fiy = fi™'^ + f2y'^'^(T = 0, /i) in a simple and compact 
analytic form, in terms of the vacuum part 



lim n™'^ 



ANFB^™(m)+ lim H 



exc,REN , 



with 



lim fi--^^^ 



N, 



2 4 

m 

1 



l^'i (0) + 



m,m^) 



(6.32) 



7 + log 



A2 



^^2 (0) + 



+ - 



7 + log 



2\ n2 



1^3 (0) + 6 a (m^) 



(6.33) 
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and the in-medium contribution 



+m u 



(6.34) 



fipf-m^ log (^^). 



where = /i^ — rn'^ and u 

Figure 6.6 shows the pressure as a function of fi for a large value of the coupling, g = n. 
As one can see both methods predict very similar results. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Figure 6.6: Pressure in the absence of tadpoles normalized by as a function of the 
chemical potential fi normalized by A^. Dots represent the OPT result and the line stands 
for PT. The fermion mass is fixed at m = 0.1 Aj^, and g = n. 

To analyze the tiny differences, let us define the quantity AP/Pp = |(Popt — Pp)\/Pp 
where Pp and Popt are, respectively, the pressures predicted by PT and OPT. Although 
the numerical discrepancies appear to be rather small. Figure 6.7 nicely illustrates that 
AP/Pp increases with higher couplings and decreases with higher chemical potential values 
as opposed to the case where the scalar tadpole is present. 



This behavior can be better understood if one analyzes how the OPT effective mass 
varies with g and fi. Figure 6.8 shows that the quantity {m + ff)/m deviates from 1 as 
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Figure 6.7: Difference between the perturbative and the OPT pressures in the absence of 
tadpoles as a function of the couphng g and the fermion chemical potential fi normalized by 
Aj^. The fermion mass is fixed at m = O.lAjjg. 

increases but, contrary to the case where tadpoles are present (see Fig. 6.5), approaches 1 
as /i increases. 

Therefore, the nonperturbative OPT results support the PT results of Ref. [47], at T = 0, 
when scalar condensates are not considered. Notice that even though we have considered 
in this subsection only the = case our numerical simulations show that the agreement 
between PT and OPT remains valid for ^ 0. For more details concerning the effects of 
the scalar mass the reader is referred to Ref. [47]. 

6.6.2 Thermal effects 

So far, we have seen through comparison with MFT and OPT that PT does not give reliable 
results when condensates are present but, as seen in the previous subsection, the situation 
improves in their absence, at least at T = 0. In principle, it is not obvious that PT will 
furnish reliable results at high temperatures, even when only exchange contributions are 
considered. The aim of this subsection is to analyze this situation. Figure 6.9 compares the 
OPT and PT results for the pressure as a function of /i in two extreme situations: T = and 
T = 10m = l.OAj^. In the first case both methods agree well, as one should expect from the 
discussion performed in the previous section, but more surprisingly is the high-temperature 
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Figure 6.8: The OPT effective mass m + f/ in units of m as a function of the couphng g and 
the fermion chemical potential /i normalized by when tadpoles are absent. The fermion 
mass is fixed at m = O.lAj^. 

result where OPT also supports PT. 

Like in the previous T = case, the quantity AP/Pp assumes very small values which 
increase with high values of g but decrease with high values of T or /i. The OPT effective mass 
in units of m, (m + f/)/m, also behaves as in the previous case, deviating from 1 a.s g increases 
and, contrary to the case where tadpoles are present (see Fig. 6.5), approaching 1 as /i and/or 
T increase. Therefore, it looks like the absence of direct terms means that m ^ m + f/ at 
high T and/or fi, so that PT behaves like OPT. Finally, Figure 6.10 suggests that even at 
relatively high T and fi OPT and PT deviate from the free gas as the coupling increases, 
while in the case with condensates (see Fig. 6.4 ) OPT (and MFT) had a better agreement 
with the free gas as opposed to PT. Therefore, the correct resummation of condensates seems 
to reduce the effects of interactions at high T and fi. 

6.7 Final Remarks 

This chapter represents a first step towards a nonperturbative picture of the phase diagram 
of Strong Interactions via (chiral) effective models. We have investigated the thermody- 
namics of the Yukawa model at finite temperature and chemical potential by evaluating its 
thermodynamic potential up to the two-loop level which includes direct (Hartree-like) as 
well as exchange (Fock-like) types of contributions. Three different methods have been con- 
sidered: the usual perturbation theory (PT) with its bare mass, the optimized perturbation 
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Figure 6.9: The pressure, P/A^ as a function of /i/Aj^ at T = (top figure) and at 
T = 1.0 when tadpoles are absent. PT results are represented by the continuous lines 
while the OPT results are represented by the dot-dashed lines. 

theory (OPT) with its effective mass given by the PMS variational criterion, and the well- 
known MFT with its self-consistent effective mass. Our comparative analysis shows that, as 
far as thermodynamic results are concerned, OPT is not only consistent with MFT when 
one-particle irreducible contributions are neglected, but also furnishes nontrivial corrections 
when the exchange class of diagrams is included. 

In this vein, OPT seems to be a well-defined framework that accommodates all the desir- 
able ingredients to obtain a robust description of phase diagram of Strong Interactions via 
effective theories. It suffers, however, from its close connection to perturbative computation 
techniques and their renormalization procedure. In the case of the chiral phase transition 
described within the LSM, for instance, two-loop vacuum contributions apparently spoil the 
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Figure 6.10: The pressure, P/A^ as a function of g ioi ji/ = T = O.OSAj^ (top figure) and 
at fi = T = 0.5Aj^ wlien tadpoles are absent. PT results are represented by the continuous 
lines while the OPT results are represented by the dot-dashed lines. 

positivity of the effective action of the model (cf. e.g. [48, 173]) ^ and further work to 
circumvent this obstacle is called for. 

Our results are nevertheless general, featuring coupling values ranging from g = to 
g » 1 in situations where the temperature and chemical potential ranged from zero to 
ten times the highest mass value, which is kept fixed. As discussed in the Introduction, 
the Yukawa model usually emerges in the description of various physical situations, ranging 
from low-energy condensed matter phenomena to extremely energetic QCD matter. In our 

"'This issue is actually the analogous of the Higgs potential instability in the Standard Model of particle 
physics, which amounts to the requirement of a finctuning of the Standard Model Higgs mass. In the case 
of Strong Interactions, however, all the physical scales related to spontaneous chiral symmetry breaking are 
known and a similar finctuning is therefore forbidden. 
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approach, the characteristic features of each physical system will be brought about essentially 
by the specific values of the coupling g and the three energy scales: the fermion and scalar 
masses m and m^, respectively, and the renormalization scale Aj^, which normalizes all 
the quantities in our plots. As observed previously using the perturbative method [47], 
variations of the masses can significantly affect the thermodynamic potential of the Yukawa 
theory, yielding extremely different thermodynamical pictures. The renormalization scale 
^Ms ^^^^ typical energy scale of the system of interest. In hadronic physics, for example, 
it is reasonable to choose Aj^ ~ 1 GeV, the confinement scale. In this chapter, we kept 
the discussion in general grounds, normalizing physical quantities by Aj^ and fixing the 
masses as m = 0.1 Aj^ and = , 0.5 m, and concentrating on the effects of Hartree- and 
Fock-like interactions and nonperturbative corrections. 

First, we have analyzed the pressure with both direct and exchange contributions, with 
the former being associated to the presence of a scalar condensate driven by the interactions 
with fermions. We have shown that OPT and MET are identical if one does not consider 
exchange terms, in agreement with many other applications [172], which is reassuring since in 
the limit of direct contributions only the MFT resummation can be considered as "exact". 
Moreover, nonperturbative effects of exchange contributions are readily incorporated by 
considering OPT consistently up to two loops. In the light of the nonperturbative approaches 
OPT and MFT, our results show that, as expected, naive PT is inadequate to deal with this 
situation since it has no ability to resum tadpoles. 

As a byproduct of this application we could see how the resummation of exchange terms 
performed by OPT corrects the MFT framework, which corresponds to the leading order 
of a 1/A^ type of approximation. Having established the reliability of the OPT, we have 
followed Ref. [47] imposing the absence of tadpoles at T = 0. We have then shown that, in 
this case, the OPT results turn out to be very similar to the ones given by PT. Finally, still 
in the limit in which condensates are not present, we have investigated the high-temperature 
case where one could expect the breakdown of PT. However, our findings have shown that 
this is not the case and the numerical differences between the OPT and PT pressures are 
very small. The outcome of this analysis also suggests that, at high T and /i, the presence of 
condensates minimize the effects due to interactions when these contributions are properly 
resummed. 
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One should notice that the direct apphcation of these results as effective-model predic- 
tions in different physical contexts is restricted to definite energy regimes, within which the 
relevant physical degrees of freedom can be translated into a Yukawa model. For instance, 
when discussing the thermodynamics of cold and dense baryonic matter, the Yukawa model 
considered here might be a suitable effective theory only at small values of the chemical 
potential, e.g. at < 350 MeV. To describe the properties of baryonic matter at higher 
values of /i, a different framework is necessary to account for the phenomenon of color su- 
perconductivity (for a review, see Ref. [13]). One alternative could be to extend this model 
by adding boson fields from another, non-singlet, multiplet of the color SU{3) group. Never- 
theless, since Yukawa-type interactions are almost ubiquitous in the description of fermionic 
matter, this analysis provides an understanding of the interplay between direct and exchange 
contributions to the thermodynamics as well as the role played by nonperturbative effects 
in the scalar Yukawa sector of any given extended model. 



Eff. models j 

beyond 
mean-field 



Chapter 7 




NONPERTURBATIVE ANALYSIS OF 
RELATIVISTIC BOSE-ElNSTEIN 
CONDENSATION OF PIONS 

In this chapter we will address the phase diagram of Strong Interactions in the temperature- 
isospin chemical potential plane through a nonperturbative approach. 

Instead of analyzing the deconfining phase transition in this type of media, as done in 
Chapter 2, we will neglect baryon density effects and concentrate on the phenomenon of rela- 
tivistic Bose-Einstein condensation of pions to investigate the influence of (nonperturbative) 
interactions in this context [174]. Besides its phenomenological motivation, the problem of 
relativistic BEC in isospin dense QCD represents a particularly interesting framework to 
develop efficient nonperturbative methods for dense systems: being Sign-Problem free, a 
robust reference from lattice simulations at finite isospin chemical potential is in principle 
available. Moreover, the role of interactions in the BEC phenomenon is a nontrivial funda- 
mental question. Even for the nonrelativistic case, some issues were only understood quite 
recently, such as the shift of the critical temperature in cold atomic systems, that was suc- 
cessfully described using the Functional Renormalization Group framework [175]. Our intent 
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in this on-going work [174], centered in the relativistic version of BEC that should occur in 
QCD matter, is to analyze in detail how the Functional Renormalization Group technique 
and the different possible approximations that must complement it perform in this case. 

7.1 Basics of BEC in isospin dense QCD 

Proposals of pion condensation in Strongly interacting matter under extreme conditions 
were made back in the 70's by Migdal and others (for a review of these early works, cf. Ref. 
[176]). They were basically concerned with astrophysical environments and concentrated on 
electrically neutral media. The imposition of electric neutrality actually lead them to another 
(region of the) QCD phase diagram than the one to be addressed in this chapter: not only 
isospin density is present, but also baryonic. In this more complicated and phenomenological 
context, the pion condensate is found to be inhomogeneous and triggered by complicated 
interactions with nucleons and external fields. 

Since we are interested in developing nonperturbative techniques for dense systems, we 
shall restrict our approach to isospin-dense, baryon-free matter for which the Sign Problem on 
the lattice is absent Pioneering work on this type of system was done by Son and Stephanov 
(cf. e.g. [23]). Pion condensation in this case will be a clean relativistic realization of Bose- 
Einstein condensation and, differently from these early proposals, the condensate occurs as 
a constant and uniform state whose onset is triggered even for weak coupling. 

Indeed, pion Bose-Einstein condensation in isospin dense environments is such a clear 
relativistic manifestation of thermodynamical properties of bosons that one can formulate a 
simple heuristic argument for its existence. Departing from the thermodynamical definition 
of chemical potential associated with a given net charge: the free energy the system spends 
(gains) when one negative (positive) charge is taken out of (added to) the system. In the case 
of isospin-dense media, charged pions vr^ may play the role of ±1 isospin charges, but, being 
massive, the relativistic free energy spent when including one of them in the system must 
account for its mass rrij^, at least. It is then straightforward to conclude that the free energy 

^Actually, one must consider two-flavor QCD with degenerate masses m„ = irid, for which the imaginary 
part of the fcrmionic determinant cancels exactly. 
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cost to add a zero- momentum vr^ decreases with increasing isospin chemical potential /i/, 
going to zero when it reaches the pion mass. At fij = an instability sets in and the new 
thermodynamically stable ground state corresponds to a macroscopic occupation of vr+'s in 
the zero- momentum level, i.e. a Bose- Einstein condensate of 7r"'''s. 



7.1.1 Modeling in a chiral effective theory 

Let us now build up an effective framework to address the issue of Bose-Einstein condensation 
of pions in isospin dense media. We will show that the heuristic argument presented above is 
indeed verified in a chiral model in the form of a nontrivial static and uniform pion condensate 
at high isospin densities. 

Departing from a standard chiral Hamiltonian, the constraint of isospin number conser- 
vation is imposed via the method of Lagrange multipliers: "H t— t- "H — fijui, where the isospin 
chemical potential /i/ plays the mathematical role of the Lagrange multiplier for the con- 
servation of the isospin number density rij. Therefore, to describe an isospin dense system 
within a chiral effective theory, it is crucial to define isospin number nj in the context of 
chiral transformations and fields. 

Isospin is the conserved charge associated with the invariance of Strong Interactions under 
a two-flavor subset of the global vectorial chiral transformations exp —i6 ■ Tap/ 2 , which are 
left unaffected by the spontaneous chiral symmetry breaking. Here, r = (ri,r2,r3) are the 
Pauli matrices, a, /3 are indices in flavor space and 9 is the parameter of the transformation. 
The specific isospin subset is the one related to the light flavor sector of quarks up and 
down. We will be interested in the projection of the isospin on a specific coordinate axis, 
say z, so that the conserved isospin current can be obtained, through Noether's theorem, by 
transforming the Lagrangian locally using 9 = (0, 0, 03(x)). In QCD, one derives in this way 
the conserved isospin current in terms of quarks 

fv = ^.7^^*/., (7.1) 

while in a chiral effective theory involving mesons (such as the LSM or, equivalently, the 
0(4) model), one obtains the isospin number rii = 113 {n = jy) related to the pion degrees 
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of freedom vf = (tti, 7r2, tt^): 

Tij = 7ri9°7r2 - TTsS^TTi =i[Ti-d'^n+ - n+d^n-] , (7.2) 

where we have used the mass-eigenstate basis for the pions, vr^ = (tti =i=i7r2)/ a/2 and = ir^, 
the latter being neutral under isospin rotations around the z — axis. A peculiar feature of 
the isospin number of pions is that it involves not only the charged pion fields but also their 
time derivatives, i.e. the momenta canonically conjugated to these fields. 

Since we are ultimately interested in thermodynamics and the phase transition to the pion 
condensed phase, we need to construct the partition function and the associated euclidean 
action for this isospin dense chiral model. For the LSM without constituent quarks, defined 
in the vacuum by the Lagrangian (the reader is referred to Appendix A for an extended 
description and references): 

>C = l{d,a){d^a) + l{d,^){d'^n)-V{[a' + n%a) , (7.3) 

{[a' + n%a) = ^{[a' + n']-v'Y-ha, (7.4) 
the result is (a detailed derivation is given in Appendix G; cf. also Ref. [149]): 

[D$] e"^^ , (7.5) 
with the euclidean action for the isospin-dense LSM being given by: 

Se = j d'^x!^^{df,a){d^a) + ^{d^7f){d^7f) + /i/fvriSoTra - tts^otti] - ^fi'j (tt^ + vr^) + 

+V,i{[a^ + n%a)]. (7.6) 



Within this model, pion Bose-Einstein condensation is realized in the form of a nontrivial 
pion vacuum expectation value at high isospin densities. To verify this and understand the 
quasiparticle degrees of freedom in both phases of this isospin-dense system, one should 
consider the theory in the presence of condensates and find the solution that minimizes the 
effective action as the isospin chemical potential fij is varied. 

Adopting, as a first step, the classical (or mean-field or Landau- Ginzburg) approximation 
to obtain the condensate solutions for different fij domains, the effective action is assumed 



CHAPTER 7. NONPERTURBATIVE ANALYSIS OF RELATIVISTIC. 



135 



to have the same form as the euchdean action in Eq. (7.6). The equihbrium expectation 
value of the fields is then given by the absolute minimum of this action. 

Considering only static uniform condensate solutions (a, n) and small perturbations 
around it {C,,if), 



a = a + C, 
n = n + f], 



(7.7) 
(7.8) 



the effective action in the classical approximation can be expanded as: 



d'^ lkd,Oid,0 + kd.ff) ■ id.ff) + fii [(vr2 + V2)doVi - (vn + Vi)doV2] + ^r''^ 



Sf 



while the expansion of the potential, defined in Eq.(7.4), reads (N.B.: n = (vti, 7r2, tts) 



V2 



V2 



(Lo-M) 
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4 1 
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+ViVj 
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{2cx' 2tt, + 2tt' 27r, + Awii-v')) - ^(1 - ^is) [2vr. 



+ 



L 4 



(4^2 + 2a^ + 2(^2 _ ,,2)) 



+ 



A^ 



A^ 



L 4 
A2 



{ATTiTTj + 2a^5ij + 27r^(5ij + 2[-v^]6ij) 



/if 



2 ^ii(l -fe) 



L 4 



(4(T27r, 



+ t{^' + ViV^VjVJ + 4ae[e' + f ] + 2e'[27f . rf + f ] + 4^ ■ r/[f ]|(,7.io) 

where: (i) the first line corresponds to the classical potential evaluated at the condensate; 
(ii) the second and third lines (linear in the fluctuations) give the equilibrium conditions: 








1,2,3) 



(7.11) 
(7.12) 



(iii) the terms quadratic in the fluctuations (4th, 5th and 6th lines) are related to the 
dispersion relations of the quasiparticles above in the condensed phase and define the stability 
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of the condensate solution; and (iv) the last line corresponds to the interaction vertices for 
the fluctuations. 

The dispersion relations for the quasiparticles in each phase are derived in Fourier space. 
Using the following field expansion for each component of the fluctuation fields {ip = C)'?*)- 



(7.13) 



Eq.(7.9) then can be written as (up to quadratic order in the fluctuation fields and using the 
equilibrium conditions, Eqs. (7.11) and (7.12), to cancel the terms which are linear in the 
fluctuations): 



a, TV 



+Pyi[m ViiQ) V2iQ) r/3(g) ) yA^'(Q) 

J (J 



+0(3), 





-Q) ^ 




-Q) 




-Q) 







(7.14) 



with: 



M\Q) 



-Ql + {q)' - 



uric 



-Ql + {(i? + 2A27ri7r2 + 2z/i/Qo 



2X'^7llTT3 



2A27f27f3 



2X^a7t3 ^ 

2\'^7llTT3 

2\^n2Ti3 
-Ql + {q)^ + TA j 

(7.15) 



where we have defined: 



rric 



X^[27r^ + P + a' 



(7.16) 



We have also used the notation 0(3) to represent the interaction between the perturbations, 
i.e. the terms that are cubic or quartic in (^,?7j). 
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The spectrum of quasiparticles is then given by the energies Qo which correspond to a 
zero of an eigenvalue of the mass matrix Ai'^{Q) in the hmit of zero three-momentum (i.e. 
a pole of the propagator). Therefore, they are solutions of the following equation: 

det[M\Qo,q = 0)] =0. (7.17) 



In what follows, from the equilibrium conditions, Eqs. (7.11) and (7.12), the conden- 
sate solution is obtained as a function of the chemical potential: ((T(/i/), -^(/i/)) . Then the 
quasiparticle spectrum is derived through Eqs. (7.15)-(7.17) in each phase. 

Consider the (physical) case with explicit chiral symmetry breaking {h = fT^m^ ^ 0). 
Since h 0, Eq. (7.11) implies that there exists a nonzero chiral condensate, in the 
(T— direction, i.e. a ^ 0, and also that the combination ((x^ + n"^ — v"^) must be nonvan- 
ishing. This result together with Eq. (7.12) implies that: 

7r3 = 7r° = 0. (7.18) 

On the other hand, the condensates vti and are not constrained similarly and there 
will be solutions of the equilibrium conditions for vfi and Ti2 both vanishing or not: 

• Normal phase, with vanishing pion condensates (vfj = 0, for all i): The equilibrium con- 
ditions in the pion directions are trivially satisfied, while in the cr— direction we get the 
isospin-independent condition: 

\^cx[cx^ -v^)-h = {). (7.19) 

This equation has in general 3 solutions, which will be real and distinct for the physical 
values^ of ~ 20, v"^ ~ 7700 MeV^ and h = fT^m^. The true minimum is associated 
with the solution a > f , while the others correspond to a metastable minimum and a 
maximum in between the minima. 

The parameter fixing of the model implies that the true minimum of the potential in 
the vacuum (/i/ = 0, in this case) occurs at a = f^^ (cf. Appendix A). Therefore, Eq. 
^The procedure of parameter fixing within the LSM is discussed in Appendix A, Subsection A. 3. 
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(7.19) becomes a condition^ on the parameter v: 



x'fAf'-v')-h = o 



J n 



J TT 



m: 



(7.20) 



The value of the classical approximation for the true minimum of the effective potential 
in this case reads: 



(Lo-M) 



cl 



A 



min,7r-^=0 



a' 



4A2 



.212 



J TT 



(7.21) 
(7.22) 



In this normal phase (/i/ < m^), the mass matrix in Eq.(7.15) becomes diagonal except 
for the terms ~ f^iQo, and the determinant equation for the spectrum becomes: 

i-Ql + ml){-Ql + ml) {{-Ql + ml)i-Ql + ml) + 4^V?Qo} = (7.23) 



Qo = \m^\ or Qo - Qo("^i + "^2 + 4/i/) + 



2 2 

mi 



0(7.24) 



Using Eqs. (7.16) and the condensate solution (cr = yf = 0), the resulting spectrum 
has the following masses: 

Qi""^ = \m^\ = .j2\^fl + ml = m,, 



(7r±) 



'm{ + m2 + 4/i^)2 — 4m^m2 



f^i T 7T\/4(m2 + /i2)2 _ 4(m2 - fij)^ 



(7.25) 



Pion-condensed phase (tt^ = (tt^ + tt"^) 7^ 0): Let us define radial and angular coordi- 



nates such that 



#2 



7f2 



(vri' + vr2^)=pV0 

pcosa 

p sin a . 



(7.26) 
(7.27) 
(7.28) 



■^Strictly speaking, one should first verify tlrat tlris condensation case (and not one witli e.g. tt^ 7^ 0) 
indeed minimizes tlie effective potential in the vacuum. This can be seen from the equilibrium conditions, 
Eqs. (7.11) and (7.12), with ^7 = 0: they cannot be solved simultaneously if /i 7^ and ^/ = unless we 
consider Tfi = 0. 
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The equilibrium conditions then become: 



X^a{(T^ + - v^) - h = (7.29) 
{\'^{a^ + -v^) - f4} pcosa = (7.30) 
{\'^{a^ + p"^ -v^) - ^j} psina = 0, (7.31) 



which are satisfied for vr^ = 7^ if and only if: 



\'^a{a^ + p"^ -v^) -h = (7.32) 
X\a^ + p^ - v^) - fij = 0, (7.33) 



so that 



^2 



(7^0) p\-ml ,2_ff_ 



A2 



4 



(7.34) 



It should be noted that the equilibrium conditions are completely independent of the 
angle a (as the remaining symmetry under rotations around tt^ requires). Physics 
does not depend on the direction of the condensate in the tti — 712 plane, which is 
defined by a choice of the coordinate system. If we choose the direction of the pion 
condensate such that sin a = (cos a = 0), the condensate occurs in a direction that 
is a linear combination of a and vti {1^2), while vr^ mass rises linearly with the chemical 
potential and the remaining pion stays massless at any pj for which this solution is 
the true equilibrium state of the theory. Independently of the value of a we follow the 
literature ([177] and references therein) and name the quasiparticles in the condensed 
phase by a continuity criterium in the spectrum. 
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The minimum of the effective potential in the classical approximation in this case is: 



(Lo-M) 



min,7r2^0 



1 



4A2 



2 H 



ml 1 /ij 1 



4A 

(7^2) .,{LaM) 



9 2 



(7.35) 
(7.36) 
(7.37) 

- 2 

(7.^^8) 



so that it is clear that it is smaller than K 



(Lo-M) 



min,7r-^=0 



for Hj > rriT,. For /i/ < m^, 



p < 0, the solution is unstable in the directions of the tti and Ti2 fields, as can be 
verified by examining the quadratic terms in the expanded potential, Eq. (7.10). 

For this phase with pion condensation and sin a = 0, the Eq. (7.17) for obtaining the 
quasiparticle spectrum becomes: 



(Qo-'^D(Qo) h^H + (2^'^/')'-'^|(2/^/)' + QoK + '^? + (2/i/)') -QS = 0, 

(7.39) 



whose solutions are: 



Qq = \/"^3 = l"^3|, 

Qo = 0, 



(7.40) 
(7.41) 



together with the solutions of the quadratic equation: 

-mlml + {2\^apf -ml{2^iif + Ql[ml + ml + {2^iif) -Ql = 0, (7.42) 



namely: 



Ql 



^ |m2 + ml + (2/i,)2 ± (-m| + m? + (2/i,)2)' + 4(2A2(Tp)2| (7.43) 
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Using the results for the condensate solutions, Eqs. (7.34) and (7.16), we arrive at: 



(a) 





-ml + X^fl ± 



(7.44) 
(7.45) 



±1 



p/ 



At the transition (/i/ = m,r), these solutions continuously connect to the particle 
spectrum in the first phase (Eqs. (7.25)). We will, for the sake of definiteness and 
connection with the literature, use the same names for the quasiparticles in the second 
phase, assuming continuity: (a) i— )■ vr", (b) H- vr"*", (c, +) H- a and (c, — ) H- vr^. 
Nevertheless, it should be emphasized that the states in the second phase are (apart 
from the neutral pion, which decouples) nontrivial mixings of the sigma and the charged 
pion fluctuation fields. 

The results for the mass spectrum in both chiral and pion condensed phases are summa- 
rized in Fig. 7.1. In the normal phase, for low densities, a and 7r° are independent of fij, due 
to their invariance under isospin rotations around the 2;-axis. The free energy cost related to 
one charged pion, however, changes linearly with fij. As discussed heuristically above, when 
Hi = m^r = 138 MeV, the vr"^ becomes effectively massless and its bosonic character causes a 
Bose- Einstein condensate to develop at higher densities. The phase transition is therefore of 
the second order. Due to the effective isospin rotation the ground state of the system suffers 
for fij > rriT^, originating the nonzero pion condensate, all the fluctuations on top of it have 
clearly a nontrivial isospin dependence, as shown by the dashed lines. 



(7.46) 
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Figure 7.1: Full lines are the masses in the normal phase, with standard chiral condensation, 
corresponding to, from top to bottom: a, tt^, 7r° and tt"^. The dashed lines are the masses 
of the quasi-particles of the phase with pion condensation. 

7.2 The Functional Renormalization Group 

Before applying the Functional Renormalization Group (FRG) to investigate nonpertur- 
bative aspects of relativistic Bose-Einstein condensation of pions, let us discuss in more 
fundamental grounds the formalism itself. This section is not intended to be a thorough 
presentation of the subject, but rather a brief introduction to the main ideas and hypotheses 
that form the basis of FRG. Along the discussion, we shall demonstrate and/or justify the 
main results that will be used in what follows. For comprehensive texts the reader is referred 
to [178, 179, 180]. 

The main aim of the formalism of the Functional Renormalization Group or non-Perturbative 
Renormalization Group is to provide a systematic framework for the inclusion of nonper- 
turbative effects in field theory calculations. This is implemented via the construction of a 
set of effective actions Tj^cj)] which interpolate between the classical action S[(j)\ and the full 
effective action r[0]. 

Within the FRG formalism, this set of interpolating effective actions is constructed in 
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the form of a flow, i.e. a trajectory in the effective action space parameterized by the scale 
K. It is determined via an exact renormalization group equation for the effective action T^, 
which actually encodes an infinite hierarchy of coupled exact renormalization group equations 
involving the ra— point functions. 

The full solution of such a hierarchy is formally equivalent to solving the full interacting 
quantum field theory and as such is prohibitively complicated in the great majority of the 
physically interesting cases. An adequate truncation or approximation is needed to render the 
exact FRG flow a practical nonperturbative method. By adequate we mean an approximation 
that contains the relevant physical features of the system under investigation. Some of the 
most popular truncations and approximations are critically discussed in subsection 7.2.2. 

In principle any k— dependent deformation of a given theory defines an exact flow in 
the effective action space. Nevertheless, for any particular flow to be useful as a tool for 
computing nonperturbative effects some special features are required (or desired). Among 
them a smooth k — > limit and simplicity are relevant for pragmatic applications. The 
deformed theory at finite k, should partially include the nonperturbative corrections, so that 
the K — > theory encodes the full nonperturbative result. 

In field theory, nonperturbative effects are usually connected to the extremes of momen- 
tum values: ultra-violet (UV) or infra-red (IR) modes. The UV regime is well-understood 
and the theory of renormalization provides a systematic method for taking these nonper- 
turbative effects into account within the perturbative expansion of a renormalizable theory 
[145]. The IR regime on the other hand seems to pose more severe difficulties for a pragmatic 
definition of a computable method, even at weak-coupling (as it happens, for instance, in 
thermal perturbation theory [149]). It is reasonable therefore to define an exact fiow that 
incorporates in a controlled and successive manner the IR modes. 

With this in mind, the interesting dependent deformation of the theory should affect 
differently UV and IR modes, including the effects of lower and lower momentum modes as 
the K — limit is approached. The freedom of choice of the form of the deformation allows 
for choosing the simplest one, affecting only the quadratic part of the theory, or, equivalently, 
its propagator. The simplicity of the modification of the theory will translate into the exact 
FRG fiow equation to be solved in the end. A quadratic deformation guarantees a one-loop 
structure in the final fiow equation, as we shall verify in the next subsection. 
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Putting these features together, we conclude that one possible interesting definition of 
an exact flow is the one deflned as the original theory with a modified propagator which 
suppresses IR modes with momentum g < k. Since the aim is that of interpolating between 
the classical action S[(j)] and the full effective action r[0], the trajectory of the FRG flow in 
effective action space should satisfy the following features at different scales k: 



K = 


all quantum ffuctuations are included {Ti^^q = F) 


Intermediate k, 


IR modes (g^ < k^) suppressed; UV fluctuations (g^ > k^) included 


K = A — i> oo 


all fluctuations are suppressed; physics is classical (rK=A = S) 



Notice that the existence of such trajectory for a flnite UV scale A relies on a nontrivial 
hypothesis: the classical action must be the full effective action at some scale A, i.e., there 
exists a scale A at (or above) which all fluctuations are indeed suppressed in the full the- 
ory. Strictly speaking, this is in general not the case in quantum fleld theories. From the 
experience with perturbative UV renormalization, we know that not even renormalizable 
theories present this exact property. When one renormalizes the theory, integrating out the 
UV modes with wavelengths k, G [A, oo}, it amounts to considering the same operators as in 
the original theory, but with dressed parameters. Stated in a more precise way, the initial 
effective action Fk=a for the FRG flow is not the naive classical theory, but rather an action 
containing the same operators as the classical one with dressed coefficients. Some remaining 
fluctuations with wavelengths k > A are encoded in the dressing of the parameters of the 
classical action. 

All the subtleties and intricacies of the perturbative renormalization procedure are there- 
fore implicitly included in the FRG formalism through its initial conditions. Once flnite 
initial conditions are given and a sufficiently large initial scale A is chosen, all results at the 
end of the flow {n = 0) are both flnite and consistent. 

7.2.1 The exact flow equation 

Let us now illustrate the derivation of the exact flow equation that generates an interpolating 
trajectory as described above by considering a scalar fleld theory. The general guideline, 
however, is universal: starting from the classical deformation of the original theory through 
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a K-dependent quadratic term, one constructs (via a standard Legendre transform, including 
a convenient modification) the full effective action for the theory at the scale n. The exact 
FRG equation is then obtained straightforwardly by taking an ordinary derivative of this 
definition with respect to the flow scale k. 

In euclidean space, a general scalar 0{N) model with ($^)^ interaction is defined classi- 
cally as: 

^ = l{d,m,^) + l^'^' + ^{^')\ (7.47) 

where the field is $ = {ipi, . . . , (Pn)- 
In momentum space. 



*(^) = yiM^ I e^'-'^q) (7.48) 



6{q) = I d'^xe"'-^ (7.49) 
so that the euclidean action reads: 



^[$] = j d'^xC 



V' 1 
V^2 



q 



+ 77i7rfTy / 5{qi + q2 + qz + qi)^i{qi)Vi{q2)^j{q-i)^j{qi)- (7.50) 



As discussed above, the modified theory is defined through the addition of a quadratic 
term to the action in order to suppress IR modes and leave UV modes unaffected, the 
threshold between these regimes being given by the flow parameter k. These features are 
implemented by adding a term to the action S Sf, = S + AS'^, with: 

r2 



= (7-51) 



in which we wrote a general form without explicit 0{N) symmetry and the regulator R]^{q) 
implements the desired suppression of IR modes with g < k. Thus, typically, R^^^q) — 
when q^ <^ K,"^, and R^i{q) — )■ when q"^ ^ k"^- It is also usually convenient to use a symmetric 
regulator: R^J{q) = 5^^RK{q). 
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The generating functional of connected Green's functions of the k— modified theory is 
then given by: 



W^[J] = log j exp - AS,[$] + j d^x^x)J{x) 



(7.52) 



Following the standard procedure, one can construct from via a Legendre Transform, 

the effective action of the k— modified theory: 



J^(x) 



f.[0] = -W,[J^] + j d^xct^J^i^x) (7.53) 



' 6J{:x) 

This functional has interesting properties, namely: 



($(x)).. =0. (7.54) 



in the absence of external currents, TfJ\(f)\ is extremized by the vacuum expectation 
value of the field = ($), where the average is defined using VF^, i-e.: 



= 0; (7.55) 

J=0,(/,= ($) 

• rfj[0] is the generating functional of one-particle irreducible correlation functions; 

• In the limiting case in which fluctuations are suppressed (classical limit; Landau- 
Ginzburg approximation), this effective action reduces to the original action of the 
K— modified theory: S + A^k- 

However, as discussed above, the effective action we are interested in is actually one 
which interpolates between the full, exact theory at k = and the original classical theory 
at large scales k = A, i.e. Tpjycj)] — )■ 5(0]. While the first property is fulfilled by (since 
AS'k=o = 0), the second is not. In this vein, it is convenient to define a Modified Legendre 
Transform, subtracting the extra term AS^ to reproduce the desired classical theory when 
fluctuations are suppressed without affecting the full quantum theory (the one at k = 0): 

TM = -W^[J^] + J d''x<j)J^{x)-AS4<P], (7.56) 

and Jfj, is still the source for which the vacuum expectation value of the field equals 0, 
determined via Eq. (7.54). 
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The exact Renormalization Group equations satisfied by can then be derived directly 
from this expression for the Modified Legendre Transform. 

Taking the derivative of Eq. (7.56) with respect to k, we obtain: 



(7.57) 



5Js 



(7.54) 



{dj^) + j d''x(l)d^{J^{x)}-d^/\SM 

(7.58) 



Using the definition of the generating functional of connected Green's functions TV^, Eq. 
(7.52), the right-hand side of the flow equation above can be rewritten as: 

(-a,A54<l>])exp {-^[<l>] - A54<l>] + jd''x^{x)J^{x)] 



Y2/d 2 



(7.59) 



where we have used the quadratic form of the (diagonal) regulator term (cf. Eq. (7.51)) and 
the average represents: 

/D$ (. ..)exp{-5[$] - A5«[$] + j d^x^{x)J^{x)] 



so that {(Pi)j^ = (f>i (cf. Eq. (7.54)) and it becomes clear that: 



(7.60) 



5Jd> 



n -1 



Tr V 



^J<f>,i 



(7.61) 



where the full propagator of the k- modified theory appears^: 



d-2 
d 



6(f)i6(f)j 



6(f)i6(f) 



+ SijR^{q) 



(7.62) 



^The functional derivative notation is built in Fourier space, in which all fields are dimensionless by 
construction. Thus the overall dimensional factor — \fV (d = 4) that appears. At finite temperature, 

V~3~ H- > in our notation. 
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The final form for the exact FRG flow equation for the effective action is therefore: 



which can be viewed as a simple one- loop bubble diagram with an insertion d^RniQ) and the 
line representing the full propagator at the scale k. 

We have stated above that in the FRG framework renormalization is encoded in the 
initial conditions of the flow. Once initial parameters at the large scale k = A are given, 
the physical results, obtained in the end of the flow, are all finite. This is realized by the 
finiteness of the flow equation itself. Indeed, in Eq. (7.63) both IR and UV divergences 
are absent. The good IR behavior is guaranteed by construction due to the presence of the 
regulator in the full propagator maintaining the modes in the deep IR with a large mass 
~ K. Another property of the regulator is responsible for the taming of UV divergences in 
Eq. (7.63). The insertion d^Rniq) in the flow equation plays the role of an effective UV 
cutoff, due to the requirement that the regulator leaves the UV modes unaffected. 

The simplistic appearance of the FRG flow equation in Eq. (7.63) is deceiving. Fully 
solving this exact flow equation is equivalent to solving the whole interacting quantum field 
theory under consideration. This is of course a prohibitively complicated task in most 
interesting cases. To apprehend the complexity encoded in Eq. (7.63) it is crucial to note 
that and are the exact effective action and propagator of the full quantum theory at 
the scale k. Therefore to solve the equation for F^ one must first know the full propagator 
exactly. One can also obtain a flow equation for from Eq. (7.63) by taking two functional 
derivatives with respect to the fields. This procedure generates an equation containing the 
full 3- and 4-point functions of the quantum theory on the right-hand side, which are of 
course not known a priori. Proceeding in this line, one concludes that the compact FRG 
flow equation shown in Eq. (7.63) represents actually an inflnite hierarchy of coupled flow 
equations for all n-point functions of the full quantum theory being analyzed. 

7.2.2 FRG as a practical tool: truncations and approximations 

As it stands, the flow equation in Eq. (7.63) is a formal result and cannot be considered a 
nonperturbative method to tackle interacting physical systems. An adequate truncation or 




(7.63) 
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approximation is needed to render the exact FRG flow a powerful nonperturbative tool. In 
practice, this amounts to transforming this infinite hierarchy of complicated flow equations 
into a small set of sufficiently simple equations. Such a drastic procedure should in principle 
be implemented through the knowledge of the relevant properties of each system under 
investigation, in order to obtain a sensible physical framework in the end. 

There are two mainstreams (cf. e.g. [178, 179, 180]) of simple-minded approximations in 
the literature that are often combined: the derivative expansion and the truncation of the 
hierarchy of equations for the n— point functions by setting them to zero for n > no, with 
uq small (usually 4 or 6). 

The derivative expansion corresponds to a specific hypothesis for the form of the k- 
modified effective action rK[{0(x)}]. It assumes that an expansion on gradients of the fields 
is convergent and well-defined. The most popular particular case of derivative expansion is 
the so-called Local Potential Approximation (LPA), which correspond to the leading-order 
term. In LPA, the functional effective action is reduced to an effective potential, which is 
an ordinary function of constant and uniform classical fields: rK[{0(x)}] ~ {(p)- The 
flow equation therefore becomes an ordinary differential equation involving derivatives with 
respect to k, and the values of the classical fields {</)}. 

In many applications in the literature one simplifies further the LPA ansatz by truncating 
the taylor expansion of the LPA effective potential with respect to the fields. In this way, 
the FRG flow is transformed into a finite small set of coupled ordinary differential equations 
for the coefficients of the potential. 

Though seemingly a very drastic simplification of the full FRG flow, the LPA ansatz has 
been shown to encode different nontrivial contributions, coming from the resummation of en- 
tire classes of diagrams in the perturbative framework. For example, it can be demonstrated 
[180, 181] that the exact diagrammatic computation in the large-N limit of n-point func- 
tions of the K— modified 0{N) field theory satisfies the corresponding FRG flow equations in 
LPA. The large-N approach is well-know as an alternative to the weak-coupling perturbative 
expansion, bringing about nonperturbative effects which are not present in the latter. The 
same is therefore true for the LPA version of FRG. 

Nevertheless, for physical phenomena in which interactions with nontrivial momentum 
dependence plays an important role, LPA may not be a good approximation. This proved 
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to be the case for nonrelativistic Bose-Einstein condensation of cold atoms. The shift of 
the condensation critical temperature has been successfully described [175] within a FRG 
analysis with a more sophisticated approximation scheme for the flow known as the BMW 
approximation [182] , whose main feature is the inclusion of momentum dependence for the 
n— point functions. 

7.3 FRG analysis of pion BEC 

The specific physical problem we are interested in is relativistic Bose-Einstein condensation 
of pions in isospin dense media. In analogy with what was done in the nonrelativistic case 
in cold-atoms physics, we would like to understand the role played by (nonperturbative) 
interactions in the second-order transition between the normal and condensed phases. 

Being a renormalization group based technique, FRG should describe well the physics 
of a second order phase transition. The IR phenomena intrinsically present near criticality 
(due to the massless zero mode that condenses) is exactly the type of nonperturbativeness 
that the flow was constructed to deal with. 

The fact that this method was successfully applied to the nonrelativistic case (where 
extensive and well-controlled experiments are available) suggests that FRG should be an 
adequate approach to this relativistic version. The usage of the same semi-analytic technique 
in both cases will allow us to identify similarities and differences between the nonrelativistic 
and relativistic BEC, so that in principle a lot can be learned about the fundamental features 
of BEC. 

7.3.1 Results within the Local Potential Approximation 

Adopting the pion chiral effective model discussed in Subsection 7.1.1, a FRG analysis using 
LPA was implemented to describe the BEC transition in the temperature versus isospin 
chemical potential plane [183]. The result is shown in Figure 7.2. 

Consistently with the heuristic argument in the beginning of this chapter, BEC conden- 
sation of pions occurs for any nonzero isospin chemical potential in the chiral limit. At the 
physical point, the pion mass is nonzero = 138 MeV, so that condensation takes place 
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Figure 7.2: [extracted from Ref. [183]] T-fii phase diagram within LPA. 

only when /i/ > m.^. Actually, as stated in the paper [183], the second-order phase transition 
does not occur at fij = m^, but slightly below it. Another puzzling result is the transition 
line in the limit of T — )■ 0: it does not appear to be perpendicular to the fij axis as commonly 
happens. 

The purpose of our work is to understand these issues and, in particular, find out whether 
they are unphysical artefacts of the nonperturbative method and approximations used. 

7.4 Going beyond LPA in a toy model 

Let us analyze the same phenomenon of relativistic Bose-Einstein condensation, but in a 
simple field-theory model: a complex scalar field n = (tt^ + i7i'^)/\/2 with real mass {m? > 0) 
and U{1) symmetry at finite density. Our goal is to study the issues encountered in the LPA 
computation of the phase diagram in the T — fij plane in a simpler theory that contains 
most of the features of interest. In what follows, physical similarities of this toy model and 
the chiral model presented above are made explicit, especially in what concerns the phase 
transition due to condensation and the fluctuations that should play an important role in the 
vicinity of criticality. This toy model will allow us to go beyond LPA in a more controlled 
fashion and address the questions concerning the quality of the nonperturbative method and 
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its approximations. 

The euclidean action Se (with the partition function given hj Z = Tr exp(— S'^;)) includ- 
ing the U{1) charge conservation constraint is: 

Se = j d^x {-(9o + i^ii)TT*{dQ - ifii)7r + (Vvr*) ■ (Vvr) + m\*7r + X^Tr*^^} (7.64) 

which can be rewritten as: 

Se = j d^x {d^Tr*d^TT - z/i/ [7r*9o7r - tt^ott*] + (m^ - iJ.j)7r*TV + \^{7r*n f} , (7.65) 

or, in terms of the real {n^/\/2) and imaginary (7r^/-\/2) parts of the field, 

/I 1 

+ i^r] + ^[{^r + i^r?} ■ (7.66) 

Comparing the latter form with the problem of the linear sigma model at finite isospin 
density, one can see that this is a simplified version of it in which we consider only the 
charged pion directions and — A^f^ i— )■ m^. 



7.4.1 Condensates and fluctuation spectrum — classical results 

With the aim of understanding the phase structure of the model and the quasiparticle spec- 
trum at finite density, we can rewrite further the effective action as an expansion around 
(constant and uniform) condensates, parameterized by the real parameters p and 6 as follows, 

n-^i^x) = p cos 9 + 7]^ (x) , (7.67) 
Tr'^ix) = psinO + r]'^{x) , (7.68) 
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where rji are the fluctuation fields around the condensates. One obtains [def.: rj^ = (77^ cos6' + 



7]"^ sin ( 



(3V\ 



2 ^+T^I + 



d x< 



[m 



+ 



2\2l 



4 L 



+ 



4 



(7.69) 



where we have omitted a term of the form oc J d'^x dot], which gives only a surface contribu- 
tion. The second line yields the equilibrium condition for the condensate solution: 



p(m^ - pj) + A>' 







p = or p = — — — 



(7.70) 



Since p G M by construction (and A^ > to guarantee the stability of the potential), the 
nonzero condensate solution is only possible when [p'j — m?) > 0. Notice that in the absence 
of self-interactions (A^ = 0) condensation is only possible at pj = m? (this reproduces 
Kapusta's result shown in the p— axis of Fig. 1 of Ref. [184]) and the theory becomes badly 
defined for chemical potentials above it (since the positivity of the action is violated). 
In Fourier space, the effective action becomes: 



Be (m^ -p]) . A2 



(3V 



-P 



d^x< 



P 



(m^ - pj)p + AV^ 



^ * V2{-Q) 



+ 



(7.71) 



where we have defined the matrix of the inverse propagator (def.: Uc = A p (1 + 2 cos 
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and Us = A^p^(l + 2sm^6')) for the fluctuations around the condensate solution: 



-Ql + (q)'^ + m'^-^i] + u^ sin(2^) - 2i^iiQ 





AV sin(2^) + 22/i/Qo -Ql + (g)^ + - fij + Us 



(7.72) 



The eigenvalues of the inverse propagator are directly obtained: 



2 



= -{jQl + ^Tf^A +(gT + m^ + 2AV(l + ||) • (7.73) 



The dispersion relations for the quasiparticles are then given by the zeros of these eigen- 
values: 



uj'^{g) = (g-)^ + m^ + /if + 2AV±y'AV + 4/i?[(g)2 + m2 + 2AV]- (7.74) 

In particular, the standard result for the dispersion relations in the absence of condensates 
(p = 0) is consistently reproduced: 

limw|(g) = \^{q)^ + ± fiiY . (7.75) 

When fij > m, the equilibrium condition for the true minimum yields = = {fi'j — 
m?)/\^ (cf. Eq. (7.70)) and the dispersion relations in this condensed phase become: 



lim a;i(g) = {q)"^ + ?,p] - ± \ Up] - m'^f + Ap?Mq)'^ + 2p?r - 



= {q)' + 3pj-m'±^J{3pj-m^y + 4p%q)\ (7.76) 

whose associated mass spectrum, obtained in the limit g — )■ 0, is: ml = 2(3pf — m^) and 
ml = 0. The mode b is clearly the massless Goldstone boson in the condensed phase. 

Fig. 7.3 displays the mass spectrum as a function of the chemical potential. Comparison 
with the spectrum of quasiparticles obtained for the chiral pion model. Figure 7.1, shows 
that not only the BEC transition appears in this simplified model as a second-order one, 
but also the spectrum of fluctuations is similar to that of the charged pions. Therefore, 
this toy model contains the relevant features of relativistic BEC, including the spectrum 
of quasiparticles directly involved in it. Moreover, the extra complication regarding chiral 
spontaneous symmetry breaking is absent, allowing for a clean study of BEC physics and 
the role of interactions in it. 
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m 



Figure 7.3: Mass spectrum as a function of the chemical potential fij. Units of mass m. 
7.4.2 FRG analysis within LPA 

Now, in order to study nonperturbative effects, we will implement a FRG flow. As discussed 
in detail in subsection 7.2.1, a standard procedure of deforming the theory yields the FRG 
flow equation, Eq. (7.63), for the k— dependent effective action: 



where Q is the full propagator matrix. 



iiv 

2 



(7.77) 



(7.78) 



/32 6n^{Q)6TTJ{-Q) 

and Rk is the regulator, which we choose to be of the Litim form [178]: R^iff) 

To proceed further and obtain a closed, explicit form for the flow equation, one has to 
assume a form for the k— dependent effective action F^, from which one can write both the 
left hand side and the full propagator itself. In practice, one must select an approximation 
or truncation scheme to work with. 

First, it is interesting to analyze the toy model within the Local Potential Approximation. 
This serves two main purposes: (z) to compare further the findings within the toy model 
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with the ones from the chiral model for pion condensation [183] (cf. subsection 7.3.1) and 
(ii) to provide a reference for comparison with the on-going investigations beyond LPA [174]. 

Let us therefore consider the LPA ansatz for the flowing effective action, i.e. in the 
form of the classical euclidean action as described above, but with k— dependent parameters 
instead of m^, A^, /i/ as follows: 

= J d'^x^^dyoy + ]^d,Tx^d,Tx^ + /i, [n'doTi^ - TT^aovri] + K («) , (7.79) 

where we have defined a = (vr^)^ + (tt^)^ and 

K(a) = V^+ ^'^^'^'K + ^a'. (7.80) 

Notice that the value of effective potential at zero field (V^) will be in general a nontrivial 
function of temperature T and chemical potential /i, being directly related to the pressure 
of the system (and its derivative with respect to fi, to the density). 
Partially transforming to Fourier space, we have: 

( .... o,.. ^ ( ^'(-Q) 



r.[vri = ^^(vri(g) n\Q)))C{Q){ } ^' + / rf^x K («) , (7.81) 
•Jq ^ ^ ^ \ n\-Q) I -J 



where 



2z;u,go -Ql + (g)' 

The full inverse propagator associated with this effective action is then: 



13^ d7r^{Q)67rb{-Q) 
= m) + / ^^^V:{a)n%x)n\x) + 2r^(«)5"^| . (7.83) 

In the case of constant fields (consistently with the LPA ansatz), it simplifies to: 
-Ql + (g-)2 + 4K'(«)(7ri)2 + 2V^{a) 4K'(«)vrV2 - 2zf,^Q, 

4K'(«)vr^vr2 + -Ql + (q)' + 4K'(«)(vr^)^ + 2K(«) 



G-\Q) 



(7.84) 
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while the effective action reduces to the potential (note that one of the chemical potential 
terms involves derivatives and goes to zero in this approximation): 



VJa) 



Using the Litim regulator RK{(f) = {n'^ — (f)9{H?' — (f), we obtain: 



(7.85) 



(7.86) 



(7.87) 



so that the flow equation becomes: 



2k9{k' - f) Tr <^ g-' go, r ^ r + («:' - qio{^^" - r) 



-, — tttK — > ir 
27r 3 3 ^ 

^ ' Qo 



-1 



|7.88) 
(7.89) 

(7.90) 



and we are left with the Matsubara sum of the trace of the full propagator. The trace, in 
turn, can be computed in a straightforward way from the expression for Eq. (7.84). 
The result is [defs.: v^ = k^ + 2aV;! + 2^^ and = y^jJ^Hn^ + 2V^) + {fil + aV^')^]: 

-Ql + 



[-Ql + v^]^-Aa^[V:]^-AfilQl 

-Ql + Vk 

[Ql - {v^ + 2/i2 + 2w^)] [Ql - {v^ + 2/i2 - 2w^)] 



(7.91) 



Ql + E^, 



Ql + E- 



2 ' 



(7.92) 



where we have defined: 



E 



± 



K 



+ 2aV': + 2V', + 2/i^ ± 2^iil{K^ + 2V>) + {iil + aV;!Y 



(7.93) 
(7.94) 
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Explicitly, the flow equation for V^, Eq. (7.90) becomes: 
'9k K 



I //,2 



(7.95) 



It is therefore clear that the Matsubara sums to be evaluated in the flow equations for 

(n) 

and its derivatives Vk will all be of the form 



Qo 



E 

Qo 



with 



and: 



Qo ^ ^' Qo 



where ni,{uj) = [1 — exp(/3a;)] ^ is the Bose-Einstein distribution. 



The final closed form for the flow of the effective potential \4 in LPA is: 



12^ 



^E'i 



1 



1 



l + 2nh{ \lEl^ 



with: 



El 



K 



The potential in our ansatz is: 



VJa) 



(7.96) 



^[n,(v^)-n,(-v^)] =-^[l + 2n,(v^)] , (7.97) 



(7.9^ 



(7.99) 



(7.100) 
(7.101) 
(7.102) 

(7.103) 
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so that its derivatives give: 

V^ia) = -^Y^ + fa, (7.104) 

= A^ (7.105) 

In this case, the parameters and above become: 

= V/i2(K2 + ml -^l + Xla) + (/.2 + aXir (7.106) 

= v//i^(«:2 + m2+3aA2) + a2At, (7.107) 

^2 = ^2 + ^2+ SaXl + fil± 2 V/i2(fi:2 + m2 + SaA^) + ^2^4 . (7.108) 



Some simple hmits are verified: 

• Free (A^ = 0) or noncondensed (a = 0) theory: 



Ei 



(^V^k2 + m2 ± /i^ 



(7.109) 



(7.110) 



The free fiow equation reduces to: 
.4 1 



127r2 ,/fi:2 + ml 



2 + 2nf,(VK2 + m2+/i,) +2725(^^2 + m2-/i,) J7.III: 



It is interesting to note that in general the field-independent part of the effective 
potential will satisfy a fiow equation of the form of the free fiow equation above. This 
is demonstrated straightforwardly by setting a = on both sides of Eq. (7.99): 
A I 



K, 



127r2 + ^2 



2 + 2 (Vk2 + ml + + 2 n, (^y/n^ + ml - ](7.112) 



Vanishing chemical potential {fi^ = 0): 



-■2 ,.2 , _2 I /Q _i_ o\ „ \2 



E^ = + mi + {3±2)aXi. 



(7.113) 
(7.114) 
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The flow at = yields: 



127r2 



2 rib 



l + 2nA \ E_ 



. (7.115) 





= 0) 


= 9.v«, 




= 0) 






= 0) 


= 



The full flow equation for the effective potential V^, Eq. (7.99), can be translated into a set 
of coupled equations for the parameters V^, "m^, /i^ and of our ansatz (def.: = m^— yU^): 

(7.116) 
(7.117) 

where the left-hand side is then written in terms of (derivatives with respect to a of) the 
right-hand side of the full flow equation, Eq. (7.99) evaluated at a = 0. 

Notice that with this ansatz for the effective potential and the assumption of constant 
fields we cannot fix both the flow of the chemical potential and that of the mass indepen- 
dently. Pragmatically, we assume a fixed /i^ = /i. 

Since Vk does not appear at all in the right-hand side of Eq. (7.99), its flow equation 
(explicitly given in Eq. (7.112)) can be solved independently, after the set of two coupled 
equations (7.117) is solved. 



Numerical results at zero temperature 

In what follows we investigate numerically the solutions of the set of coupled flow equa- 
tions in Eq. (7.117), with their left-hand side being given by the a— derivatives of Eq. (7.99). 

We use the following procedure to fix the parameters of the model (notice that, as dis- 
cussed above, /i^ is fixed at /i, which will play the role of the control parameter to investigate 
the phase structure of the theory at zero temperature): 

1. The initial conditions of the flow will be given at a scale k = A: fh\ = — /i^ and 
\\ = g. We choose: 

A = 1000 a.u. 

= -1000 a.u. 
9 = 2 (7.118) 
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2. Consider the flow ai jj = 0. Using the initial conditions above, the physical values, 
obtained in the end of the flow (at k = 0), are (cf. full k dependence in Figs. 7.4): 

ml{fi = 0) = 38487.1 a.u. (7.119) 
A2(/i = 0) = 1.47035. (7.120) 
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II 
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0- 




200 400 600 800 1000 
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Figure 7.4: Flow in the vacuum (T = yU = 0): 




200 400 600 800 1000 
K [a.u.] 

versus k (left) and versus k (right). 



We use rholfi = 0) = M to set the unit of energy in the problem. As shown in the 
analysis of the spectrum of excitations in the classical approximation above, we expect 
condensation to manifest at /i = mphys(/i = 0), i.e. in our case: /i = rfiQ^fi = 0) = M. 

It is interesting to note the saturation of the flow when n < M. This is a generic 
feature of the FRG flow that reflects the fact the below n = M the IR modes of the 
theory are naturally suppressed and the regulator i?^ plays no role in that anymore. 

With the scale set, we turn on the chemical potential in the flow equations and look for 
the critical value /icrit at which rho{fi = ficnt) = 0. The vanishing effective mass characterizes 
the second-order phase transition corresponding to Bose-Einstein condensation. 

The flow for different values of chemical potential is shown in Figs. 7.5 and 7.6. The 
effective mass rfiK^) decreases as expected when one increases the chemical potential fi and 
goes to zero when /i — )■ M. We find: 

^2^2 ^ 0.996 M^) 

^ ^ = 0.0000693186. (7.121) 

However, above ^ 0.996 M^, numerical instabilities appear, signaling the criticality before 
the expected value /i^^j'^ = M. Notice that this feature was already encountered in the LPA 
FRG analysis within the chiral model [183] (cf. Fig. 7.2). 
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>^ 
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KfM 



Figure 7.5: Flow for different /i's (T = 0): ml versus n for fi'^/M'^ = {0.026 (top), 0.260, 
0.650, 0.996 (bottom, dashed, critical)}. 

Actually, this can be understood as the effect of the running of the physical mass (m^ = 
ml + yU^) at zero temperature and finite /i, as illustrated by Fig. 7.7. 
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Figure 7.6: Flow for different /I's (T = 0): versus k for /x^/M^ 
0.650, 0.996 (bottom, daslied, critical)}. 
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Figure 7.7: Physical mass at finite fi, obtained at the end of the flow (k — )> 0). 



This running mass at zero-temperature and finite isospin chemical potential is also re- 
sponsible for generating a small but nonzero net density in the normal phase. This can be 
verified as follows. The LPA result for the pressure and the density of the system can be 
obtained from the flow of Vk(T, /i) = —Vk{T,^), given by Eq. (7.112). At zero temperature 
and low /i (below the condensation point), we have: 



9«V40,/x) 



A 



1 



67r2^K2 + m2(0,/i) 



(7.122) 
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In this case, we can fix the initial condition for the flow by setting the pressure to zero in 
the vacuum (T = /i = 0): 



A 

dn - 



(7.123) 



/o 67r2 ^«:'2 + ^2^(0,0) ' 

and assuming that, at the scale k = A, fluctuations which are T—, yU— dependent are totally 
suppressed and the potential Vk=a is independent of T and /i. Using this result, the solution 
of the flow equation (7.122) can be written as: 

r-A ,JA 1 



67r2 ^«'2 + ^2^(o,/i) 



(7.124) 



which can also be done with the zero-point energy term in the thermal case. 

The density is then given by the derivative with respect to rii^iT^jj) = '^'^"gj'^^ = 
— At zero temperature this is directly related to the /i— derivative of the mass 



dpi 



(which is clearly nonzero in Fig. 7.7): 



djji 



A 

dn' - 



127r2 (^K'2 + m^,(0,/i))3 
The numerical results are given in Fig. 7.8. 



(7.125) 
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Figure 7.8: Density versus chemical potential /i [left; n/K = O(top), 0.1, 0.5, 1 (bottom)] 
and flow variable k [right; /x/M = 0.88(top), 0.51, 0.16, 0.0005 (bottom)]. 
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Numerical results at finite temperature 

Using the same procedure and initial conditions, we analyze temperature effects. Fig. 7.9 
shows the ffow evolution of the parameters of the effective potential at finite temperature. 




Figure 7.10: T x /i phase diagram for BEC. The darker region is the condensed phase. 



A previous result for this line^ was obtained [184] within mean- field approximation 

^In Ref. [185], the author obtains Tciit{n), where n is the density, within dimensional reduction (with 
lattice input), so that the comparison is less direct. 
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(MFA), 

T^i = if^'-m^). (7.126) 

As discussed for the results at zero-temperature, the shift of the critical line in LPA 
as compared to the MFA is due to the effective running of the mass (i.e. the quadratic 
coefficient of the potential) in a medium. The hotter the environment the larger is this 
effect, as depicted in Fig. 7.11. 




Figure 7.11: Mass spectrum as a function of the chemical potential /i for different temper- 
atures: T/M = {0.30584(bottom),0.611679, 1.01947, 1.42725,2.03893 (top)}. The thick red 
line is the "critical" curve 171^=0 = /i, while the dashed one is the reference 771^=0 = M. 

The comparison between MFA (with = 2 and m = M) and our LPA results is shown in 
Fig. 7.12. This plot explicitates the fact that the critical line within the LPA FRG analysis 
is for most values of fij, but not all, below the one within MFA. Physically, however, one 
expects that, at fixed /i/, the treatment of fluctuations - as implemented by FRG - should 
in principle always enlarge the disordered phase, i.e. the noncondensed one. This might be 
yet another peculiar feature of the LPA FRG result. 

The second-order critical line obtained within LPA for this toy model is qualitatively very 
similar to the one in Fig. 7.2, computed within the pion chiral model, even though we have 
used a simplified LPA ansatz (neglecting O([0^]^) terms). Once again, this simple field-theory 
appears to be a good toy model for the description of relativistic Bose- Einstein condensation. 
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Moreover, the issues encountered in the analysis of Ref. [183] (namely, fif^^{T = 0) < m^r 
and the transition line being nonperpendicular to the //-axis at zero temperature) are also 
present in this toy model computation and may be scrutinized in this cleaner framework. 

1.0 

0.8 
^ 0.6 
^ 0.4 

0.2 

0.0 

0.90 0.95 1.00 1.05 1.10 1.15 1.20 1.25 

Figure 7.12: Critical lines obtained within LPA (red, solid line) and MFA (black, dashed). 

7.4.3 LPA': introducing the Z— factor and momentum dependence 

To investigate whether the issues encountered in the result for the phase diagram of rela- 
tivistic Bose- Einstein condensation in isospin-dense systems within LPA (Figures 7.2 and 
7.12) are indeed artefacts of this approximation, we would like to implement a FRG flow 
going beyond the zeroth order in the derivative expansion. 

In order to incorporate the effect from momentum-dependent interactions in our com- 
putation of nonperturbative corrections, we go one step further in the derivative expansion 
for the effective action, including the k— dependent wavefunction renormalization in our 
ansatz. This improved truncation scheme is the so-called LPA'. 

Our LPA' ansatz for the effective action is therefore (recall that a = (vr^)^ + (tt^)^): 

= /" d^xl ^d^Ti'd^Ti' + ^d,Tx^d,Ti^ + /i, [Tr^Sovr^ - tt^^oTtI] + K (a) \ (7.127) 
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Following Ref. [186] (section 5, in particular), one can obtain the flow equation for the 
wavefunction renormalization by expanding the fleld content in the flow equation for the 
effective action in the following way: 



7r(x) = ip + A(x) 



(7.128) 



with If being a constant background fleld (that corresponds to the full Q = contribution 
to the fleld) and A(x), a small inhomogeneous fluctuation, i.e. satisfying A (a;) <^ ip and 
A{Q = 0) = 0. One can then show that the wave function renormalization is given by: 



Zf- = — lim 



/32QCo9g2 5Ai(-g)5Ai(g) A=o ■ ^^'^'^^^ 

Indeed, if we make the substitution (7.128) in our LPA' ansatz for the effective action, Eq. 
(7.127), we obtain: 



(7.130) 



where we have deflned = {^^f + {^^f and 6 = 2ip^A^ + 2ip^A^ + (A^)'^ + {A^ = 0(A). 
Expanding up to quadratic order in the inhomogeneous fluctuation A, we have: 



d^xi -A' 



A^' + 



+ U^if'-^do - fi.^'-^'do + 2^XK))^' + K(«^) + 0{A 



(7.131) 



or, in momentum space, 

r. = Y^PV)^ A\Q) [6^^ [Z,[-Ql + (f] + 2K'(«^) + 4a^K'(«^)) + 

^ Q 

+6''6''fi^{-zQo)-S''6^'fi,i-zQo)]A\-Q) + 
+PVV^{a^) + 0{A'), 

where we have used the standard Fourier expansion for the fluctuation fleld: 

A(x) = y|(W^e^^-A(Q), 



(7.132) 



(7.133) 
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and A((5 = 0) = in the last line. The functional derivatives of with respect to the 
inhomogeneous fluctuation field A then give: 



5Ab{-Q)6A-{Q) 



6'^'[Z,[-Ql + (f] + 2V;'(a^) + 4«^K"(«^)) + 
+(5'^i5''V.(-^Qo) - 5'^'5'V«(-^Qo)l + 0(A) , (7.134) 



so that (notice that -Ql + (f = Q"^) 



5Ai(-Q)5Ai(g) 



(7.135) 



A=0 



and the expression (7.129) for the k— dependent wavefunction renormalization is verified. 

In our case, we can similarly obtain the k— dependent chemical potential (cf. Eq.(7.134) 
for a = 2 and 6 = 1): 

d 5^T^ 



lim 



(7.136) 



A=0 



P^QoM^odQo 6A^{-Q)6A^{Q) 

whose running is not accessible within LPA, as discussed above. 

The results in Eqs. (7.129) and (7.136) define special projectors^ that, if applied to the 
flow equation for the effective action (Eq. (7.77)), yield flow equations for the k— dependent 
wavefunction renormalization and chemical potential /i^, respectively. 

The derivation and solution of these projected equations for the LPA' ansatz is work in 
progress [174]. 



7.5 Remarks and perspectives 

In this chapter a nonperturbative description of the phase diagram of Strong Interactions in 
the temperature versus isospin chemical potential plane was presented. We have discussed, 
using effective models, the physics concerning the formation of a new state of matter: the 
relativistic Bose-Einstein condensate of charged pions. Nonperturbative aspects were taken 
into account via an FRG flow. 

^We followed the prescription by Wetterich in Ref. [186], but it seems to be completely equivalent not to 
expand the field. The same result is obtained for <^ i— > 0, except that we would not be allowed in general to 
assume A(Q = 0) = and would have to evaluate tt = (p in the end, instead of A = 0. 
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Throughout the chapter it becomes clearer how one may render, through sensible approx- 
imations, the FRG formalism a powerful tool to address in a nonperturbative fashion the 
phase structure of in-medium field theories. The adequacy of approximations and trunca- 
tions is, however, a subtle system-dependent issue. Our formal aim in this on-going work is 
to scrutinize the results provided by FRG and its approximations within a sufficiently simple 
theory containing a physically motivating phenomenon (BEG), gaining understanding of how 
this nonperturbative flow implements nontrivial contributions and what are the limitations 
of the approximations used. 

In the near future, we expect to have results beyond the leading order in the derivative 
expansion [174] that should eventually answer whether the puzzling features found in the 
description of the T — fij phase-diagram are physical or rather artefacts of the LPA ansatz. 



T 



Chapter 8 




Conclusions 



The study of the phase structure of Strong Interactions is an intricate enterprise. Being 
essentially an open problem, it generates an intense demand for theoretical development, 
since features of the different phases and the transitions between them influence the physics 
of observables in relevant contexts as, for instance, the evolution of the early universe, 
compact star structure and supernovae explosions. Moreover, the critical region of the 
QCD phase diagram is currently within the reach of ultra-relativistic heavy-ion collision 
experiments, especially via the on-going program Beam Energy Scan at RHIC-BNL and the 
facility FAIR at GSI which were designed to probe media with intermediate values of baryon 
chemical potential. This intense demand comes up against technical difficulties, so that it 
is inescapable to resort to various alternative/complementary methods, which range from 
nonperturbative, effective and out-of-equilibrium approaches to the study of cousin theories 
(with different masses and numbers of colors and flavors) and various macroscopic control 
parameters (temperature, chemical potentials, magnetic field, system size, etc). In this way 
the genesis of the different phase diagrams of Strong Interactions which are the subject of 
this thesis takes place. 
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The complexity of the problem is reflected in the plurality of projects and formalisms 
involved. The frameworks adopted included from universal statistical mechanics and simple- 
minded phenomenological models to intricate in-medium gauge theory and sophisticated 
nonperturbative methods. These were the tools used in different contexts to approach ex- 
perimental and observational conditions or the domain of applicability of lattice QCD, with 
the aim of gathering complementary information about the physics of Strong Interactions in 
a medium. 

The projects presented ranged from model-building with lattice input and HIC data 
analysis to the theoretical study of nonperturbative methods. The guideline underneath these 
various directions is unique, however: posing a well-defined question about the behavior of 
Strongly interacting matter, identifying the relevant ingredients in this context and adopting 
a sufficiently simple framework that still has the capability of accommodating most of these 
features, keeping close connection whenever possible to experiments, observations and lattice 
QCD. Frequently drastic simplifications are needed, but in answering these, sometimes simple 
and isolated, questions, the hope is to advance in the puzzle of constructing the phase 
structure of Strong Interactions as it is realized in Nature. 

Chapter 2 illustrates how considering extra axes (mass and isospin chemical potential, in 
this case) in the phase diagram of Strong Interactions enhances sensibly the potential of con- 
straining effective models. We have provided a unified description of the functional behavior 
of the deconfining critical temperature Tc{mT^,fij) in reasonable agreement with available 
lattice QCD data and also predictions for the cold and dense critical line /xg'*(m,r, Z^/)- Even 
though our model is assumedly crude (it fails completely, for instance, in predicting the 
crossover nature of the transition), the task of describing Tc{mT^, fij), as far as we know, is 
currently not fulfilled by the most popular effective models for the QCD phase transitions. 
If not a striking coincidence, the results suggest that (i) the location of the critical region 
is reasonably well determined by the consistent description of high and low energy extreme 
domains and (ii) our crude model probably encodes the relevant ingredients in both high 
and low energy sectors. Moreover, since chiral models fail completely to describe Tc(m,r), 
it seems that the critical temperature is an observable driven dominantly by confinement 
rather than chiral symmetry. Further analyses in this direction might help understanding 
the interplay between chiral and confinement features and hopefully building better models. 
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Approaching much closer the physics of HIC experiments, Chapter 3 presents an analysis 
of the pseudocritical phase diagram of Strong Interactions in a finite system. Finite-size 
effects on the experimental CEP search are extensively analyzed, both as a drawback of the 
known signatures and as the possibility of a new alternative signature based on FSS. First, 
it is shown that the sizes of the systems created at RHIC are comparable with the scale of 
the QCD critical region, causing the pseudocritical behavior to suffer significant shifts in the 
T — fiB plane. Estimates for these shifts within the LSMq were obtained and consequences for 
the signatures of the second-order CEP were discussed, mainly representing a caveat in the 
search for nonmonotonic signals. The fact that HIC data can be organized as coming from an 
ensemble of systems of different sizes, however, brings a bright side of having significant finite- 
size corrections to the critical correlations: FSS, a universal statistical mechanics technique, 
can be used as a tool for the CEP search at HICs. Of course, a lot of caveats should be kept 
in mind, as discussed in the beginning of the chapter, due to the very complicated nature 
of the physical process occurring in a HIC. Nevertheless, we have pragmatically shown, in 
collaboration with an experimentalist from the STAR collaboration at RHIC, that a FSS 
analysis of data is feasible, has the capability of furnishing information about the location 
of the CEP and provides falsifiable predictions. 

In Chapter 4 we turn our attention to the cold and baryon-dense region of the phase 
diagram. In this regime, some observable astrophysical phenomena are connected with the 
microscopic dynamics of Strong Interactions, via the possibility of nucleating deconfined 
bubbles in extremely dense systems such as compact stars and core-collapse supernovae. 
The key physical quantity coming from theory as an input to the modeling of astrophysical 
processes is the surface tension of the QCD phase transitions in the cold and dense regime. 
This is, however, not known via first-principle methods, due to its genuine nonperturbative 
nature and the Sign Problem on the lattice. We have provided estimates within the LSMq 
for the nucleation parameters of the cold and dense chiral phase transition, showing that the 
values obtained are small enough to allow for the proposed observable signals of quark matter 
formation in ultra-compact astrophysical objects. Finally, it is interesting to notice that 
astronomical observations might be entering a precision era in which it will start constraining 
the physics of Strong Interactions in the cold and dense domain, so that further theoretical 
predictions in this direction are called for. 
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DjTiamics of quark bubble nucleation in supernova explosions is an example among sev- 
eral of the importance of nonequilibrium approaches in connecting theory and observations. 
In this work, we have concentrated on estimating time and energy scales associated with 
the dynamics of phase conversion by assuming the stationary scenario of homogeneous nu- 
cleation. Applying full out-of-equilibrium quantum field theory in a Strongly interacting 
medium is a crucial but extremely hard task. One alternative direction is to try to describe 
the approximate Langevin dynamics of the order parameters of the QCD phase transitions, 
in which memory effects should play an important role, especially since we are dealing with 
a relativistic system (as occurs in relativistic hydrodynamics). 

In Chapter 5 a hot QCD medium exposed to a constant and uniform external magnetic 
field is considered in an on-going project. The temperature versus external magnetic field 
plane of the QCD phase diagram combines three desirable features: experimental relevance, 
capability of being investigated on the lattice and a rich phenomenology. The extra mag- 
netic axis opens, therefore, a new channel of comparison between theory and numerical 
approaches. Our aim in this on-going project is to make proper use of this opportunity 
by providing predictions within the fundamental gauge theory for the thermodynamics of a 
hot QCD environment immersed in an intense magnetic field. Questions like "what are the 
adequate quasiparticles?" and "how does the non-Abelian interaction affects the pressure 
in this extreme regime?" are investigated analytically within this well-defined context. The 
background classical field is treated nonperturbatively as a magnetic dressing of the quarks, 
resulting in an effective dimensional reduction of the dynamics for very large fields. We have 
constructed the building blocks of the framework and discussed its subtleties at very large 
B. Preliminary results for the two- loop pressure of a gas of gluons and magnetically dressed 
quarks are presented. We are looking forward to comparing our findings to recent lattice 
QCD simulations and to gaining, in this way, further information on the quasiparticles in 
this regime of the QCD phase diagram, but also on the behavior of thermal perturbation 
theory in this domain. 

Towards having a robust description of the phase diagrams of Strong Interactions within 
effective field theories, it is crucial to go beyond the mean-field approximation and test for 
consistency. As discussed in Chapters 6 and 7, this is an involved quest, especially because 
effective descriptions of a QCD medium usually require several ingredients, such as finite 
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temperature and various chemical potentials, finite masses and condensates. The adequate 
nonperturbative treatment will ultimately be dependent on the specific system treated and 
its main ingredients and shall be in practice considerably restricted by technicalities. In this 
thesis, we have used two nonperturbative frameworks to address completely different physical 
systems: the Optimized Perturbation Theory was applied to describe the thermodynamics 
of a general massive Yukawa theory in the presence of condensates in Chapter 6 and the 
Functional Renormalization Group is applied, in the on-going project of Chapter 7, to the 
relativistic condensation of pions in hot isospin-dense media. 

We have shown that a two-loop implementation of OPT consistently reproduces the 
mean-field result when only Hartree-like (but not Fock-like) contributions are included and 
further improves over MFA when the class of exchange diagrams are resummed. Nevertheless, 
this is an intricate framework, especially due to the complications of renormalization at higher 
loops, and, as discussed in the final remarks of Chapter 6, seems to encounter difficulties to 
describe a stable effective potential in the physical case of a chiral model. 

Finally, in Chapter 7 the formalism of the functional renormalization group was presented 
and preliminary results on the T — fij phase diagram of a scalar toy model for relativistic 
Bose-Einstein condensation are shown. We believe that within this simplified toy model, 
we may apprehend the physics driving the phase transition to the condensed phase and, in 
particular, the role played by nonperturbative interactions. Our expectation is, as suggested 
by the analogous nonrelativistic problem, that the LPA (i.e. the zeroth order of a derivative 
expansion) will prove not to be a good approximation to describe BEC physics, in which 
momentum-dependent interactions seem to be important. 

Of course, the problems addressed in this thesis represent only a restricted exploration 
of a subset of the various phase diagrams of Strong Interactions. There is still much to be 
explored and learned in these planes and others. Not to mention the, much less investigated, 
phases of in-medium Electroweak Interactions and Beyond the Standard Model theories. In 
any of these rich physical scenarios, we believe that the development may be accomplished 
through an optimized utilization of complementary frameworks, as done throughout this 
thesis. 



Appendix A 



Effective theory for the 
chiral phase transition 



The Linear Sigma Model with constituent quarks (LSMq), also known as the quark- meson 
model, is very suitable for the study of the chiral transition [187]. As argued in Ref. [188], 
QCD with two flavors of massless quarks belongs to the same universality class as the 0(4) 
Linear Sigma Model (LSM), exhibiting the same qualitative behavior at criticality. Besides, 
the LSM is renormalizable [189] and reproduces correctly the phenomenology of QCD at low 
energies, such as the spontaneous (and small explicit) breaking of chiral symmetry, meson 
masses, etc. In spite of the fact that an effective theory does not require renormalizability 
to be well posed, this attribute is highly desirable if one is interested in investigating the 
behavior of physical quantities as the energy scale is modified, which can be accomplished 
via renormalization group methods. 

Since its proposal [79], the LSM has been investigated in different contexts, from the low- 
energy nuclear theory of nucleon-meson interactions to ultra-relativistic high-energy heavy- 
ion collisions, and also in different varieties, e.g. including explicitly constituent quark 
degrees of freedom or not. The thermodynamical aspects of the model were first considered in 
the very early days of finite-temperature field theory [190] , and systematic approximations for 
the study of the chiral transition started to be implemented soon afterwards [191], producing 
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an extensive literature. 

When used to mimic and study the chiral transition of QCD, the emphasis was generally 
on thermal effects [81, 82, 158, 191, 192, 193, 194, 195] rather than considering a cold and 
dense scenario, although chiral symmetry restoration at high densities was predicted quite 
early [196]. This choice was, of course, well justified by the stimulating experimental results 
coming from high-energy heavy- ion collisions [25], and by the possibility to compare to 
numerical output from lattice QCD [55]. Usually the gas of quarks is treated as a thermal 
bath in which the long-wavelength modes of the chiral field evolve, the latter playing the 
role of an order parameter in a Landau- Ginzburg approach. The standard procedure is 
then integrating over the fermionic degrees of freedom, using a classical approximation for 
the chiral field, to obtain a formal expression for the thermodynamic potential from which 
one can compute all the physical quantities of interest. In this case, the sigma field is 
approximated by the condensate, and the functional integral over sigma fluctuations is not 
performed. The fermionic determinant is usually calculated to one-loop order assuming a 
homogeneous and static background field [149]. 

In a theory with spontaneous symmetry breaking, the presence of a condensate will 
modify the masses. In particular, in the case of the LSMq the masses of quarks and mesons 
will incorporate corrections that are functions of the chiral condensate, which is a medium- 
dependent quantity. Therefore, contributions from vacuum diagrams can not be subtracted 
away as trivial zero-point energies since they contain medium- dependent pieces via effective 
masses. So, although the presence of the medium brings no new ultraviolet divergence, 
in principle one has to incorporate carefully finite vacuum contributions that survive the 
renormalization procedure. These contributions have been very often discarded in studies 
of the LSMq, but were shown to play an important role by the authors of Ref. [173] who 
incorporate scale effects phenomenologically. Vacuum contributions were considered at finite 
density in the perturbative massive Yukawa model with analytic exact results up to two loops 
in Ref. [47, 48] and, more specifically, in optimized perturbation theory at finite temperature 
and chemical potential in Ref. [152], also comparing to mean-field theory. More recently, 
this issue was discussed in a comparison with the Nambu-Jona-Lasinio model [197], in the 
Polyakov-LSMq model in the presence of an external magnetic field [131], and in the quark- 
meson model, with special attention to the chiral limit [198]. 
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To describe the chiral phase structure of strong interactions at finite density and vanishing 
temperature we adopt the LSMq, defined by the following lagrangian 

C = ipf[i'y''d^-mf - g{a + i'j5T-if)]ipf 

+ ^{d^ad^a + d^n-d^'n)-Uia,n) , (A.l) 

where 

U{a,n) = ^{a^ + 7i^-v^f-ha (A.2) 

is the self-interaction potential for the mesons, exhibiting both spontaneous and explicit 
breaking of chiral symmetry. The Nf = 2 massive fermion fields ipf represent the up and 
down constituent-quark fields tpf = {u,d). For simplicity, we attribute the same mass, 
TTif = rriq, to both quarks. The scalar field a plays the role of an approximate order parameter 
for the chiral transition, being an exact order parameter for massless quarks and pions. The 
latter are represented by the pseudoscalar field vf = {711,112,113). It is customary to group 
together these meson fields into a 0(4) chiral field </> = (cr, tt). For simplicity, we discard 
the pion dynamics from the outset, knowing that they do not affect appreciably the phase 
conversion process [82], and focus our discussion in the quark-sigma sector of the theory. 
Nevertheless, pion vacuum properties will be needed to fix the parameters of the lagrangian 
later. As will be discussed below, the parameters of the lagrangian are chosen such that 
the effective model reproduces correctly the phenomenology of QCD at low energies and in 
the vacuum, such as the spontaneous (and small explicit) breaking of chiral symmetry and 
experimentally measured meson masses. 

Since chiral symmetry is spontaneously broken in the vacuum, it is convenient to expand 
the a field around the condensate, writing a{x, r) = (cr) + ^(x, r), where r is the imaginary 
time in the Matsubara finite-temperature formalism. Although we assume the system to be 
at temperature T = 0, it is convenient to compute the loop expansion at finite temperature, 
taking the limit T — at the end. For a dense system, the chiral condensate, (cr), will be 
a function of the quark chemical potential, fi. Given the shift above, the fiuctuation field ^ 
is such that {^) = and ^{p = 0,u = 0) = 0. From the phenomenology, one expects that 
((T)(/i — !■ 00) ^ (being equal in the case of massless quarks). 
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Keeping terms up to 0(.^^), one obtains the following effective lagrangian 



(A.3) 



where we have defined the /^-dependent effective masses 



(A.4) 



Linear terms in ^ can be dropped in the action because of the condition $,{p = 0, = 0) = 0. 

In the effective lagrangian, the medium-dependent masses will entangle vacuum and 
medium contributions in the loop expansion, rendering the renormalization process more 
subtle. This is a consequence of the presence of a nonzero //-dependent condensate in the 
broken phase. Of course, the ultraviolet divergences are the ones of the theory in the 
vacuum, and renormalization is implemented in the standard fashion by adding medium- 
independent counterterms to the original lagrangian (A.l). In our case, though, it is more 
convenient to introduce counterterms in the effective theory, so that there are contribu- 
tions from pure-vacuum, pure-medium and vacuum-medium pieces in the cancellation of 
ultraviolet divergences. Renormalization is then implemented by using standard methods of 
finite-temperature field theory [149]. 

A.l Cold and dense 1-loop effective potential and vac- 
uum corrections 

The effective potential can be computed exactly and in closed form following the steps 
detailed in Ref. [47] (cf. also [48]), a procedure that can be even generalized to optimized 
perturbation theory [152], as presented in Chapter 6. Keeping only contributions to one loop 
order, the effective potential is given by 



where f2^™(o") is the fully- renormalized thermodynamic potential for the fluctuation effective 
theory. The latter corresponds to a Yukawa theory for massive fermions and a massive scalar, 
with masses given by (A.4) as functions of a. The thermodynamics of this theory was fully 



(A.5) 
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solved analytically to two loops in the cold and dense regime in Ref. [47], where all details 
can be found. In the MS scheme, can be written as a sum of a medium contribution 



247r2 



fipf - Ml log 



and a vacuum contribution 



647r2 



647r2 



log 



Ml 



(A.6) 



(A.7) 



where u = fipj — log (^^-jp-^, P/ is the Fermi momentum and A is the MS scale. The 
latter can also be fixed by vacuum properties, as will be discussed in subsection A. 3. 



A. 2 Thermal effects 

We can also incorporate thermal effects in the calculation of the effective potential within 
the LSMq. The inclusion of the temperature dependence allows for testing the validity of the 
cold (T = 0) approximation at low temperatures and to investigate if the thermal corrections 
to nucleation parameters can play an important role in the phase conversion process. 

At one loop, the well-known temperature- and density- dependent correction to the ther- 
modynamic potential is that of an ideal gas of massive sigma particles and constituent quarks 
{cul = k^ + Ml and El = + M^): 

(1) - T I '^^^ W h - e-^^/^1 - 
^Vcd.Th ~ ^ J (27r)3 ^ ^ ^ 

-2TNfN, I ^ {log [1 + e-(^'-^)/^] + log [l + e-(^'+^)/^] } . (A.8) 

A. 3 Parameter fixing 

As stated above, the LSMq is adopted as an effective model for QCD at low energies, so 
that the parameters g, A, m^, v, h and A are fixed in order to reproduce QCD properties 
either measured in the vacuum or calculated numerically via lattice QCD. Therefore, the 
conditions for fixing the parameters are imposed on the vacuum effective potential. Since 
we aim at comparing results from cases with different vacuum effective potentials (namely 
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U and U + riiic), and even with different parameter sets^, it is useful to make explicit the 
parameter fixing procedure, consistently. 

The conditions for fixing the model parameters are the following: 

• The chiral condensate in the vacuum is the pion decay constant, /,r = 93 MeV, or, in 
terms of the minimum of the vacuum effective potential. 



OT/vac 



= 0; (A.9) 



d{a) 

The partial conservation of the axial current yields 

h = f^ml = (93 MeV) (138 MeV)^ , (A. 10) 

where is the pion mass; 

The current quark masses calculated in lattice QCD (cf. e.g. Ref. [39]) provide: 

rrig = 4.1 MeV , (A.ll) 

which we neglect^, setting rUg = 0, for the sake of comparison with previous model 
calculations; 

Using the constituent quark mass in the vacuum as 1/3 of the nucleon mass {rriN = 
938 MeV), we can fix the Yukawa coupling, g: 

Mg{{a){pf = 0)) = M^^^ ^=1 Qm;v-m,^ =3.32; (A.12) 

The value of the dressed mass of the a field is given by the experimental value of the 
mass of the a meson: 



a2T/vac 



= {Ml'^y ^ (600 MeV)2 ; (A.13) 



^The case with vacuum logarithmic terms has an extra parameter: the MS scale A. 

^The current quark mass also does not bring any extra qualitative feature to the model nor changes 
significantly the quantitative results concerning the chiral phase transition [48] . 
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The quark condensate is fixed by tiie lattice result (including only quarks up and down, 
we have [39]: (V'/V^/)vac = -2 (225 MeV)^), so that 



GT/vac 



drrig 



= {^f^fUc ■ (A.14) 

vac; {a)=fTr 



The conditions (A.10)-(A.12) above fix directly the parameters h, rrig and g, indepen- 
dently of the inclusion of quantum corrections to the vacuum thermodynamic potential. 
On the other hand, the Eqs. (A. 9), (A. 13) and (A.14) are coupled equations to determine 
the parameters v and A (and A, if quantum corrections to the vacuum are considered) and 
depend on the explicit form of the vacuum effective potential. In the case without quan- 
tum corrections, i.e. with the vacuum effective potential being purely the classical potential 
(V™^ = ^7), we find ^ 20 and ^ 7696.8 MeV^ With the addition of the 1-loop vacuum 
term rilac to the vacuum effective potential, the solution of Eqs. (A. 9), (A. 13) and (A.14) 
yields A^ ^ 16.65, ^ 3296.89 MeV^ and A ^ 16.48 MeV. 



A. 4 Influence of interactions: higher-loop corrections 

So far, we have only included the effects of interactions indirectly in the construction of the 
effective model itself, through dressed masses and the presence of a nonzero quark condensate. 
However, the interaction of the sigma meson with the medium constituent quarks could 
in principle alter the predicted dynamics for the chiral transition. The incorporation of 
such corrections in the calculation of the effective potential can be implemented via the 
perturbative technique order by order. 

The first interaction correction in the present case appears at the two-loop order. The 
contribution to the effective potential corresponds then to that of the thermodynamic po- 
tential of an interacting Yukawa theory with massive scalars and massive fermions, which 
was obtained and analyzed in detail in Ref. [152] (cf. Chapter 6), including the vacuum con- 
tribution and nonperturbative effects within the optimized perturbation theory framework. 
Therefore, one has in principle all the machinery to investigate the influence of interactions 
on the phase diagram of the LSMq and the associated process of homogeneous nucleation. 
It should be noted that the full case, including quantum corrections in the vacuum effective 
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potential, requires the (nontrivial) solution of the coupled equations (A. 9), (A. 13) and (A. 14) 
with the 2-loop result plugged in. Being a more technical analysis, it is out of the scope of 
this Appendix, but some preliminary results and discussion can be found in Ref. [48]). 



Appendix B 

Vacuum thermodynamic potential 
and renormalization at two loops 



In this appendix, we address the details involved in the explicit derivation of the vacuum 
contributions to the two-loop thermodynamic potential of the Yukawa theory. In particular, 
we concentrate on the calculation and renormalization of the 1-loop bubble diagrams and 
the exchange diagram in the vacuum. The vacuum and in-medium direct contributions (the 
third diagram in Fig. 6.1), which were not considered in Ref. [47], are left for the next 
appendix. 

After computing the Matsubara sums (cf. Ref. [47] and the next appendix), the pieces 
which are not explicitly dependent on T and/or n correspond to the vacuum contribution in 



The first two terms in Eq. (6.10) can be expressed in terms of the following UV-divergent 
function 



Those divergences are cancelled by a field-independent counterterm in the Lagrangian, com- 
monly known as a vacuum expectation value subtraction or a cosmological constant. Within 
the MS subtraction scheme, the tridimensional momentum integral above is renormalized to 



Eq. (6.10). 




(B.l) 
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[157]: 

A/f4 / A2 M 

(B.2) 



■3 , f 

- + log — 



2 \My_ ■ 



647r2 

The two-loop T, /i- independent exchange contribution to the thermodynamic potential 
Qy, the last term in Eq. (6.10), can also be written in terms of UV-divergent vacuum 
integrals: 

oexc = ^r I. f d'Pd'Q 4K + P-Q) .R^. 

2 7 (27r)8 _ ^2) (p2 _ ^2) [(Q _ p)2 _ ^2j ' ^ ' ) 

corresponding to the vacuum exchange diagram ^, as shown in Fig. B.l. 



^.^s. V 



VAC 



Figure B.l: Vacuum contribution of the exchange term written in terms of the associated 
vacuum diagram. 

The renormalization is then implemented through the usual procedure, with the addition 
of the appropriate diagrams containing counterterms, as represented in Fig. B.2. 



N. I r ^^'^ ^^VAC ^.X.VAC 



"vac 2 ^ T T "V / '* ' f ""vac ' ""vac 



Figure B.2: Renormalized exchange contribution to the vacuum thermodynamic potential. 
The crosses indicate counterterm vertices. 

The counterterm vertices are defined in Fig. B.3 in terms of wavefunction and mass 
counterterms. At this order within the MS subtraction scheme, these counterterm vertices 
cancel exactly the 1-loop vacuum self-energy poles, yielding (d = 4 — e): 

~ 2(47r)2 e ' (Any e ' ^ ^ ^ 

5(2) 2 . (,)__12m^2 
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Figure B.3: Definition of counterterm vertices. 

Tlierefore, tlie second and tliird contributions in Fig. B.2, tlie vacuum bubble diagrams 
witli counterterm insertions, yield respectively: 



fCT 



Np^ ^ f d^P ^ 



(27r)^ 
:-2)z(z/) 4 



m 



bCT 



2 

N, 



(27r)4 V^^ "™^V K^-ml 



wliere we liave defined the following divergent integral 

d^P 1 



Km' 



(27r)4 P2 _ ^2 



whose regularization yields 



1 + loe 



A2 



(B.6) 



(B.7) 



(B.8) 



(B.9) 



with 



aim 



log — 



1 / 



(B.IO) 



Using the results (B.4), (B.5) and (B.9) in (B.7) and (B.6), we arrive at the following 



^Throughout the appendices, whenever we refer to vacuum diagrams, we adopt the Feynman rules from 
Ref. [145], with factors {—Np) associated with fcrmion loops excluded from the diagram definition. 
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final regularized expressions for the second and third terms in Fig. B.2, respectively: 



fCT 



-12 7 TT 



bCT 



2 ^ (47r)4 



e 
mi 



^ ^log('4)+2«(m2) + 0(6)'>,(B.ll) 



2 2 



e e 



log 



+ 2 a(ml) + 



-0(e) 



(B.12) 



In Fig. B.2, the contribution still to be calculated explicitly in a regularized form is Q 
Using the identity 



exc 
vac' 



m^ + P-Q = 2m^-^ml-^[{Q-Pf 



we can rewrite fi^ac' ^'i- (B-3), as: 



-ml]+-{P'-m') + -{Q'-m' 



in terms of the following integrals: 

d^Pd^Q 



I^ 

-'2 



1 



(27r)8 _ ^2) (p2 _ ^2) [(g _ p)2 _ ^2j ' 

d'^Pd^Q 1 

(27r)8 _ ^2) (p2 _ ^2) ' 

rf^prf^g 1 



(27r)8 (g2 - m2) [(g - P)2 - ml] 



(B.13) 
(B.14) 

(B.15) 

(B.16) 
(B.17) 



Defining 



J(a, 6) 



d^Pd^Q 



(27r)8 (P2_a)(g2_5) 



(B.18) 



we have = J{m ,m ) and /| = J{m ,ml). The dimensional regularization of J{a,b) is 
straightforward, yielding: 

where d = 4 — e. 



1 /e^A2\^ 



(47r)^ V 47r 



-| 2 



(B.19) 
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On the other hand, the evaluation of the integral If is extremely involved, essentially due 
to the absence of factorization of terms containing different masses. This calculation was 
performed in Ref. [199] and the result is 



b,REG 



An J (2 



TT 



\2d 



4(l + 2^) + -[42 log(42)] 



-2z [log{Az)f + 2(1- z)^z)+0{e)^ 



(B.20) 



With z = ^, 



A 



/3i 



1 + e A + + 0(e^ 



[r(l + e/2)f 
(l-e/2)(l-6) 
3 

2"^ 

7 3 1 2 TT^ 

7 + -7 H 

4 2 ' 2 ^ 24 



Az { [2-log(42;)] 2i^i( 1,1,2; ' 



d 

— 2F1 ( a, 1, -; z 



^ 2^,(1, l,c; z) 



(B.21) 
(B.22) 
(B.23) 



a=l 



(B.24) 



and 2F1 is the hypergeometric function, defined by: 



{a)k{b)k z^ 



k=0 



ic)k k\ 



(B.25) 



where {a)k = ^[."^^'^ is the Pochhammer symbol. In Ref. [199], simplified expressions for 
valid in the regions 2; > 1 or z < 1 were also derived. 
Taking the results in Eq. (B.19) and in Eq. (B.20) into the Eq. (B.14) and expanding 
around e = 0, one obtains, after a long algebra: 



Qexc 



2 4 

g m 



2 647r4 



poles in e = 



1 

+ 2 



7 + log ^ 



V3 



rrij 



Am? 



Am? ^ 

+ 0(6) 



7 + log ( 



m" 



V2 



mz 



Am^ 



(B.26) 
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where 



2(70 - 4/32) - 8^ (270 + /^a) + 16/32 + 4(1 - zf ^z) + 

+ \og{4z) {8z [2(l-z) /3i + 7i]} + [log(4^)]' {4^2-6^} 
2(7i - 4/3i) - 82(271 + /3i) + 16/3i + 8^(3 - 2z) hg{Az) , 



-6 - 24 2 + 16 r 



(B.27) 
(B.28) 
(B.29) 



with 



3 7(7 - 2) 
70 = -T+ ^ + 



7i 



4 
1 - 



24 



7' 



(B.30) 
(B.31) 



Finally, the 2-loop contribution to the vacuum expectation value counterterm, X, is then 
defined in order to cancel the poles in Eqs. (B.ll), (B.12) and (B.26). 

Therefore, collecting the finite parts of Eqs. (B.ll), (B.12) and (B.26), we obtain the 
final expression for the renormalized exchange vacuum contribution to the thermodynamic 
potential: 



exc,REN 



m 



viiz) + 



7 + log — 



V2{Z) + 
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[2 a(m' 
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Appendix C 



Vacuum and in-medium direct 
parts of the two-loop 



The third term in the diagrammatic expansion of the perturbative thermodynamic po- 
tential, shown in Fig 6.1, is the direct term. In this appendix, we concentrate on the exphcit 
evaluation of both vacuum and in-medium parts of this contribution. As usual, the full 
renormalized form of the direct term of the thermodynamic potential is obtained through 
the addition of the appropriate counterterms, as shown in Fig. C.l. 

Figure C.l: Renormalized direct contribution to the thermodynamic potential, represented 
by in-medium diagrams. As before, the crosses indicate (vacuum) counterterm vertices. 

-7-REN ^ _ <^ 

-^tadpole — 

Figure C.2: Definition of the renormalized tadpole integral. 
Defining the renormalized tadpole integral, as in Fig. C.2, we can rewrite the renormal- 
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ized direct contribution to the thermodynamic potential as: 



dir.REN 
Y 



1 (Npf 



PV 2 



/t-REN \ ' 
V tadpole / 



2ml 



/t-REN \ ' 

V tadpole / 



(C.l) 



Following the MS prescription, the counterterm is defined through the renormalization of 
the 1-loop fermionic self-energy, cancelling exactly the pole in e = of the tadpole integral, 
yielding: 

REN 

T-REN ^ 

tadpole 
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f — m 





-g) 4m / S^^{-V 



REN 



+ (—(?) 4m 



(27r)3 2Ep 
rf3p 1 
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(27r)3 2Ep 

d 

tadpole -'tadpole ; 



_exp [Ep/T - /i/T] + 1 exp [E^/T + /i/T] + 1 



7-vac,REN I 7-med 



(C.2) 



with = + m^. 
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tadpole 



-g 2m 



d^p 1 
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_exp[Ep/T - /i/T] + 1 exp[Ep/T + ^/T] + 1_ 



,(C.3) 



and 
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(C.4) 



Finally, taking these results into Eq. (C.l), we obtain the following expression for the 
renormalized direct term of the thermodynamic potential: 

\2 „ ( KT ^2 



dir,REN 



2ml 



/t-REN \ ' 
V tadpole^ 



I 7-vac,REN I rmcd 
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where 
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exp[E, - /i/T] + 1 exp[E, + /i/T] + l_ 



(C.9) 



in terms of the dimensionless quantities E\^ ^ m^ /T"^ and z = p/T. 



Appendix D 



QUASIPARTICLES IN HOT, ISOSPIN-DENSE 
MEDIA: MODIFIED PION PROPAGATORS 
AND PHYSICAL MASSES 



In what follows we sketch the calculation of the physical masses in the dressed theory in 
the regime in which fij < m.^ following Refs. [31]. The physical masses in the regime with 
fJ'i ^ iTT'TT can be computed in a similar fashion, as was also discussed in Refs. [31], as well 
as the one in which fij ^ m,r- 

The tree-level physical masses ^: 

which are associated with the following propagators (omitting the thermal piece) 

Ddj{p)oo = — , Ddj{p)+- = ^ 7^ , (D.2) 

where = 6^^. On the other hand, the dressed propagators (denoted by D) are given in 

^Here we adopt the notation for the vacuum pion mass in the isospin symmetric case. In Refs. [31], 
this mass is denoted by m. 
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terms of the self-energies as 

DDjip)oo 

whereas the self-energies up to one loop can be written as 



i 



(D.3) 
(D.4) 



00 



am. 
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where 



0^00 



-k- 327rhlJr-2Io + 4I 



e(/i/)8m^ J , 
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(D.7) 
(D.8) 
(D.9) 



with a = (m^/47r/^)2 , = 2.9 ± 2.4, h = 0.005 (cf. Table B.l of the last Ref. in [31], i.e. 

Ph.D. Thesis of C. Villavicencio), and the following additional convenient definitions: 

mu — md 6m 
m ' 



Jo = /o(m^,T) 
/ = I{m^,T, fii) 
J = J{m„,T,ni) 



mu + md 

l/^/l ' 
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dx — 1 2nB(m^x) 



(D.IO) 
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dx y . 
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^^(m^x - + ^^(m^x + , (D.13) 



dx xa/x^— T 



Therefore, we can rewrite the dressed propagators as: 

i i 
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where we have defined the dressed masses and the parameter h: 



mo 



acToo = 



m± = m, 



1 + 2"^oo + 0{a^) 



1 + -aai + O(a^) 



h = fij — -acTo + O(a^) 



(D.18) 
(D.19) 
(D.20) 



the latter playing the role of a dressed contribution to the zero-component of the 4-momentum. 
From Eqs. (D.15) and (D.17), it is clear that the dressed propagators have the same struc- 
ture as the original ones, Eq. (D.2). Therefore, the physical masses in the dressed theory 
are given by 



1 + 2*^^00 



1 + -aai ± -a — 
2 2 



(D.21) 
(D.22) 



up to 0{a) and with 



o"oo 



k- 327ihlj7-2Io + 4I 



Ox 



--h + 2/o 



(D.23) 
(D.24) 
(D.25) 



which reproduces the results of Refs. [31]. In our effective theory, we use the dressed masses, 
m± and mg, and the function h for the pion quasiparticles. 

Following a similar procedure for the second phase for ^ vfi-^^ the equations remain 
the same as a first approach, but replacing by 



Appendix E 



Some typical integrals and traces 
in a magnetic background 



In this appendix, we derive and collect some results for integrals and Dirac traces ap- 
pearing in the analysis of QCD thermodynamics in a magnetic field. 

• Properties of Dirac gamma matrices and the projector Vq = {1. + «7^7^)/2 over the 
LLL degrees of freedom: 







Vo 











^J.=o,3 




Tr [-f^-f^YY] 


u,p=0,3 





Tr [Y^Y] 


u=0,3 





Tr [YY] 




0, 



so that, for any given matrix M. that is a linear combination of 1,7^ and 7^7*^, with 
/i, z/ 7^ 1, 2, we have: 

Tr[7S^-M] = (E.2) 
Tr po-M] = ^Tr [M] . (E.3) 
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From these properties, one concludes that the LLL traces appearing in Chapter 5 (e.g. 
in Eq.(5.55)) are essentially occuring in a (l+l)-dimensional subspace (/i = 0, 3). This 
reduction shall be denoted by Trj^g- 

Integrals in the phase (l)x,y, defined in Eq. (5.27): 
Separating each of the 8 coordinate integrations, we have: 



exp 



Kpo - qo)ixo - yo) + i(p2 - q2)ix2 - y2) + iiPs - q3){x3 - ys) 



yi 



Xi 



eB ( eB 
— exp< — 

TT I 

eB I 

— Ji{k[)Ji{ki) Jx{kQ,po,qo) Jx{k2,P2,q2) Jx{k3,P3,q3) 



P2_Y 

eB) 



eB) 



yi 



eBJ 



TT 



Jy{ko,Po,qo) Jyik2,P2,q2) Jyiks^Ps^qs) 



where the integrals over Xi and yi are of the form 



(E.4) 




e5^"Pl AeB 



[k^ + 2i{p2 + q2)k + {p2 - q2y] 



(E.5) 
(E.6) 



while the integrals over Xj and y^ for i = 0, 2, 3 give directly momentum conservation 
delta functions in these directions and are respectively of the form: 



/oo 
dx = 27r5{k-p + q), 

-oo 

/oo 
dx e^('=+^'-«)^ = 27r 6{k + p - q) . 
-oo 

Thus, the phase becomes: 



(E.7) 
(E.8) 



= i27rfS^'^'''\k'-p + q)i27r)'S^'''^''\k + p-q) 

exp{ - ^ [k'l + kl + 2i{p2 + q2){k[ + k^) + 2{p2 - ga)'] } , (E.9) 



where we have defined the 3-dimensional delta functions, excluding the i = 1 direction. 
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Integral I{ki,k2,eB), defined in Eq. (5.61^ 



I{h,h,eB) 



dxidy 



lg«fci(j/i-a:i 



eB 



X exp 



2 ^\ eBJ 



TT ./ (27r 

2 

' Xi 



-{27r)d{p2 -q2- h) 



eB 



P2 V 



yi 



(E.IO) 



We can further simplify this expression as follows. Using (def.: x = {xi + X2)/2 and 
A = xi - X2) 

1 2 



P2_\' 

eBJ 



dxidyi 



P2 _ xi + yi 
eB 



{xi+yif ,2,2 
^ + xi+yi 



P2_ 

eB 



— X 



2 A2 



A 2 


eB 








— = 




[(- 


2 


2" 






eB. 



dxdA , 



(E.ll) 

(E.12) 
(E.13) 



we get: 
I{k^,k2,eB) 



dx eB , , . 
I dAe 

Li TT 



eB 
2 



ik^ \2 , fef 



dp2dq2 

w 



{2ti)5{p2 -q2-k2 



X exp 



\eB 



(E.14) 



dxeB ^j4. [2tt f dp2 ( r^y f P2 

e 2es w — / exp< — eis — 

Li TT \l eB J 27r ^\eB 



P2 - fcs 

e5 



(E.15) 



where we have displaced P2 P2 + C-Bx and q'2 >■ q'2 + before using the S function 
to evaluate the q2 integral. Completing the square in the final integration, over p2, one 
obtains \pl + (ps - ^2)^ = 2(p2 - k2/2y + k^/2]: 



I{ki,k2,eB) 



dxeB /27r 1 eBn 



27r 



-e . 



(E.16) 
(E.17) 



Appendix F 



Hot pure gauge QCD pressure 
at two-loops 



In this appendix, we collect results for diagrams contained in the gluonic part of the 
thermodynamic potential, ^'qcdi defined in Eq. (5.44). 

These contributions do not contain fermions, coming directly from hot Yang-Mills theory, 
and were computed long ago (see, e.g., [150, 200, 149]). Picking the renormalized results for 
each piece from the appendix of Ref. [200], we have 




A-Z 



, X A-Z 




— ^ — 



(F.l) 
(F.2) 
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A-Z 



A-Z 



1 1 1 

2 Wq 




{N'i - l)N, g^T' 



9 



4 X 144 



(iV^ - l)iV, /T^ 



144 



(F.4) 



where the superscript A-Z labels the diagrammatic representation adopted in Ref. [200]. 
Summing all the terms, we arrive at 



1 



90 ^ " ^ ^ 144 



(F.5) 



Appendix G 



The EUCLIDEAN ACTION FOR 
ISOSPIN-DENSE CHIRAL MODELS 



In this appendix the derivation of the well-known result for the euclidean action in the 
presence of a nonzero conserved charge involving momenta canonically conjugated to the 
fields is sketched, following the first Ref. in [149]. 

In the presence of a chemical potential, the partition function includes the constraint of 
charge conservation explicitly: 



[DU][D<^] exp{- / d'^x 



d 

-Hz— $ + n] - /in[$, n] 

or 

where 11 is the momentum conjugated to the field $, being defined by (r = it): 



(G.i: 



n = = (c 21 

Note that the presence of a conserved charge that depends explicitly on the momentum 11 
will modify the result of the functional integration over 11, rendering a modified euclidean 
action. 

The isospin number ri/, Eq. (7.2), indeed depends on the momentum conjugated to the 
pion fields (let $ = ($0, $1, $2, $3) = (c^, v^i, vra, TTg)): 
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where we have used a Lagrangian of the form : 

'(^ = - V $o) , (G.4) 

which is compatible with chiral models, since $j<I>j has the form of a scalar under standard 
chiral rotations and the $o (c") dependence of the potential represents a possible explicit 
breaking, commonly included. For example, in the LSM (cf. Appendix A neglecting con- 
stituent quarks) one has: 

y($,$„$o) = V{[a' + n%a)=^{[a' + n']-v'y-ha. (G.5) 
The isospin number written in terms of the momenta Ilj reads: 

m = $i9°$2 - $ia°$2 = $in2 - $2ni . (G.6) 

The exponent in the integrand of the partition function, Eq. (G.l), becomes then: 

- y"rf'^x[-n,9o<f. + ?/[$,n]-^i,(<i>in2-<f2ni)], (g.7) 



with 



7/[$,n] := n,9o<f. - c ^='^ Jn^ + l(v$,)(v<i'.) + y($.$.,0o) • (g.s) 



Although modified, the functional integration over the momenta Ilj remains gaussian and 
can be done exactly once the terms linear in Ilj are absorbed via completing the squares in 
the exponent above. Finally, the result for the partition function reads: 



e"^^ , (G.9) 
with the euclidean action for the isospin dense system being given by: 

+V($,$„<l>o)|. (G.IO) 
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For the particular case of the LSM, we have $ = ($0) •^•2, ^3) = (c", tti, 712, tts) and 
Se = j (fx\^{d^a){d^(j) + ^(9^7f)(9^7r) + /i/fvriSoTTa - vra^oTTi] - ^/ij (tti + Til) + 

+r([a2 + 7r2],a)|. (G.li: 
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